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RESEARCH AND DEVELOPMENT FOR A 
GROUND-BASED HYDROGEN-MASER SYSTEM 


Final Report 
Contract NSR 09-015-098 


1. INFORMATION AND BACKGROUND 

On July 1, 1969, development of the atomic hydrogen maser conducted by Hewlett 
Packard (HP) under Contract NASW 1337 was transferred to the Smithsonian Astrophysical 
Observatory (SAO), Cambridge, Massachusetts, where it continued under Contract 
NSR 09-015-098. The objective was to continue the development of hydrogen masers 
and clock systems for terrestrial and satellite use and to perform research on the wall 
coatings of the maser storage bulb. 

The work statement was divided into the following four tasks, which are discussed 
in subsequent sections: 

I. Establish a laboratory to continue hydrogen-maser research and development. 

H. Set up and test the existing autotuner. 

HI. Evaluate the PTFE -coated maser storage bulb and test it for wall-relaxation 
rate. 

IV. Assemble engineering model of the spacecraft maser and perform limited 
mechanical testing on it. 
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2. TASK I 


Six maser oscillators were transferred to Contract NSR 09-015-098 with SAO on 
July 1, 1969: 


NASA-1 
NASA -2 
H-10 #6 
SR-1 MK n 
SM-1 
ONR-1 


Originally built under Contract NAS 8-2604 
Originally built under Contract NAS 8-2604 
Originally built under Contract NAS 8-2604 
Spacecraft maser developed under Contract NASW 1337 
Small-maser system developed under Contract NAS 8-2604 
Original prototype maser built in 1960 


NASA-1 has been modified under Contract NAS 5-11598 to permit adiabatic fast- 
passage-state experiments. It has also been equipped with a storage bulb carefully 
evaluated by the Lyman Laboratories at Harvard University. It serves as the SAO 
lab reference. 


NASA-2 is a source of spare parts and electronics to keep NASA-1 going. 

H-10 #6 is used in experiments on wall-shift and wall -relaxation rates. 

SR-1 MK n is under continuing development (see Task IV). 

SM-1 has been extensively modified at SAO; it serves as a traveling frequency 
standard. 


ONR-1 has been lent to Professor Ramsey at the Lyman Laboratories for flexible- 
bulb experiments. 

Other equipment and supplies were transferred to SAO as part of the contract trans- 
fer. New benches, cabinets, hand tools, and wiring supplies were rapidly obtained, and 
the lab became operational by the middle of July. As plans developed, major capital items 
were purchased or obtained by transfer from other labs at SAO and from NASA/Marshall 
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Space Flight Center (MSFC); these items included a bulb-curing oven capable of tem- 
peratures to 750° C, a Leeds and Northrup millipoint strip-chart recorder, and 
amplifiers. 

2 2 

The original lab space of 641 ft was expanded by 434 ft in an adjacent room 
when the oven and its controlling equipment were obtained. 

Requirements for electronic tests have been very nearly met by equipment from 
NASA/MSFC. Except for a vector voltmeter, an L-band spectrum analyzer, and a 
low-frequency spectrum analyzer, the requirements for present plans have been met. 

The SAO/Harvard College Observatory (HCO) model shop located in the same 
building has very competently provided services. This shop is well run and has personnel 
familiar with vacuum requirements and atomic clocks. Glassblowing facilities continue 
to be provided by the same outside vendor in Cambridge and by the Massachusetts 
Institute of Technology (MIT) glass shop. High-vacuum welding facilities are available 
at the nearby Cambridge Electron Accelerator Laboratory. 

Because electronic components are readily available from nearby supply houses, 
it is unnecessary for us to maintain an inventory. The purchasing department at SAO 
has quickly procured small quantities of electronic components. 

Finally, our maser lab and others at SAO, Harvard, and MIT have lent each other 
equipment very infrequently used, such as a residual-gas analyzer for vacuum studies. 
This cooperation does much to alleviate budget constraints and often serves to broaden 
the capabilities and outlook of the people involved. 

In contrast to the Varian-HP labs where the maser work was previously carried 
on, the present setup offers substantially improved capability for research, develop- 
ment, prototype assembly, and testing. 
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3. TASK II 


The existing autotuner, which uses binary logic, has been tested with the SAO 
lab maser as reference. Its performance has been up to expectation. To test the 
autotuner-maser combination, we decided to include it with SM-1, which was being 
compared with the clock ensemble at the National Bureau of Standards (NBS), Boulder, 
Colorado, to obtain a more definitive hydrogen-cesium frequency relationship. To 
make the device more compact, the time sequencer was rebuilt. 

To the time of writing this report, the performance of the autotuning device com- 

4 5 

bined with the maser has been excellent. Long-term stability (10 to 10 sec) has been 

-14 

measured at 1. 5 X 10 in fractional frequency stability. 

The results of the experimental redetermination of the hydrogen-cesium hyperfine 
transition frequency are given in the following paper by Hellwig, Vessot, Levine, 
Zitzewitz, Allan, and Glaze. These results are related to the wall-shift data taken in 
Task III . 

The performance of the autotuning system is described in the following NASA 
Tech Brief 70-10616. 



.- * i» 
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4. TASK HI 


4. 1 Introduction 


The effect of wall collisions on the phase of the hydrogen atom has long been the 
cause of inaccuracy in the output frequency of the maser. Further, the relaxation 
rate of the atom as a result of such collisions, whether due to recombination into 
molecules (or chemical loss) or to dephasing with the oscillating r.f. field, has been 
a limitation on the linewidth of the stored aggregate of atoms. These wall processes 
limit the accuracy of the maser and affect its stability. 

The properties of polytetrafluoroethylene (PTFE) films are being studied to 
determine whether reduced relaxation rates are obtainable and whether it is practical 
to operate the maser at a temperature where the wall shift can be made small. The 
H-10 #6 maser has been modified to allow more experimental flexibility and to permit 
oscillation at temperatures in excess of 100°C. 

An important aspect of this work relating to the wall-shift measurements is know- 
ing the unperturbed frequency of the hyperfine separation of hydrogen. This is usually 
determined by measuring a series of bulbs that have widely differing diameters but are 
coated in the same manner. Assuming that the wall-collision rate (hence the average 
phase shift per bounce) is proportional to the reciprocal of the diameter of the bulb, 
one arrives at a plot whose extrapolation to infinite bulb diameter (zero collision rate) 
gives the desired hydrogen frequency. The difficulties involved in this technique are 
due to the variability of the surface texture and the consequent error in the collision 
rate. The experiments give this value; however, much greater weight to this deter- 
mination is possible by comparisons with other laboratories where different types of 
teflon and different coating techniques have been used. . 

As mentioned in Task n, the small maser has been compared with the clocks of 
Harvard and NBS to establish the hydrogen frequency and its relation to the best 
available embodiment of the cesium frequency. 
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4.2 Wall-Shift Measurements 


Wall- shift data as a function of temperature for the 4" and 7" bulbs have been 
obtained. These bulbs were coated with tetrafluoroethylene (TFE) teflon and cured in 
a rapid cooling cycle. The data show that the temperature at zero wall shift is much 
lower than anticipated. The temperature has been determined to be 83° ± 1°C. The 
4” bulb gave the same crossover temperature (within 1°C) as did the 7" bulb. With 
these two bulb diameters, we can extrapolate the wall shift to zero collision rate 
(infinite bulb diameter). The extrapolation of the frequency- shift plot as a function of 
(diameter) gives a value for the unperturbed atomic hydrogen hyperfine frequency 

_3 

within 2 X 10 Hz of the determination made with the clocks of Harvard and NBS and 
the SM-1. This new value lies within the error limits of the previous experiment. 

The wall-shift data have been reduced in terms of phase shift per wall collision. 

The reduced data points, when plotted against temperature, lie on a straight line 

-7 

with slope 1. 31 X 10 rad/°C. The intercept of zero phase s hi ft per collision with 
temperature is, of course, the same as before. It is very interesting to report that 
the phase shift for collision data for both the 4" and the 7" bulbs lies well within 
1X10 ° rad over the entire temperature range. This confirms that, at least for 
these two bulbs, the coatings are similar. 

Having successfully obtained the data on the 7" and 4" bulbs coated with PTFE 
under quenched conditions, we annealed the 4" bulb. This process changes the crystal 
structure of the PTFE coating and should affect the wall shift and reduce the wall- 
relaxation rate. 

The PTFE we have been using has a very high molecular weight, in the order of 
10 or greater. This is desirable, since with large molecules it is possible to pre- 
sent a surface with fewer nonfluorine end groups than smaller molecules. There is, 
however, a difficulty in annealing the polymer of high molecular weight; it is extremely 
viscous and requires a very long time to cool through the crystallizing range, 340 to 
290°C. The rate should be between 0. 1 and 0. 05°C/min. 
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A special oven controller with a clock-driven potentiometer was built. This 
controller was used in place of the set -point adjustment. We provided for a 
0. 05°C/min rate; roughly 18 hr was required to cool from 340 to 290°C. 

The 4" bulb was annealed in this 18-hr cycle. It was then placed in the maser, 
and wall-shift measurements were made. We found that the bulb showed entirely 
different properties in its annealed state. The data were discussed in "Studies of 
Hydrogen Maser Wall Shift for High Molecular Weight Polytetrafluoroethylene, " a 
paper read by Vessot and Levine at the 24th Annual Frequency Control Symposium. 

A copy of this paper follows. 

This work led us to the invention of a method for eliminating the wall s hif t in the 
hydrogen maser. The principles of the system are discussed in a paper published in 
Metrologia and reproduced at the end of this section. 

A method for making an automatic temperature-control system to ensure that 
the bulb is maintained at the temperature of zero wall shift was disclosed to NASA. 
Subsequently, this disclosure was filed by NASA on behalf of the inventors. The 
NASA case number for this invention is HQN-10654. A copy of the original disclosure 
is appended to illustrate the principles of the automatic system. 
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Measurement of the Unperturbed Hydrogen 
Hyperfine Transition Frequency 

HELMUT HELLWIG, member, ieee, ROBERT F. C. VESSOT, MARTIN W. LEVINE, 

PAUL W. ZITZEWITZ, DAVID W. ALLAN, and DAVID J. GLAZE 


Abstract — The results of a joint experiment aimed primarily at the 
determination • of the frequency of the H 1 hyperfine transition 
(F = 1, m? =0) *-* (JF = 0, m F = 0) is reported. In terms of the 
frequency of the Cs 133 hyperfine transition (F = 4, m F — 0) «-> (F = 
3, m F = 0), defined as 9192 631 770 Hz, forthe unperturbed hydrogen 
transition frequency the value 

v u = 1420 405 751.768 Hz '; >• -* 

is obtained. This result is the mean of two independent evaluations 
against the same cesium reference, which differ by 2 X 10~* Hz. 
We estimate the one-sigma uncertainty of the value vu also to be 
2 X 10 -1 Hz. One evaluation is based on wall-shift experiments at 
Harvard University; the other is a result of a new wall-shift meas- 
urement using many storage bulbs of different sizes at the National 
Bureau of Standards. The experimental procedures and the applied 
corrections are described. Results for the wall shift and for the fre- 
quency of hydrogen are compared with previously published values, 
and error limits of the experiments are discussed. 

I. Introduction 

A SERIES of measurements involving hydrogen 
masers, cesium beam frequency standards, and 
the NBS clock system was made at the National 
Bureau of Standards (NBS) in Boulder, Colo., during the 
period from November 1969 to February 1970. The 
original purpose of these experiments was an evaluation 
of the stability performance of the involved systems and 
devices. A report of these results will be made elsewhere [1]. 

A unique opportunity thereby existed for a remeasure- 
ment of the frequency of the unperturbed hydrogen (H 1 ) 
hyperfine transition (F = 1, m, = 0) «-» (F = 0, m F = 0) 
in terms of the cesium (Cs 133 ) hyperfine transition (F = 4, 
m F — 0) <-» (F = 3, m F = 0), defined as 9192 631 770 Hz. 
This paper reports the results of our frequency measure- 
ments. 

The hydrogen maser frequency differs from the un- 
perturbed atomic transition frequency because of a variety 
of effects including cavity pulling, spin exchange, magnetic 
fields, second-order Doppler (related to the temperature of 
the storage vessel), and collisions with the walls of the 

Manuscript received June 5, 1970; revised July 13, 1970. This 
paper was presented at the 1970 Conference on Precision Electro- 
magnetic Measurements, Boulder, Colo., and the experiments at 
Harvard University were supported by NASA. 

H. Hellwig, D. W. Allan, and D. J. Glaze are with the Atomic 
Frequency and Time Standards Section, National Bureau of Stand- 
ards, Boulder, Colo. 80302. 

R. F. C. Vessot and M. W. Levine are with the Smithsonian 
Astrophysical Observatory, Cambridge, Mass. 02138. 

P. W. Zitzewitz was with the Lyman Laboratory of Physics, 
Harvard University, Cambridge, Mass. 02138. He is now with the 
Department of Physics, University of Western Ontario, London, 
Ont., Canada. 


storage vessel. A detailed discussion of these effects is 
found in [2]. They can be measured and accounted for 
with fractional uncertainties of less than 10 13 with the 
exception of the wall-collision effect (wall shift). Thus, the 
measurement of the unperturbed hydrogen transition 
frequency involves primarily a measurement of the wall 
shift. We performed two different and independent 1 
measurements, which are referred to hereafter as experi- 
ment 1 and experiment 2. 

Experiment 1 relies on a previous determination of the 
wall shift at Harvard University [3]: as a result of this 
wall-shift evaluation the Harvard reference maser has a 
known wall shift. A comparison between the frequency 
of this maser and a cesium primary frequency standard 2 
(NBS-III), using a different maser as a transfer standard, 
gave a value for the unperturbed hydrogen hyperfine 
transition frequency. Experiment 2 is a separate and new 
measurement of the wall-shift correction with, however, 
an important difference from any previous wall-shift 
measurement: the reference oscillator that was used in 
the measurement of the frequency <?f the maser equipped 
with storage bulbs of different sizes was directly related 
to the frequency of a cesium primary frequency standard 
(NBS-III). Wall-shift correction and hydrogen-cesium 
comparison in experiment 2 were thus not separate steps 
but were integral parts of the same measurement. 

The published values of the unperturbed hydrogen 
hyperfine transition frequency disagree among themselves 
considerably more than is expected from the published 
values of the accuracy. In addition, there has been an 
almost exclusive usage of one particular wall-shift correc- 
tion in the various publications of the hydrogen frequency. 
(See Tables I and II. ) These were further motivations for 
our measurements. 

The value for the hydrogen hyperfine transition fre- 
quency may be written as 

v H = (1420 405 751 + v£) Hz (1) 

and we shall discuss in the following only v{,. 

The wall-shift correction may be written as approxi- 
mately 

tow « (K/D)[ 1 + ai (T - r,)] (2) 

1 They are not completely independent in the sense that a common 
cesium reference (NBS-III) was used. 

* For our purposes a cesium primary frequency standard is an 
apparatus that has been experimentally evaluated in relation to 
all known perturbing effects [4J. 
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TABLE I 

Wall-Shift Results 


Author and Reference 

Approximate 
Year of Measure- 
ment 

Number of 
Bulbs 

Teflon 

Arw-D 
(mHz -cm) 

(A" 1 ) 

Crampton et al. [5] 

1963 

2 

TFE 

852-201 

470 ± 47 

not 

Vanier and Vessot [6] 

1963-1964 

4 

Du Pont 
FEP-120 

at 35 °C 
528 ± 5 

measured 
-(5 ± 1) X 10-* 

Mathur et al. [7] 

1964 

4 

Du Pont 
1962 

FEP-120 

at 40 °C 
513 ± 76 

not 

Elkina et al. [8] 

1968 

2 

Blend 33 

53.9 percent solids 
1964 
Du Pont 
Teflon D 

at 35°C 
1993 ± 75 

measured 

not 

Zitzewitz et al. [3] 

1968-1969 

18 

FF 

(tempera- 
ture not 
given) 

386 ± 8 

measured 

-(12 ± 1) X 10-* 

Menoud and Racine [9], [10] 

1968-1969 

2 

B1 

5'* u Aids 

. U 

A .1 

FEP-120 

at'40°C 
515 ± 28 

not 

Hellwig et al. [this paper] 

1969-1970 

11 

Du Pont 
TFE 

at 40°G 
528 ± 17 

measured 

not 


1 

Blend 42 

34.3 percent solids 
Du Pont 
1969 

at 24°C 

measured 


where 7\ is the reference temperatu: usually chosen at 

40°C, T is the operating temperature of the storage bulb, 
«. is the wall-shift temperature coefficient, D is the 
lame ter of the storage bulb, and K is the wall-shift 
coefficient. 

To our knowledge only five independent wall-shift 
measurements have been published in the past; one 
reason for this must be sought in the rather tedious and 
time-consuming experimental effort necessary. The six 
published results are summarized in Table I together 
with the results of this paper. Table I gives the author(s) 
and bibliographic reference, the year(s) of the actual 
experiment, the number of bulbs different in coating 
and/or size, the type and year of purchase of the Teflon, 3 
and the values A v w -D and a, of (2) together with the 
accuracy claims. 

In Table II we list and compare values for v„. Given 
are the author(s) and published reference, the year(s) of 
the actual measurement, the Teflon that was used in the 
storage bulb of the hydrogen maser, the applied wall- 
shift correction listed by the corresponding author from 
Table I, the hydrogen maser type that was used, the 
cesium reference standard, and the value for v{,. The 
values for v/, are rounded in the last digit. 

A discussion of Tables I and II will be postponed until 
the end of this paper where we will take a critical look at 
them in connection with a discussion of our new values. 

* Trade names and the names of manufacturers are used in this, 
paper for the sole purpose of conveying scientific and technical 
'■(.formation, and their citation is not to be construed as an endorse- 
ment or approval of commercial products or services by the authors’ 
organizations. 


II. Measurement Procedures 

Experiment 1 involved a transportation of the portable 
hydrogen maser of the Smithsonian Astrophysical Ob- 
servatory (SAO) first to the Lyman Laboratory of Physics 
at Harvard University (HU) and then to the Atomic Fre- 
quency and Time Standards Section of the National 
Bureau of Standards (NBS). The SAO maser is described 
in detail in [20]. At HU its frequency was compared with 
the HU reference hydrogen maser. The wall shift of the 
HU maser was known from a previous evaluation [3]. 
Both masers were timed and the appropriate bias correc- 
tions (see Section III) were applied. From the thus 
corrected difference frequency of both masers the wall- 
shift correction of the SAO maser could be inferred as 
Aj >w = +20.19 mHz. This experiment was done in De- 
cember 1969. The SAO maser was brought to NBS, 
where its frequency was referenced to cesium on January 
22 and 23, 1970. 

Experiment 2 involved the NBS experimental hydrogen 
maser NBS-H2. Its general configuration is similar to 
tHe'Varian H-10 hydrogen masers. During the period from 
November 1969 to January 1970 this maser was suc- 
cessively equipped with bulbs of different sizes ranging 
from 7.5 to 20 cm in diameter. A total of 15 frequency 
determinations was made of which 11 were used for 
the final analysis. Each individual measurement was 
made with a reference frequency source whose frequency 
was known in terms of cesium (NBS-III) [19]. Thus 
the measurement of the wall shift was' simultaneously 
the measurement of the hydrogen-cesium ratio, and the 
value for zero wall shift was obtained by extrapolating 
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TABLE II 

Hydrogen-Cesium Measurements 
v R = 1420 405 751 Hz + va 


Author and Reference 

Approximate 
Year of Meas- 
urement 

Teflon Used 

Wall Shift 
From Reference 
(See Table I) 

Hydrogen Maser 

Cesium Reference 

*h 

(mHz) 

Crampton et al. [5] 

1963 

<— Rela 

TFE 
852-201 
Du Pont 

ted — > 

Crampton 
et al. 
(1963) 

Harvard 

laboratory 

type 

NC 2001 

800 ± 28 

Vanier el al. [11] 

1963 

FEP 
Du Pont 
Lot 10144 

Vanier-Vessot 

(1963-1964) 

Varian 

laboratory 

type 

several Cs 
clocks 

via Loran C 

827 ± 20 

Peters el al. [12] 

1964 

FEP 
Du Pont 

Vanier-Vessot 

(1963-1964) 

Varian H-10 

HP 5060 A 

778 ± 16 

Peters and 
Kartaschoff [13] 

1964 

FEP 
Du Pont 

Vanier-Vessot 

(1963-1964) 

Varian H-10 

LSRH 

laboratory type 

785 ± 16 

Johnson and 
McGunigal [14] 

1965 

FEP 
Du Pont 

Vanier-Vessot 

(1963-1964) 

Varian H-10 

HP 5060 A 

781 ± 16 

Vessot el al. [15] 

1965 

FEP 
Du Pont 

Vanier-Vessot 

(1963-1964) 

Varian H-10 

NBS-III 

786 ± 2 

Becker and 
Fischer [16] 

1966 

TFE 
852-201 
Du Pont 

Vanier-Vessot 

(1963-1964) 

PTB-H2 

laboratory 

type 

HP 5060 A 

756 ± 3 

Becker and 
Fischer [16] 

1967 

TFE 
852-201 
Du Pont 

Vanier-Vessot 

(1963-1964) 

PTB-H2 

laboratory 

type 

PTB laboratory- 
type C2 and 
HP 5061 A 

758 ± 2 

Mungall et al. [33] 

1968 

FEP-120 
Du Pont 
Drum 10020 

Vanier-Vessot 

(1963-1964) 

NRC 

laboratory 

type 

NRC laboratory 
long cesium 
beam standard 

778 

(an explicit ac- 
curacy claim was 
not made) 

Chi el al. [17] 

1968 

FEP-120 
Interchemical 
vanaflex (1967) 

Vanier-Vessot 

(1963-1964) 

NASA-NX-1 

laboratory 

type 

4 units HP 5061 A 
1 unit HP 5060 A 

777 =fc 3 

Menoud and 
Racine [9], [10] 

1969 

<— Relat 

FEP-120 
Du Pont 

;ed — > 

Menoud- 

Racine 

(1968-1969) 

LSRH-H2 
laboratory 
type * 

Oscillatom 3 
(Ebauches) 

778 ± 4 

Bangham [18] 

1969 

PTFE 

Vanier-Vessot 

(1963-1964) 

NPL 

laboratory 

type 

NPL 

laboratory 

type 

734 ± 25 

Hell wig et al, 
experiment 1 
[this paper] 

1970 

<— Relat 

FEP-120 
Blend B-107 
Du Pont (1967) 

;ed — * 

Zitzewitz 
et al. 

(1968-1969) 

portable 
SAO maser 
transferred 
from HU 

NBS-III 
via NBS clock 
system and 
NASA-NP3 H-maser 

769 ± 2 

Hell wig et al., 
experiment 2 
[this paper] 

1970 

<— Relat 

TFE 
Blend 42 
Du Pont 
(1969) 

,ed — > 
this 
paper 

(1969-1970) 

• NBS-H2 
laboratory 
type 

NBS-III 
via NBS clock 
system and 
NASA-NP3 H-maser 

767 ± 2 


the measured data to 1/D = 0 (D s bulb diameter). 

Fig. 1 shows the measurement system used in both 
experiments. The receiver for the output frequency v M 
of the hydrogen maser is basically a triple-superheterodyne 
system. The first intermediate frequency is in the MHz 
range, the second in the kHz range, and the third is 
adjusted by means of a synthesizer (output frequency 
»s) to a slow beat frequency v B of the order of 1 Hz. 
If the external reference frequency v R is chosen to be 
nominally 5 MHz we have at the output of the second 


stage an intermediate frequency of (v M — 284v«) and 
for the final beat frequency 

v B - O'M — 284v fl ) — v a . (3) 

By comparison with (1) we see that v s ss 405 kHz. The 
maser frequency v u in terms of the reference frequency 
v R can then be calculated from (3). If v R is given in terms 
of the cesium transition frequency the hydrogen-cesium 
ratio can be computed using the appropriate corrections. 
A measurement precision of a few parts in 10 13 was 
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V » {v — 284y v 
B M R S 


Fig. 1. Measurement system used for determining frequency 
ratio between two frequency standards. 

desired. To take v B directly from the output of a cesium 
standard would have meant averaging times of many 
hours to attain this precision. This would have been 
awkward in experiment 1 and prohibitive for experiment 
2 since the NBS-H2 maser had no temperature control, 
a limitation that necessitated a fast measurement tech- 
nique (30 minutes maximum duration). 

Fortunately we also had available the hydrogen maser 
standard NP3 of the National Aeronautics and Space 
Administration (NASA). Its design is described in [21]. 
The most important features of this maser for the present 
purpose are its stability and the provision of synthesized 
standard output frequencies. We used one of its outputs 
on nominally v R = 5 MHz. However, the NASA maser 
not only served as a highly stable frequency source 
allowing short- and long-term measurements with a 
precision of a few parts in 10 13 but also as a calibrated 
transfer standard. As indicated in Fig. 1, the frequency 
of the NASA maser was constantly monitored by the 
NBS clock ensemble, which consisted of six commercial 
cesium standards. The clock ensemble generates the 
atomic time scale AT (NBS) with the frequency f at <nbs); 
this frequency is set by calibration with the primary 
frequency standard. The most recent calibration had 
been made in May 1969 with the NBS-III cesium beam. 
A coordinated universal time scale is generated from 
AT (NBS) and is called UTC(NBS) with a frequency 
of fctcinbr) = ■'at(nbs) 4 " Av UT c- During the course 
of the measurements we had Aj> utc = —299.995 X 10 -1 ° 
vat(nbs)- The nominal output of the NASA hydrogen 
standard NP3 was synthesized to reflect a coordinated 
universal time scale (UTC) [ 17 ]. Fig. 2 shows a plot of 
the fractional frequency of NP3 referenced to r axc(NBB) 
(zero on the vertical scale) versus time. From Fig. 2 we 
can derive the correction Av N as the offset of NP3 from 
"utccnbh) evaluated at the nominal hydrogen frequency. 
Several additional corrections are necessary to relate the 
frequency of the hydrogen maser to the unperturbed 
transition frequency in hydrogen. They are summarized 
in Table III together with those already mentioned. 



Fig. 2. Fractional frequency offset of NP3 from itjtcinb si- Derived 
from a 4-day average of the time error of NP3 with respect to 
UTC(NBS): 6 clocks. Zero on. the vertical scale corresponds 

to VUTC(NB8>- 


TABLE III 

Summary of Bias Corrections 


Offset of putc(nbs) from pat(nbs) 

Acutc 

Offset of NASA-NP3 frequency from itjtcinbs) 

A«w 

Wall shift 

Ayjp 

Second-order Doppler 

Av d 

Thermalization 

A vt 

Magnetic field (Zeeman) 

A vz 

Spin-exchange collision 

Apse 

Cavity tuning 

A rc 


The wall-shift correction was discussed previously in 
Section I, and (2), a semiempirical equation, was given. 
For (2) it is assumed that the wall shift is proportional 
to the wall-collision rate of the hydrogen atoms, i.e., 
to the inverse diameter of a spherical storage bulb. 

The second-order Doppler effect in the hydrogen maser 
due to the thermal motion of the atoms necessitates a 
correction of [15] 

Avd — CiTs ( 4 ) 

where C\ = +1.9557 X 10 -4 Hz -A -1 . We assume here 
that thermal equilibrium between the kinetic energy 
of the hydrogen atoms and the walls of the storage bulb 
is established. T s is then the temperature of the storage 
bulb. 

If the kinetic energy of the atoms is not fully in equi- 
librium with the temperature of the storage bulb,, we 
will commit an error in using T s - We should then introduce 
a correction Av T , which relates to the difference between 
the bulb temperature and the effective temperature of 
the stored atoms. Preliminary theoretical calculations 
of the thermalization process indicate [22] that Av T 
should be sufficiently small to warrant Av T = 0 within 
our uncertainty limits. Some preliminary experimental 
results based on velocity selection in an atomic hydrogen 
beam support this assumption [23]. 
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A magnetic field H causes a Zeeman effect on the 
energy level structure of hydrogen, and the corresponding 
frequency correction is given by [15) 

A v z = CjP (5) 

where C 3 = —2750 H z-G~ 2 , and H 2 denotes the average 
of the squared magnetic field over the volume of the 
bulb. It is possible to measure the average magnetic 
field ft by inducing transitions among the magnetic 
sublevels of the F = 1 state and observing a change in 
the amplitude of the maser signal. This Zeeman frequency 
is given by [15] 

u = c 3 \n\ (6) 

where C 3 = +1.4 X 10 6 Hz G~ l . The error that on e is 
forced to commit by substituting ft 2 from (6) for H 2 
in (5) in order to obtain Av z is negligible because very 
low and homogeneous magnetic fields (|i/| < 1 mG) 
are typically used. 

The last two corrections from Table III, Av aE and 
Av c , must be discussed together. The effects of spin 
exchange in combination with cavity tuning have been 
treated extensively in the literature [24]-[27]. The net 
result is that spin-exchange frequency shift and cavity 
pulling have a similar functional relationship to the 
atomic-resonance linewidth. As a consequence, cavity 
ling and spin-exchange shifts cancel each other when 
laser is “tuned,” or 

(A Vc + Av 8E ) tun ,d = 0. (7) 

, “tuned” condition is here defined as a setting of the 
maser cavity such that the output frequency of the 
maser v M is independent of changes : the hydrogen 

pressure (beam intensity). 

In experiment 1 the SAO maser was tuned in separate 
runs, both manually and automatically using its own 
automatic cavity servo, until the output frequency was 
unaffected by changes in the beam flux. The actual 
tuning element was a varactor diode mounted' inside 
the cavity. Details of this procedure can be found in 
[20]. Averaging times of about 10-100 seconds were 
used to measure ve- 
in experiment 2 the beat frequency v B for the tuned 
condition was calculated from four measurements of 
the maser frequency at two cavity settings and two 
beam intensities using a linear interpolation. 4 The two 
cavity settings corresponded to maser frequencies that 
were typically a few parts in 10 12 above and below v M 
(tuned). In this procedure, the cavity settings did not 
have to be known in an absolute sense but had to be 
reproducible to Av M /v M 10“ I3 . This was assured in 
separate experiments. The actual tuning element was 
also a varactor diode mounted inside the cavity. One 
complete measurement of v B took typically 15 minutes 
at averaging times of about 10 seconds. 


From Table III and (3)— (7) we are now able to calculr.t 
the unperturbed hydrogen transition frequency referenced 
to NBS-III as 

Vg = 2&4vg + Vs + Vb + AvpTC 

+ Avg + Av,e + Av d + A V Z , (o) 

where 284 v B = 1420 000 000.0000 Hz with v B = 5 MtTz. 
All corrections Av are taken at the nominal hydrogen 
frequency. 

III. Experimental Results 

We will now discuss the experimental results that we 
obtained in experiments 1 and 2 together with the un- 
certainty contributions associated with each of the 
frequency corrections. We begin with experiment 1; 
its results are summarized in Table IV. 

The synthesizer was set at v a = 405 794.4200 Hz; 
no uncertainty is associated with this value. 

The period of the beat frequency was measured; the 
beat frequency v B is listed in Table IV. In order to obtain 
this value the maser had to be tuned. This can be done 
only with a certain precision. From the experimental 
results we estimate the associated uncertainty as ±3 X 
KT 13 . 

The fractional frequency offset of UTC(NBS) against 
AT(NBS) is nominally Av UTC = -299.995 X 10 _, ° 
vatinbs) as listed in Table IV; no uncertainty is involved. 

The actual synthesized output frequency of the NASA 
hydrogen maser standard deviated by a small amount 
Av N from the frequency v UTC(NBS) . For the date of our 
experiment we can l:e the corresponding bias correction 
directly from Fig. 2 and obtain Av N = —0.0015 Hz. 
There is some uncertainty in this value due to the estimate 
of the accuracy capability of NBS-III (l<r = 5 X 10~ 13 
[19]) and due to the frequency dispersion of AT(NBS) 
since it was calibrated (May 1969). The actual uncertainty 
(one sigma) during the time of the measurement was 
±9 X 1(T 13 . 6 

We have already discussed, the fact that the wall-shift 
correction was obtained by a comparison between the 
portable SAO maser and the HU reference maser. The 
wall shift of the HU maser was measured in early 1969; 
a one-sigma confidence interval of ±0.5 mHz was assigned 
to this measurement [3]. However, we have to consider 
the possibility of secular changes of the physical pro- 
perties of the wall coating since that time. Preliminary 
measurements at HU indicate that the associated un- 
certainty could be as large as 2 mHz. A correction of 
—2 mHz seems most likely [28] and would put the value 
of experiment 1 very close to that of experiment 2. We 
are not applying this correction, however, because we 
do not yet consider this result conclusive. The combination 
of the original measurement uncertainty and this bulb- 
aging effect thus leads to a total uncertainty of ±2 mHz. 


6 This value has been confirmed by comparisons between the 

* This technique was also used in the comparison between the rates of AT(NBS) and AT(U.S. Naval Observatory) as well as 

SAO maser and the HU reference maser. AT(NRC) over the past year. 
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TABLE IV 

Summary of Results From Experiment 1 


■ 

vs 

VB 

Afutc 

Ax// 

Apw 

Ap d 

Ap z 

Observed 

quantity 

synthesizer 

setting 

electronic 

counter 

measurement 

-299.995 X 10- 10 

Fig. 2 

beat frequency 
at HU 

Fhulb = 

47 “C 

(320.15°K) 

fz = 1020 Hz 

Bias 

corrections 

(Hz) 

+405 794.4200 

-0.1193 

-42.6115 

-0.0015 

+0.0202 

+0.0627 

-0.0015 

Bias 

uncertainties 

(Hz) 


±0.0004 

— 

±0.0013 

±0.0020 

±0.0001 

±0.00002 


v„ = 1420 405 751.7691 ± 0.0024 Hz 


The second-order Doppler correction for the SAO maser 
is calculated from the bulb temperature using (4). The 
temperature was (47 ± 0.5) °C maintained by the auto- 
matic temperature control of the SAO maser. This leads 
to the values for Av D and its uncertainty in Table IV. 
The last correction in Table IV is the frequency change 
due to the Zeeman effect. A Zeeman resonance frequency 
of 1020 Hz was measured, and Av z can be calculated 
from (5) and (6). The Zeeman frequency was typically 
stable to about ±5 Hz. This leads to a negligible un- 
certainty. 

From (8) and Table IV we can now calculate the 
result for v H , the unperturbed hydrogen hyperfine fre- 
quency. Its uncertainty is calculated as the square root 
of the sum of the squared individual contributions since 
they are independent of each other. From 

Experiment 1 : 

v„ = 1420 405 751.7691 ± 0-0024 Hz. 

We now discuss the results of experiment 2. A total 
of 15 measurements of different coatings and bulb sizes 
were made. As coating material we used Du Pont’s TFE 
Teflon, Blend 42, Lot 8842, 34.3 percent solids, purchased 
in 1969. 

The actual coating solution consisted of 50 percent 
Teflon, 40 percent distilled water, and 10 percent Triton 
X-100 (Rohm and Haas). The bulbs were thoroughly 
cleaned with diluted sulfuric acid, and the coating solution 
was applied and then dried by circulating dry nitrogen 
gas through the bulb. The baking was done at 380°C 
for one-half hour while oxygen was circulated through 
the bulb. The bulbs were heated up rapidly in about 
15 minutes; the cooling after baking took about 1 hour 
although the initial drop in temperature to below 300°C 
occurred rather rapidly within a few minutes. The bulbs 
were heated up to above 300°C immediately before 
setting them into the maser. They were put hot into 
the maser, which was then pumped down immediately. 
Some bulbs (indicated by a dagger in Table V) received 
a second Teflon coating on top of their first one. The 
procedure was identical to the one described previously, 
including the cleaning of the first coat with acid. No 
problems were encountered in any of the many coating 


procedures. All coatings, nrst and second ones, passed 
the “water drop test” (tested after their use in the maser). 
The coatings were transparent and appeared only slightly 
milky. All bulbs, except bulb 1, were new and never coated 
before. 

The first column in Table V gives the identification 
numbers of the bulbs. The second column lists the dates 
of their measurements in the hydrogen maser. In the 
third column we find several remarks; in the fourth, the 
nominal diameter is given; and in the fifth column we 
find the accurate inverse diameter that was determined 
from a volume measurement. The sixth column gives 
the temperature of the bulb (actually the temperature 
of the cavity was measured) at the time of the measure- 
ment. The NBS-H2 maser used in this experiment had 
no temperature control of its own. Therefore the tem- 
peratures for the various measurements differ slightly. 
The seventh column lists the Zeeman resonance frequency. 
The eighth column gives the beat frequency v B of the 
tuned maser. In the ninth column the beat frequency 
v B (corrected) is corrected for the offsets Av N of the 
NASA hydrogen standard according to Fig. 2. The last 
column lists the standard deviation of several (more 
than ten, typically) independent tunings of the maser. 
A plot of the beat frequency v B (without the Av N cor- 
rection) as a function of the inverse bulb diameter is 
shown in Fig. 3. No oscillations were possible with bulb 
no. 4 since this bulb caused a significant reduction of 
the cavity Q. Bulbs 1*, 5*, 6*, and 10* were excluded 
from the final analysis (indicated by the asterisk). 

They all showed obvious differences from all other 
bulbs, as we now describe. Bulb 1* was a previously 
used one, stripped of its old coating with hydrofluoric 
acid, and recoated. It showed up in the graph with a 
far too high v B , i.e., a larger wall shift than the other 
bulbs. This may be explained by its surface roughness 
leading to an effectively larger surface area. A second 
coating (If) apparently smoothed out this roughness 
and brought this bulb in line with the others. 

Bulbs 5* and 6* received extremely thin, almost 
invisible coatings, but they passed the water drop test. 
They responded with somewhat lower values of v B , i.e., 
less wall shift than the rest of the bulbs. A second coat 
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TABLE V 

Determination of the Unperturbed Hydrogen Frequency 
Experiment 2 





Nominal 

Inverse 




VB 


Bulb 

Date of 


Diameter 

Diameter 

T 

fz 

VB 

(Hz) 

1 o 

No. 

Measurement 

Remarks 

(inch) 

(cm -1 ) 

(°C) 

(Hz) 

(Hz) 

(corrected) 

(mHz) 

1* 

Nov. 25, ’69 

old bulb, 
glass etched 

6 

0.0608 

24.0 

252 

-0.7191 

-0.7206 

1.4 

It 

Dec. 16, ’69 

old bulb, 
glass etched 

6 

0.0608 

23.7 

256 

-0.7094 

-0.7107 

0.9 

2 

Dec. 5, ’69 

6 

0.0599 

24.1 

233 

-0.7114 

-0.7129 

0.9 

2 f 

Jan. 16, ’70 


6 

0.0599 

24.1 

215 

-0.7105 

-0.7124 

1.0 

3 

Dec. 10 ,’69 


6 

0.0597 

23.7 

271 

-0.7095 

-0.7109 

0.8 

4 

very heavy bulb 

5 


no oscillations possible 



5* 

Nov. 28, ’69 

very thin coating 

5 

0.0786 

23.8 

255 

-0.7172 

-0.7186 

0.7 

5* 

Jan. 9, ’70 

5 

0.0786 

24.0 

262 

-0.7210 

-0.7228 

1.5 

6* 

Dec. 1, ’69 

very thin coating 

4 

0.0975 

24.0 

265 

-0.7226 

-0.7241 

2.S 

G f 

Jan. 5, ’70 

4 

0.0975 

23.8 

277 

-0.7313 

-0.7328 

1.3 

7 

Dec. 19, ’69 


4 

0.101 

24.0 

255 

-0.7319 

-0.7336 

1.0 

7 1 

Dec. 24, ’69 


4 

0.101 

23.8 

250 

-0. 7299 

-0.7316 

0.9 

8 

Jan. 8, ’70 


4 

0.103 

24.1 

263 

-0.7330 

-0.7347 

0.6 

9 

Jan. 13, ’70 


6 

0.0655 

23.8 

252 

-0.7106 

-0.7125 

1.1 

10* 

Dec. 29, ’69 

no spin-exchange 
tuning possible 

3 

0.134 

23.9 

248 

-0.741 

-0.743 

7.0 

11 

Jan. 28, ’70 


8 

0.0504 

23.7 

244 

-0.7052 

-0.7069 

1.0 


* Indicates zero weight given in final analysis, 
t Indicates a second Tenon coating. 



Fig. 3. Wall-shift measurement and hydrogen-cesium comparison 
from experiment 2 using the NBS-H2 hydrogen maser. As plotted, 
the beat frequency v B does not contain the Apy correction. Bee 
Table V for both v B and v B (corrected). See Fig. 2 for the At.v 
correction. 


(5f and 6f) brought them to values compatible with the 
majority of the bulbs. We have no good explanation 
for this effect. 

Bulb 10* was so small that oscillations were barely 
possible. For this reason the spin-exchange tuning pro- 
cedure, which required a reduction in beam intensity, 
could not be used. The uncertainty in its v B value was 
therefore prohibitively large for its inclusion in the final 
analysis. 

The “good” bulbs were used in a linear, least-squares 
fit whereby each measurement point was weighted ac- 
cording to its variance ( a 2 from Table V). The intercept 
at 1/D = 0, which corresponds to Ai/, r = 0, was obtained 
as v B (corrected) = —0.6798 Hz. The standard deviation 


of the n = 11 points from the straight line was calculated 
to be a = 0.00108 Hz. We can now estimate the confidence 
interval w of the intercept at 1/D = 0 from [29], 

' w 

where a is the confidence level, t is the student t, and 
x = 1/D. For one-sigma confidence we calculate w = 
0.0012 Hz. 

Table VI summarizes this result and all other pertinent 
data analogously to Table IV. The synthesizer was set 
at v s = 405 795.0000 Hz, and the beat frequency v B 
(corrected) in Table VI is the frequency at the intercept 
and is thus the frequency for Av w = 0. It contains the 
correction Av s that was already applied in Table V. The 
uncertainty associated with Av N is the same as in Table IV 
and was discussed there. For the Doppler correction Av D 
we use the mean value for the temperatures listed in 
Table' V. The uncertainty of the temperature measure- 
ment is again estimated at ±0.5°C, which well includes 
the systematic differences of the temperatures of the 
individual measurements. The mean Zeeman frequency 
of all measurements (Table V) is f z SB 245 Hz (175 nG), 
which leads to a correction of only 0.1 mHz. The indi- 
vidual measurements of / z spread around this mean by at 
most ±30 Hz, which corresponds to a negligible variation 
of the bias of ±0.02 mHz. This value is listed as bias un- 
certainty of Av z in Table VI. 

From (8) and Table VI we can now calculate v H , the 
unperturbed hydrogen hyperfine frequency. The un- 
certainty is the square root of the sum of the squared 
individual contributions, analogously to experiment 1. 
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TABLE VI 

Summary of Results From Experiment 2 



VS 

VB 

(corrected) 

AfuTC 

Aw 

Ar^ 

Ar c 

A p 2 

Observed 

quantity 

synthesizer 

setting 

Fig. 2 
and 

Table V 

-299.995 X 10“ 10 

Fig. 2 

Fig. 3 

Tbulb = 
23.9 °C 
(297.05 °K) 

fz = 245 Hz 

Bias 

corrections 

(Hz) 

+405 795.0000 

-0.6798 

-42.6115 

included 
in v B 

(corrected) 

included 
in v B 

(corrected) 

+0.0581 

-0.0001 

Bias 

uncertainties 

(Hz) 

— 

±0.0012 


±0.0013 


±0.0001 

j ±0.00002 


v „ = 1420 405 751.7667 ± 0.0018 Hz 


Experiment 2: 

v„ = 1420 405 751.7667 ± 0.0018 Hz. 

Fig. 3 and Table V allow a calculation of the wall- 
shift coefficient A'(at 24°C), where 

K' = K[ 1 + a (24 - 40)]. (10) 

For experiment 2 we find (for TFE-Teflon) 

K' = +0.528 Hz -cm. 

From Table I (experiment 1) we calculate for the value 
of Zitzewitz et al. (FEP-Teflon) at the same temperature 
using K' = +1.192K 

K' = +0.460 Hz • cm. 

The difference between these values is only about 13 
percent despite the fact that two different kinds of Teflon 
were measured: the homopoly mer TFE (tetrafluoroe- 
thylene polymer) and the copolymer FEP (TFE plus 
hexafluoropropylene) . 

IV. Conclusions 

The values for v„ from experiments 1 and 2 are well 
within each other’s error limits. The uncertainties assigned 
to the two values are not much different. We therefore 
take the mean of both as the most probable value of the 
unperturbed hydrogen transition frequency 

v„ = 1420 405 751.768 ± 0.002 Hz. 

This result is based on the two most extensive wall- 
shift evaluations that have been carried out so far, as 
evidenced in Table I. Furthermore, the hydrogen-cesium 
comparison was done while a direct link existed to the 
wall-shift correction data. In experiment 1 this link was 
provided by the preservation of the results of the Harvard 
wall-shift experiment [3] in the HU reference maser, 
and the calibration of the SAO maser against this ref- 
erence. In experiment 2 the link was established by the 
simultaneous, integral wall-shift/hydrogen-cesium com- 
parison experiment. In addition we used one of the best 
evaluated cesium standards currently in existence as 


our reference (NBS-III). We feel that our quoted error 
limits, as documented in this paper, are realistic and 
meaningful. 

The only question not yet completely settled is the 
problem of thermalization in the bulb. Although there 
is evidence [22], [23] that a somewhat incomplete ther- 
malization will not cause an observable bias correction 
today, we encourage experimentation towards clarifying 
this point in view of the possibility of future, more accurate 
measurements. 

A comparison of our results with those obtained pre- 
viously is made in Table II. The range of values for v B 
totaling about 6 X 10" 11 (90 mHz) far exceeds any quoted 
uncertainties and therefore requires some critical dis- 
cussion. We surely can not explain deviations of this 
magnitude by differences in the cesium reference standards 
alone. Cesium beam tubes are extensively studied and 
intercompared devices. The uncertainty to be attributed 
to the cesium standards that have been used could be as. 
high as 2 X 10~ u (28 mHz) before 1964, 1 X 10" 11 in 
1964, 3 to 5 X 10~ 12 (4-7 mHz) between 1965 and 1966, 
and 1 to 3 X 10~ 12 (1.4-4 mHz) since then. A larger 
uncertainty in the wall-shift correction as compared to 
those quoted in the past is the most likely explanation 
for most of the discrepancies. It is evident that the wall 
shift is not a constant of nature, but largely depends on 
the kind of Teflon (including differences from year to 
year due to the manufacturing process) and on the many 
complex parameters involved in the actual coating pro- 
cedure. This is evidenced in Table I (compare also [3]), 
in the exclusion of bulbs from the wall-shift evaluation 
as indicated by the asterisks in Table V, and in [3] where 
we observed wall shifts that were up to 50 percent different 
from the wall shifts of “good” bulbs. Two statements 
can be made in consequence: 1) using only a few bulbs 
(e.g., two) for a wall-shift measurement can lead to 
highly erroneous results, and 2) the use of a wall-shift 
correction not directly related to the hydrogen maser 
that is used in a v„ determination may also lead to errors 
far exceeding those quoted in Table II. If we check the 
values ir-tnd in Table II against the above two state- 
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ments, we find that the majority of the values do not 
have the direct experimental link required by statement 
2) between the bulb used in the hydrogen maser and 
the applied wall-shift correction. Furthermore, they are 
correlated through the almost exclusive use of one par- 
ticular wall-shift evaluation that we believe to be a 
valid one at the time of its measurement but that should 
not be used in the sense of a universal correction (state- 
ment 2)). Only the values of Crampton et al. and Menoud 
and Racine are based on a direct relationship between 
the hydrogen-cesium measurement and a wall-shift 
determination; unfortunately, however, only two bulbs 
were used in each case and statement 1) applies. The 
only values not affected by either statement are those 
of this paper. 

The discussion of discrepancies in past values of v H 
that is given above should not, however, be interpreted 
to imply that a frequency standard based on hydrogen 
storage is necessarily inferior to a cesium standard. 
First, the uncertainties due to the cesium comparison 
and those due to the wall-shift correction contribute 
about equally to the uncertainty that we assign to our 
measurements of v„. Furthermore, a more recent meas- 
urement of v jf has been reported [30], which was based 
on an independent wall-shift evaluation and which was 
also referenced with transfer standards to the cesium 
standard NBS-III. This new value of v H is almost iden- 
tical with the results of this paper. 

Second, new experimental developments that promise 
a reduction in the wall-shift limitation have evolved 
recently. They include the big storage vessel [31], the 
deformable storage bulb [28], and the operation at elevated 
bulb temperatures [30], [32], Thus, hydrogen storage 
devices should be considered serious competitors for 
the primary standard of the future, a place now held 
by the cesium beam standard. 
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System Automatically Tunes Hydrogen Masers 
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An automatic tuning system has been developed 
that permits frequency synchronization between two 
hydrogen masers. Although the system was originally 
designed to match a space-borne clock performance 
with that of a ground-based clock to test the so-called 
“red shift” (or effect of gravity on time predicted in 
Einstein’s general theory of relativity), it should be 
of interest to organizations concerned with industrial 
and educational research programs. Additionally, 
this system, used in conjunction with radio astronomy 
for long-baseline interferometer experiments, should 
interest the astrophysical community as a new tool 
for investigation of distant phenomena in the universe. 

A typical tuning run is illustrated in the figure. The 
tuning varactor voltage for maser 1 is recorded on the 
upper track of a strip chart recorder, and that for 
maser 2 on the lower track. Full scale for both tracks 
is 2.0 V and each major time increment in the time 
axis is 10 minutes. The center track displays an analog 


record of a digital measurement of the beat frequency 
between the two masers with a full-scale range of 
2.5 x 10 ‘ 1 * . Synthesizers in phase-lock loops in sys- 
tem electronics are offset by 0.6 Hz and 100-period 
averages are taken to provide an observation time for 
each measurement of 166 seconds. Maser 1 is delib- 
erately detuned by arbitrarily setting the tuning diode 
voltage to 5.0 V. The corresponding initial frequency 
error was about 3.8 x 10' 12 . Similarly, maser 2 was 
offset about 5.5 x 10' 12 by setting its tuning diode 
voltage to an arbitrary 4.0 V. 

The loop gain on each tuner is deliberately reduced 
to effectively demonstrate the tuning operation. The 
correction process operates smoothly, as shown by the 
recording, alternating from maser 1 to maser 2. After 
150 minutes, the initial total error of 8 x 10 12 is 
reduced to less than 5 x 10' 13 . Of particular interest 
are the fluctuations in the beat period, due to pres- 
sure changes, which can be seen clearly for the first 
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S umma rv 

The use of high molecular weight Poly- 
t e t ra f luor oe t hy lene (PTFE) as a wall coat- 
ing for the atomic hydrogen maser appears 
to offer some real advantages. By operat- 
ing the maser bulb at a temperature where 
the wall shift is zero the effects of bulb 
size, shape and the surface texture ofthe 
coating can be eliminated. Coatings of 
PTFE can be applied in two quite distin- 
guishable conditions; . quenched and 
annealed. It is believed that the quenched 
coatings offer good reproductivity of the 
zero wall shift temperature. The annealed 
coating has a higher zero wall shift 
temperature and appears to offer somewhat 
lower relaxation rate than the quenched 
coa ting . 

The results from some initial experi- 
ments are presented to support these 
contentions. Further work on wall coatings 
is in progress . 

Introduction 

Recent studies- 1 of the temperature 
dependence of the wall shift in the atomic, 
hydrogen maser have shown that the average 
phase shift per collision can be made 
positive or negative with respect to the 
phase of the free hydrogen oscillating 
dipole moment. We expect that the temper- 
ature at which a zero average phase shift 
per collision occurs is a property of the 
chemical nature of the walls and not of the 
mechanical surface conditions. By operating 
at this temperature, masers with bulbs of 
any size and surface condition should 
oscillate at the same frequency, provided 
that their surfaces are of the same type of 
material and free of contamination. The 
surface conditions will affect the relaxa- 
tion rate and for this reason the wall 
collision rate should be kept as small as 
possible. 

High molecular weight Po ly t e tr af luor o- 
ethylene (PTFE) films cast from water 
dispersed resin are of particular interest 
since these materials have, on the average, 
about ten times the molecular weight of the 
t e t r af luo ro e t hy lene (FEP) copolymer most 
often used as a wall coating material. The 
PTFE material should have fewer end groups 
exposed to the atomic hydrogen in the bulb. 


The end groups, chiefly CF 2 H can contrib- 
ute a significant area of sites contain- 
ing bound hydrogen that result in larger 
negative values of phase shift per col- 
lision than the fluorine sites and can 
also contribute substantially to recom- 
bination processes (and consequently to 
relaxation) since the end-group hydrogen 
atom is much less energetically bound 
than the fluorine atoms. The tetra- 
f luoroethylene homopolymer can be almost 
completely crystallized by the very care- 
ful and slow annealing through the range 
340°C - 295° C. This should result in a 
30-fold reduction in porosity to gases 
and the possibility of much reduced 
hydrogen collision rates. 

Studies of PTFE^ have shown that the 
change in the ordering of the amorphous 
PTFE teflon known as the alpha transition 
is at a higher temperature than the 
temperature for zero wall shift making 
this temperature a viable operating point. 

Experimental Techniques 
The ' Reference Oscillator 

A reference standard was used to 
obtain the wall shift data and relate 
the data taken over a time interval of 
several months. The reference oscillator 
used to provide the frequency datum was 
compared with the Harvard reference maser 
and with the N.B.S. Maser using a travel- 
ing maser frequency standard. This 
measurement is described by Hellwig et al 
(1970) 3 and is part of an experiment to 
determine the unperturbed Hydrogen Hyper- 
fine Transition Frequency. The wall shift 
of the reference oscillator was determined 
to be -.022 + .002 Hz at 319°K. This 
value has been used in all the measurements. 

Working Maser 

The working maser used for measuring 
the effect of temperature on wall shift 
is a modified Varian H-10 series maser. 

The r.f. power supply to the dissociator 
discharge was made to be controllable 
so as to permit variation of the 

* This work is supported by NASA under 
Contract NSR-09-015-098 . 



dissociator efficiency and consequently the 
flux rate of atoms into the bulb. By using 
automatic cavity tuning techniques 4 much of 
the effort required to tune the maser has 
been avoided. 

Serious problems with magnetic quench- 
ing were encountered at high cavity and 
bell jar temperature due to thermal gradi- 
ents across junctions where different 
metals came into contact. All the oven 
cylinders, the solenoid and the bell jar 
components were electrically insulated 
from each other. The problem has been 
almost completely eliminated and it is 
possible to operate over a range of mag- 
netic fields from 0.7 to 100 mill ioers ted s 
with very little variation of signal ampli- 
tude when the maser is oscillating near 
threshold. For wall relaxation measure- 
ments the residual magnetic quenching is 
still not entirely negligible. 


The maser bulb and cavity were both 
enclosed in the maser vacuum system. The 
bulb temperature was measured by measuring 
the cavity temperature by means of a ther- 
mocouple after allowing a 4-day period for 
thermal equilibrium of the bulb-cavity- 
bell jar assembly. 

Maser Bulbs 


bulb during the curing cycle. The oxygen 
was of high purity and was passed through 
a cold trap at -78.5°C (dry-ice) prior to 
entering the bulb. At the end of the 20 
minute curing cycle the bulb was quickly 
removed from the oven and cooled by a 
draft of room temperature air. In this 
condition, the bulb is considered to be 
quenched . 

The subsequent annealing cycle that is 
required for crystallizing the PTFE coat- 
ing consists of placing the bulb in a 
specially controlled oven with the same 
oxygen flow and cooling the oven, in a 
steady rate, from 340°C to 290°C. The 
bulb was removed from the oven at 2 9 0 ° C 
and cooled in air to room temperature. 

Experimental Data 

Relation of Wall Shift to Phase Shift per 
Collision 

The quantity of interest is the 
average phase shift per collision. 

This quantity is related to the wall 
shift by the following expression: 



Bulbs of 7" and 4" diameter were used 
for taking wall shift data as a function of 
temperature . 

These bulbs were equipped with integral 
collimators that were coated with the bulbs 
in order to be assured that hydrogen atoms 
would encounter only one kind of teflon 
surface. Figure 1 shows a bulb of 4" 
diameter with integral collimator. 

Coating Procedures 


is the wall shift 

and yft, is the wall collision rate. The 
collision rate, is related to A? , 

the average velocity and the mean free 
distance by 








For a sphere and v ' 

The resulting expression relating average 
phase shift with bulb radius, wall shift 
and temperature is 


The bulbs were chemically cleaned by 
using a water solution of nitric acid and 
rinsed with several changes of hot distilled 
water. No attempt was made to roughen the 
interior surface of the bulbs. 

The coating material consisting of a 
mixture of 10% Triton X-100, 50% Teflon 42 
(duPont), and 40% distilled water was care- 
fully poured into the bulb to be coated, 
swirled around to cover all the interior 
surface and drained out. The wet bulb was 
then placed on a rack with a small diameter 
tube extending upwards into the collimator. 

A flow of approximately 1 cm3/minute of dry 
nitrogen, at room temperature was allowed to 
flow into the bulb to dry the coating. 

The bulb was then cured in a horizontal 
air flow oven for 20 minutes at 360°C . A 
very small flow of oxygen was piped into the 
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The data are plotted in figure 2 in 
terms of average phase shift per wall 
collision. B: is seen that the 4" and 7" 
quenched bulb data lie very nearly on the 
same straight line. The least squares fit 
on each set of points is as follows: 

7" bulb. 

= -53.36 + 1.37 + (0.1503+.004 
= 0 at 35 5 ° K 

4" bulb 

= -56.53 + 1.5 + (0. 159+. 005)T* 
■ 0 at 355. 1°K. 





The same 4" bulb was then annealed and 
data on it are shown in the same figure. 

Since we wanted to get data in time for the 
conference we did not have time to allow the 
bulb and cavity to come to thermal equilib- 
rium and consequently the scatter is larger 
than for the previous data. We have for the 
4" annealed bulb 

&<*> = -39.09 + 3.67 + (.0903+. Oil) 

= 0 at 432. 9°K. 

In both cases the data excluded the 5 
circled points in figure 2. 

Summary of Data 

Table 1 is a summary of these prelimin- 
ary data and a comparison with other data. 

It should be noted here that the Teflon 
42 used in all the data shown in this table 
is from the same sample of Teflon 42 ob- 
tained by the Smithsonian As trophy s ica 1 
Observatory from E.I. duPont de Nemours and 
Company on July 29, 1969. This sample, which 
was donated to SAO by duPont, is identified 
as Blend #8892, % solids 34.3. It should 
further be noted that the coating procedure 
used by Hellwig"* was similar to the quenched 
procedure described above except for the 
very rapid cooling employed in the present 
experiment . 

It is possible to extrapolate to zero 
wall shift from data shown in the "quenched" 
plot of Figure 2. By using the least 
square fit to a straight line for the 4" and 
7" bulbs quenched data, recalculating the 
wall shift at 4 temperatures and plotting the 
results vs. the reciprocal of the bulb diame- 
ter one obtains the plot shown in Figure 3. 

Figure 3 illustrates an interesting 
point and leads us to a new method of de- 
termining the wall shift. This figure is, 
in fact, a plot of wall shift versus col- 
lision rate. Collision rate depends on the 
bulb geometry and on the collimator geometry; 
the effect of the latter is still not well 
unders tood . 

In making wall shift vs collision rate 
plots using several bulbs we make the follow- 
ing two assumptions: 

1. That the walls of the two bulbs are 
exactly the same; 

2. That we know the collision rate of 
the hydrogen atoms in each bulb. 


it appears that we are making a consis- 
tent error in the prediction of the 
collision rate. However, if we use our 
knowledge of the teimerature at which 
the observed wall shifts are equal we 
have no longer to depend on our know- 
ledge of the collision rate. We depend 
only on the similarity of the teflon 
used in the two bulbs. The data is in 
the form 


- Atk + 


(®‘ # 


where AVe is an unknown offset error 
and A>>u) is the wall shift. X is mean 
free path of atoms in the bulb and col- 
limator. By changing X, either by 
changing bulbs (or by usinga deformable 
bulb as suggested by Brenner’) we see that 
the only condition making '')('T) the same 
for the two values of X is to make 


= O. 


<j>(v) = 

To illustrate this technique we have 
reduced the data discussed earlier in a 
different way. 

We find that »/ 7- \ 

if? . 

rpA fpf*- (/rW A- 

is quite linear with temperature when Ah)e 
is small. For the two quenched PTFE 
bulbs, we have that — 

- (?, -toil T 

T* W'j Xj (1) 

>) f CP) * 4\k -f/^4, r Ci.+ c< *' r 

T**- 7^ (2) 

Under the assumption that A^)e is the 
same foi^bothj[l) and (2) and that 

djQr) =fycr)- 

crossover occurs where (pQPj^O at a 
temperature y _ Cx'Ct/ur^x. - 

Recalculating the calibration error, we 
have 




From the same data used to plot Figure 2, 
we obtain a temperature of intercept of 
355. 3°K and a 4'^(error)t0 . 2 x 10"3 Hz. 

Wall Relaxation 


The plot of wall shift vs 1/0 at 
different temperatures shows us that the 
convergence point is not at 1/0 = 0 but 

somewhere near I/O = .008. Here the 
convergence is good and lies within 2 mHz 
of our calibration. Judging from the data 


So far, we have made wall relaxation 
measurements only using 4" bulbs. Our 
results confirm the data obtained by 
Zitzewitz for PTFE . At temperatures up 
to 100°C we see no increase in wall re- 
laxation rate. Relaxation rate measure- 



merits are being continued at the Smithson- 
ian As trophy s ical Observatory using larger 
'> lbs with very long bulb storage time. 

Cone 1 us ions 


(1) From the data taken so far we con- 
clude that it is possible to operate a 
maser with its bulb at 355°K (82°C) and 
substantially reduce the wall shift and 
its contribution to the uncertainty of 

the maser output frequency. 

(2) It is possible to eliminate the 
uncertainty due to collision rate in the 
determination of wall shift by measuring 
wall shift as a function of bulb temperature 
using several bulbs of different sizes. 
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Table I 


Parameter 


Present Data 


Previous Data 


Average phase shift 
Temperature 


Zitzewitz 1964 

O. 15 + .004 Rad / °K 0.1294 

P. T.F.E. quenched 0.1248 


O. 09 + .011 Rad/°K 

P. T.F.E. annealed 


Zero phase shift 
Temperature 


355°K quenched 375°K 

4 33°K annealed 369°K 


Wall Shift. K 522 Hz-Cm 

mHz Cm at 297°K 


630.6 Hz-Cm 

640 . 7 
570.9 

544.8 


Hellwig 1970 
528 Hz-Cm 



0.040 
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AVERAGE PHASE SHIFT PER COLLISION 
4) IN UNITS OF 10 ' 6 RADIAN 
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Abstract 

The wall shift of the hydrogen maser depends on the 
collision rate and the average phase shift per collision of the 
atoms in the storage bulb. By operating at the temperature 
where the average phase shift per collision is zero it is possible 
to remove the wall shift. The determination of this temperature 
is made by measuring the wall shift as a function of temperature 
for two similarly coated bulbs offering different collision rates 
to the stored atoms. The intercept of the wall shift vs. tem- 
perature curve for the two bulbs occui'3 at the zero phase shift 
per collision temperature. To eliminate the possible differences 
in wall coatings the use of a flexible storage bulb is recom- 
mended. 

An experimental example of the technique is given using 
two separate bulbs. The data is in good agreement with a 
previous measurement of the unperturbed hydrogen hyper- 
fine frequency by H. Hellwig et at. 

One of the most vexing problems affecting the 
accuracy of the output frequency of the atomic 
hydrogen maser is the shift in frequency due to 
collisions of the hydrogen atoms with the storage bulb 
walls (wall shift). This shift results from the defor- 
mation of the electron cloud of the hydrogen atom 
during the collision process and can be represented as 
an average advance or retard per collision of the phase 
of the oscillating magnetic dipole moment of the atom 
with respect to the phase of the r.f. magnetic field in 
the maser cavity. 

The usual method of measuring wall shift [1] has 
been to coat a number of spherical storage bulbs of 
different diameters with the wall coating material and 
to measure the maser output frequency associated 
with each bulb. 

Plotting the output frequency as a function of 
collision rate (inversely proportional to bulb diameter 
for spherical bulbs) and extrapolating to zero collision 
rate gives the zero wall shift output frequency of the 
maser. 

The chief difficulties and uncertainties in making 
these measurements lie (1) in estimating the collision 
rate of the atoms in the bulb and (2) in the assumption 
that all the bulbs are similarly coated. Since the micro- 
scopic texture of the wall coating can affect the 
collision rate, these two uncertainties are not always 
independent. Furthermore, the effect of the bulb 
Collimator on the collision rate is still not well under- 
stood. 

Recent studies of wall shift versus temperature by 
Zitzewitz [2] and Vessot and Levine [11] have shown 
that there is a temperature dependence in the average 
phase shift per collision that goes from negative, to 
positive values. In particular, there is a temperature 

* This work is supported by the National Aeronaut ics and 
Space Administration under contract NSU-0!>-uir>-0!>S. 


at which the effect is zero. For polytetrafluoroethylene 
(PTFE, Du Pont Teflon 42) applied by rapid sintering 
and quenching, we have found that this temperature 
is about 82 °C. 

We wish to report the following procedure for 
determining the zero wall shift frequency which 
eliminates the error due to the uncertainty in the 
collision rate. 

The output frequency of the maser is measured 
over a range of temperatures usually from 23 °C to 
150 °CJ using two bulbs providing substantially 
different collision rates that are coated in the same 
manner with the same material. (A better method 
might be to use a single bulb having a flexible wall 
allowing a change in volume without a change in 
surface, as suggested by Brenner [4], This has the 
advantage of removing the uncertainty associated 
with the chemistry of the wall surfaces of two separate 
bulbs.) 

The wall shift, A v„ (T) is given by A r w (T) = 
{T) w r here v c is the collision rate and (T) is 
the average phase shift per collision. In the case of 
spherical bulb, one can replace v c by where v = 
8 fcT VI ' 

1 ‘is the average velocity of the atom, R is the 

7i m j 

bulb radius, k is the Boltzmann constant, m is the 
atomic mass, and T is the absolute temperature. 

The only way that the output frequency of the 
maser can be made independent of collision rate is to 
make <f> (T) zero. If the output frequency of the maser 
is plotted against temperature under two (or more) 
different bulb diameters, or two different configura- 
tions of a flexible bulb, the intercept of the tw r o plots 
will be at the zero wall shift maser output frequency. 

We offer the following example of this technique 
applying it to two bulbs coated with PTFE, and 
quenched rapidly after sintering the coating. The 
wall shift data are taken in the form: v (T) = At’ e + 
At’ w ( T) where Ar c is the frequency difference 
between the two masers after application of the 
magnetic field, and second order Doppler corrections 
to both masers and a wall shift correction to the 
reference maser. A»v is the residual offset of the 
reference oscillator, and Ai\ v (!’) is the temperature 
dependent wall shift. Since the effect of the second 
order Doppler effect in the temperature varying 
maser lias been accounted lor, A>v is a constant. 
When A iv is small we find that 

■•('A A iv . i U T \% 1 

'[■'!> ' " T' r -~ \ j 7"' ; ~ 'A 
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where X is the mean free path between collisions. The 
x’tT) 

expression varies approximately linearly with T, 

and thus calculation of the intercept is readily 
accomplished. 

We have that for two different bulbs 1 and 2, 


v x (T) 

- + 1 
J’»/2 r \ 

f J*J\ 

'1, fail') 

= c. 

+ oq 

T 

(1) 

ynt% 

{ nb/i j 

" fa 


and 








V 2 (T) 

-- + 1 
yi/3 1 \ 

f 2k ) 

‘h fa(T) 

— n 

+ *2 

T . 

(2) 

/pi/a 

{ nhn J 

1 fa 

— L/ 2 


If we assume that Av e is the same for (1) and (2) 
and that fa ( T ) = fa (T), cross-over occurs where 
4> (T) ~ 0, at a temperature T g — — — — . 

fti OC g 



We can thus determine the offset Av e , 

Av e = TJU (C\ + *T 0 ). 

Fig. 1 shows the values of plotted as a 

function of T for two bulbs of IS and 10 cm diameter. 
The lines drawn are the least squares fit to the data. 
The temperature at which the lines cross is 355.28 ± 

0.2 K and the value of i >(T) at the crossing point is 
(0.2 ± 1.1) x HP 3 Hz. 

The quantity r(T) had been corrected for the 
following oll'sets in frequency: 

1. Magnetic field effect for each maser. 

2. The second 1 order Doppler elicit for each maser. 

3. The wall shift of the reference maser. 

The first two corrections are well known functions 
of temperature and magnetic field |(i|. However, 
considerable pains were Liken to calibrate the wall 
shift of the reference maser. This work involved 
making direct frequency comparisons with wall shift 
determined at Harvard by Xit/.ewitz and at NI>S by 


Hell wig by using a traveling maser frequency standard. 
From this determination the wall shift of the SAO 
reference was determined to bo —.0222 ± .002 Jlz. 
(This maser is equipped with a 16.30 cm spherical 
bulb coated with DuPont FEP-120 teflon and 
operates at 45.5 °C.) 

The value of v c derived from the crossing point of 
the lines in Fig. 1 is less than 0.2 x 10~ 3 IIz and 
provides a further, independent, confirmation for the 
value of the unperturbed hydrogen liyperfine frequency 
given in reference 1 as 1420405751.768 ± .003 Hz, 
with respect to the NFS III cesium standard. 

Using our knowledge of the radius of the 18 and 
10 cm bulbs and calculating the average phase shift 
per collision we obtain the following equations based 
on least squares fits to a straight line. 

For the 18 cm bulb 

fa (T) = - 53.36 ± 1.37 + (0.1503 ± .004) T . 

For the 10 cm bulb 

fa (T) = - 56.53 ± 1.5 + (0.159 ± .005) T . 

The similarity of these two temperature functions 
gives us a measure of the similarity of the wall 
coatings in the two bulbs and allows us to make a 
measurement of the wall shift parameter K ( T) which 
is defined here as 


At 297 Iv we obtain K (T) = 522 Hz • cm. Hell- 
wig [1] obtains a value K (T) = 52S Hz • cm for the 
same PTFE Teflon and finds it to be consistent for a 
large number of bulbs. Note that a limitation on the 
measurement of K is due to the uncertainty in the 
effective diameter, D. 

From these data we conclude the following: 

1. That we can operate a maser with very small 
wall shift. 

2. The zero wall shift temperature can be deter- 
mined as the temperature at which variations in 
collision rate have no effect on the output frequency. 
A flexible bulb maser is ideal for this determination 
since the question of the similarity of the chemical 
nature of the walls of two or more separate bulbs is 
eliminated. 

3. High molecular weight polytetrafluoroethylene 
homopolymer (PTFE) may be a more reproducible 
wall coating material than the tetrafluoroethylcno- 
hexafluoropropylene co-polymer (FEP) that is most 
often used as a maser wall coating material. 

4. The present data confirm the unperturbed 
hydrogen hyperfinc frequency given in reference [1] 
by Hcllwig et al. 

References 

1. Hcllwig, H., Vessot, R.F. C-, Levine, M. \V„ Zitzcwitz. PAY., 
Allan, 1). \\\, Glaze, D. J.: IEEE Trans, lnstr. Meas.. (in 
press). 

2. Zitzewitz, P. W., Uzgiris, E. E.. Ramsey, X. F. : Rev. Soi. 
Instr. 41 (1970). 

3. Vessot, R. F. 0., Levine. M. W.: Proe. 24th Annual Sym- 
posium on Frequency Control, U.S. Army Elect ronics Com- 
mand, Ft. Monmouth. X..I. 

4. Premier, D.: Hull. Am. Rhys. Soe. 14. 943 (19(>9). 

5. Kleppner, J).. Rerg, H. C\. Grumpton. S. IE. Ramsey. N. l'\. 
Vessot, IE, Peters, H. E., Vauier, .1 . : Plivs. Rev. I3S, A 972 
( 1 9tio). 

(i. Vessot. IE. Peters. II., Vauier. .1.. Rcchler. IE. Ilallord. !>.. 
Harrarh, IE. Allan, D., Cla/.e. 1).. Snyder. <\. I’-arncs. •!.. 
fuller. I... Ruddy. L.: IEEE Trans. Inslr. .Meas. I *»• • 

(I'Hili). 






INVENTION DISCLOSURE 


To the attention of the NASA patent Department: 

The herein described invention is submitted in pursuance of Contract 
NSR 09-015-098. 

The invention relates to the Atomic Hydrogen Maser Frequency Standard 

The purpose of the invention is to remove the effect of the shift 
in the maser output frequency resulting from the collisions of the 
hydrogen atoms with the storage bulb walls. 


Brief Description 

Recent work on bulb coatings by zitzewitz at Harvard and by 
Vessot & Levine at Smithsonian Astrophysical Observatory have shown 
that there is a temperature at which the wall shift crosses through 
zero from negative values to positive values as temperature is in- 
creased. More recently (see appended preprint of paper to be pub- 
lished in the Proceedings of the 24th Annual Symposium on Frequency 
Control, USAECOM, Fort Monmouth, N. J. ) . Vessot & Levine show that 
there is a lower temperature zero wall shift point than that found 
by Zitzewitz and they prove that this zero wall shift point is in- 
dependent of the collision rate of the atoms in the bulb. 

This forms the basis for the invention when combined with 
D.' Brenner's idea (see reference in appended preprint) of using a 
flexible bulb. A servo system is used to keep the temperature of 
the maser bulb at a value where no change in frequency is observed 
or changing the collision rate by flexing the bulb so as to change 
its surface to volume ratio. 

The diagram shown as fig. (1) describes the servo system. The 
cycling timer moves the bulb's flexible end in and out, say, every 
5 minutes. The frequency of the maser output is measured when the 
bulb is at low volume and compared automatically with the frequency 
obtained when the bulb is at high volume. The frequency difference 
is used to control the temperature of the bulb. If the bulb is too 
cold the frequency difference (high vol.-low vol.) is positive. 
If the bulb is too hot the frequency difference is negative. 

Figure (2) shows the wall shift vs. temperature as a function of 
the bulb in its two positions. 
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A good way to implement this servo technique is by using a 
bidirectional counter to measure the period of a fixed number, 

(say 10) , of cycles of the slow heat between the maser and the 
tuning reference. The beat' signal operates a gate that allows a 
number of steady "clock" pulses to be counted by the bidirectional 
counter. The number of these pulses will depend on the frequency 
of the beat signal. When the bulb is in its low volume configura- 
tion the bidirectional counter will count "u$ when the gate is open 
during 10 cycles of the beat. The bulb is then changed to high 
volume and the pulses are again gated for 10 cycles, this time the 
counter counts 'down. The residue- in the counter is a measure of 
the frequency shift between the high and low bulb volume config- 
uration and can be used to operate a digital-to-analog device con- 
trolling the bulb temperature. 

In this system the motion of the bulb will affect the tuning 
of the maser cavity "pulling" the frequency of the maser. This 
effect will be taken out by the use of the automatic cavity tuning 
servo system invented by vessot & Levine ("Automatic Tuning of 
Hydrogen Masers" - NASA Tech Brief 69-10452) . The cavity mistuning 
due to the bulb motion can be programmed out by using pre-set tun- 
ing points at the two bulb positions. 

By the use of the temperature tuning servo the effect of the 
frequency shift due to wall collisions can be eliminated. 

This invention was conceived on April 20, 1970. First entry 
and sketches were made on May 27, 1970, in notebook of R.F.C. Vessot, 
Smithsonian Astrophysical Observatory - Used from October 20, 1969. 

This disclosure is written June 1, 1970. 



R. Vessot 

334 Ocean Avenue, Marblehead, Mass. 



M. Levipre 
Big Rock Road, Manchester, Mass. 


This invention was explained by the above identified inventors 
to me and understood on June 1, 1970.. 

‘t y. 

"H m • k . _ . * ./ / > 

Signature of Wi^ne^s liiari 


Signature 



5. TASK IV 

5. 1 Pump and Dissociator Assembly 

With parts transferred to Contract NSR 09-015-098 from the original NASA Contract 
NASW 1337 with Varian Associates, the small-maser pump section, including the r.f. 
dissociator hexapole state selector and the gas-handling system, was assembled and 
tested under simulated operating conditions; a copper oxide target at the position of the 
bulb collimator verified the alignment of the atomic hydrogen beam. The lightweight 
pump with its reduced magnetic field operated successfully; the target showed clearly 
that the alignment of the magnet and its stopping disk is satisfactory. 

5. 2 CER-VIT Cavity Design 

A significant innovation in this task was a new design of the CER-VIT cavity 
resonator for the maser. The original resonator was constructed of five pieces of 
CER-VIT clamped together with sufficient force to withstand the "g" loading effected 
by launching the maser into a 24-hr orbit with a Titan 3C booster. Although this 
design survived sine-sweep testing on three axes at the levels prescribed by NASA, 
the cavity was difficult to tune and delicate to assemble. A new design simpler in 
concept was deemed necessary. 

An invention disclosure, which has subsequently been written up as a patent applica- 
tion, has been filed by NASA. It describes in detail the cavity and its mounting 
system and is reproduced at the end of this section'. 

5. 3 Assembly of Maser Oscillator 

The maser was completely assembled, as shown in Figure 1. Before assembly, 
the r.f. cavity -bulb subassembly was cold-tested and checked for tuning, with a 
helium-filled enclosure to simulate the dielectric qualities of vacuum. Tests were 
conducted of the r.f. system by means of the 0.070" coaxial r.f. lines to be used in 
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Figure 1. Hydrogen maser. Item 1, r.f. discharge tube; item 2, ion pump; item 3, 
hexapole magnet; item 4, magnetic shield; item 5, solenoid; item 6, tuning 
loop; item 7, low-frequency transition coil; item 8, output coupling loop; 
item 9, cavity structure; item 10, quartz storage bulb; item 11, oven 
heaters. 
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the maser. The cavity was then placed on the bell-jar base and fastened with 90 lb 
of preload from three springs at the top of a canister that enclosed the cavity and 
connected it in compression to the bell- jar base. 

Three 0. 070" copper coaxial lines, one each leading to the coupling loop, the 
tuning diode, and the transverse r.f. coil (Zeeman coil), were fastened in good 
thermal contact with (but electrically insulated from) the neck heat station. This 
heat station lies midway between the pump and the bell-jar base, which are connected 
with phosphorus bronze bellows. The heat station is equipped with two thermistor 
sensors, one thermocouple, and, of course, a heater winding. A second heat station 
was installed at the end of the pump nearest the neck; it consists of a heater, a 
thermocouple, and a thermistor. We believe that, with the thermal control on the jar 
and with a single oven surrounding the jar, plus the two-stage control on the neck, 
the temperature of the cavity system can be controlled with a stability of ± 0.0 1°C 
over changes in ambient temperature of ± 10° C. With the present fused-silica cavity, 
the frequency change resulting from a temperature variation of 0. 01°C is about 
±5 Hz. The resultant pulling of the output frequency, assuming a ratio of cavity Q 
to line Q of 2 X 10 5 , is about 7 parts in 10 14 . 

The atomic hydrogen source, hexapole magnet, and beam stop were separately 
assembled to test the alignment of the beam. The cavity bell-jar base and neck 
assembly were joined to the bell-jar assembly, which had been previously made into 
a single potted unit including all the oven solenoids and magnetic shields. The oven 
and shield canisters were closed layer by layer over the bell-jar base. Each oven 
and shield canister in the potted assembly was insulated electrically to prevent thermo- 
electric currents that would cause magnetic -field gradients. The final end cap 
enclosing the nested assembly consists of a layer of moly-permalloy to ensure that 
the magnetic fields from the pump and hexapole magnet are well isolated from the 
three-layer magnetic shield enclosing the cavity. The wires leading to the solenoid, 
three ovens, and neck were all brought out together near, but not touching, the neck 
and led to a groove in the main mounting plate to which the cavity-oven-shield assembly 
is bolted. 
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This plate, with the maser neck protruding through, was then mounted. Again, 
all cables were led to the same location on the mounting plate. The neck vacuum joint 
was secured to the pump, and the hexapole atom-source assembly was placed in the 
pump. 

The main vacuum flange serves as an aligning jig for the hexapole magnet. When 
secured, it provides about 100-lb pressure on the stacked array by deforming a large 
Beleville spring. The maser r.f. discharge electronics were then installed, along with 
the newly designed pressure -control system. Finally, the magnetic-shield caps were 
installed, thus completely enclosing the pump. 

The redesigned gas-handling system includes a gas cylinder, a pressure gauge, 
and a valve; it is mounted to the main vacuum flange. Hydrogen at 1000 psi is led to 
the palladium valve. The cylinder and gas-handling system extends 2" from the main 
vacuum flange and remains within the diameter of the outer pump magnetic shield. 

5.4 Electronics 


The thermal controllers differ from previous practice in three respects: 

A. Standardized thermistors are used as oven sensors, eliminating the need for 
adjustment potentiometers for setting the oven temperatures. The thermistor R-T 
curves are predictable to within ± 0. 5°C, so that the bridge resistors can be preselected 
for a given oven temperature. Since no temperature trimming is required, the refer- 
ence resistors can be mounted inside the thermally controlled enclosure. This tech- 
nique eliminates the resistor temperature coefficients and thermal EMF as a source of 
thermal drift. 

B. Very low-drift operational amplifiers are used to sense the resistance changes 
of the oven thermistors. The amplifier circuitry is designed to reduce the tempera- 
ture variations associated with input-voltage offset to insignificant proportions, so 
that temperature control of the electronic environment is not necessary. One oven 
and its associated electronics and power supplies are eliminated, and the mechanical 
layout and the accessibility for maintenance of the electronics are greatly simplified. 
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C. Pulse 7 <iuration controllers are used on the outer ovens to increase overall 
efficiency and reduce power dissipation in the output transistors. The pulse -duration 
controller can supply up to 50 w to the outer oven, with a maximum transistor dissipa- 
tion of 2. 5 w. ' This very low-power dissipation permits the use of ordinary techniques 
for mounting printed-circuit boards, without any necessity for massive heat sinks. 

The pressure controller is a distinct departure from past practices in that tem- 
perature-compensation techniques will be used to eliminate the Pirani block oven. 

Two carefully matched thermistors, one exposed to hydrogen flow and the other in 
vacuum, will be used in a special self-balancing bridge configuration to minimize the 
effect of ambient temperature changes on the pressure sensor. 

The magnetic-field controllers for the masers are factory-assembled digital 
Kelvin-Varley voltage dividers. The overall voltage-temperature coefficient of 
this circuitry is sufficiently low that temperature stabilization is not necessary. 

•..••• i ■ ' 

A high-efficiency switching controller supplies 15 w to the bell -jar heater, with 
a total power dissipation of approximately 2 w. No heat sinks are required, and all 
controller components are located on the same chassis; thus, many of the intercon- 
necting cables formerly required are now eliminated. The thermistor resistance 
changes are sensed by very low-drift operational amplifiers, which require no ambient 
temperature control over the normal variations of room temperature (± 10 °C). 

Figure 2 is a circuit diagram of the switching controller. 

The pressure controller is an extremely stable, temperature-compensated 
Pirani gauge. One thermistor senses the hydrogen gas flow; a second is pumped 
directly by the ion pump but is in thermal equilibrium with the gas sensor. The 
second thermistor controls the "automatic vacuum zero" circuit shown in Figure 3. 

The new pressure control does not require a temperature-controlled environment 
and has been operating with good stability over normal room -temperature fluctuations. 

The maser is connected to all its controllers — thermal, pressure, field, source, 
and varactor - through a single connector. The cabling has been designed for sim- 
plicity and ease of maintenance; a complete wiring diagram for the thermal control 
assembly is shown in Figure 4. 
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Figure 2. Switching oven controller. 
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Figure 3. Pressure controller. 



















18 


Figure 4. Thermal and pressure control assembly. 





5. 5 Tests of Frequency Stability 


Results of the short-term frequency stability between the new maser package 
and the small maser formerly built for NASA under Contract NAS 8-2604 are shown 
in Figure 5. These data were obtained by locking a 5-MHz oscillator to the NAS 8-2604 
maser by means of a conventional receiver system consisting of a multiplier to 
1400 MHz and a double heterodyne first to 20. 405 MHz and then to 0.405 MHz. The 
final phase-detection and phase -correction signal to the crystal oscillator was obtained 
by phase comparison with a digital synthesizer at 405794. 3 Hz. This puts the maser 
time scale very close to universal time (UT). 
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Figure 5. Stability (two-sample variance) versus averaging time t of new NASA 

maser system. The data have been normalized to one system. The data 
points show the beat period, the number of beat periods per reading, and 
the number of readings required to estimate < 7 . 
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The 5-MHz signal is then led to the new receiver, which employs a different 
system to obtain a 1400-MHz local-oscillator signal. 

The primary power for the maser is + 28 v DC to permit operation on standby 
battery supplies. Two DC-DC converters are required for this mode of operation. 

The first, to supply ± 15 v, will be purchased from Arnold Magnetics. The second 
converter is to provide 2500 v for the ion pump. 

A converter to supply 2500 v for the ion pumps has been constructed. At the 
normal maser operating pressure, the converter requires 28-v output at 0. 8 A, so 
that a standby battery can be used to provide pump power in the event of AC power 
failure. The converter is extremely compact; all the circuitry is enclosed in a shielded 
housing 7" X 4" X 2". 

5. 6 Phase -Lock Receiver 

Figure 6 is a block diagram of the receiver, showing all major subassemblies. 

This receiver differs from earlier designs in that the source for the 1400-MHz first 
local oscillator is a transistor oscillator phase-locked to the 5-MHz precision quartz- 
crystal oscillator. The previous receiver designs used frequency multipliers to 
generate 1400 MHz from 5 MHz; these multipliers have been shown to be a source of 
phase noise and phase drift. The subsystems, less the crystal oscillator and syn- 
thesizer, are housed in a shielded aluminum box 12" X 17" X 4", which can be 
installed within the maser cabinet. 

To make the maser comparison and to evaluate the new receiver, the maser 
phase-lock loop was opened (see "time" position in Figure 6). The digital synthesizer 
was offset to produce a beat signal, which was fed to a counter and printer. Several 
combinations of frequency offset and averaging times were used; the data points in 
Figure 5 show the beat period, the number of periods per reading, and the number of 
readings used for the statistics. The data to 300 sec roughly follow a 1 /t behavior, 
indicating, as expected, that the dominating noise process is additive white-phase 
noise. Some indication of flattening out is observed at 1000 sec; this is due to sys- 
tematic effects. 
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The data are remarkably good and show that the receiver-synthesizer is at least 
as good as any other system yet tried. Two such receiver systems will be required 
to obtain more definitive data. Considering that there are two multiplier chains from 
5 to 1400 MHz and two digital synthesizers in the system, we can only conclude that for 
averaging times to 1000 sec, this system is at least comparable in stability to the 
existing NAS 8-2604 maser system. We believe that the data point (taken twice) at 
1000 sec may be a new low bound for instability measurements between masers. 

In testing the maser, we found some real weaknesses in the magnetic shields. 

These shields are 6 years old and have been reused many times in different configura- 
tions of the maser. No doubt they have suffered work-hardening and should be reannealed 
or replaced. This problem is not apparent in the SM-1 maser, which has not been 
disassembled dozens of times, as has the present one. 

The maser is still equipped with a fused-silica cavity. Since the new CER-VIT 
cavity was not obtained in time for assembly in this maser, it will be tested in a 
subsequent contract. 

The packaging of the maser into a single cabinet is shown in Figures 7a, b, and c. 
Additional space in the cabinet is available for a digital synthesizer and a 5-MHz 
oscillator. The base of the cabinet has room for a battery supply that could operate 
the maser for about 4 hr. With the addition of the digital synthesizer, oscillator, 
and batteries, this maser system can be a completely self-contained source of highly 
stable signals at 5, 60, and 1400 MHz. 
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Figure 7a. Packaged maser. 
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Figure 7b. Packaged maser with power supply and front panel removed. 
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Figure 7c. Packaged maser with side panel removed. 
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INVENTION DISCLOSURE 


To the attention of the NASA Patent Department: 

The herein described invention is submitted in pursuance of 
Contract NSR-09-01 5-098 . 

The invention relates to the Atomic Hydrogen Maser Frequency 
Standard . 

The purpose of the invention is to provide a new type of cavity 
resonator for the maser. This cavity is made of three pieces of 
CER-VIT ©(a glass-ceramic material having an extremely low thermal 
coefficient of expansion) and is made tuneable in an extremely 
stable mechanical manner. 


Brief Description 

Cavity resonator mistuning causes a shift in the output fre- 
quency of the Hydrogen Maser proportional to the extent of the 
mistuning times the ratio of the cavity Q to the line Q. Normally 
this ratio is about 10" s so that a 10 HZ shift causes a 10“ 4 HZ 
shift (or 7 parts in 1 0 1 4 ) in the output frequency. Hence the 
mechanical and thermal stability of the cavity are of great im- 
portance to the overall stability of the maser. An extremely low 
thermal expansion material called CER-VIT having thermal expansion 
coefficient a <|0.5 x 1 0 ~ 7 1 and having excellent mechanical pro- 
perties is extremely suitable for the cavity. Heretofore this 
material has been used for the top and bottom plates and for the 
cylinder of an earlier cavity. The adjustment of frequency was 
done by inserting metallic parts and problems of their mechanical 
stability in view of the desired micro-inch positional tolerances, 
contributed to the instability of the cavity in both a thermal 
and mechanical manner. 

A new type of maser cavity has been designed and constructed 
using only three CER-VIT components. Figure 1 shows the cavity 
assembly that includes 1) a top plate with three integral stiffen- 
ing webs, 2) a cylinder with a three-fold interrupted helical ramp 
or thread at its bottom, and 3) the base which also contains the 
bottom plate. The cavity is shown resting on an aluminum disc. 
Three inclined ramp lugs on the base engage the cylinder inter- 
rupted thread. The cavity is tuned by rotating the cylinder and 
thereby raising or lowering it on the base. Figure 2 shows the 
cylinder and top plate up-ended. The web structure of the base 
can be seen through the transluscent CER-VIT. Note that the base 
legs are rigidly connected to the bottom plate. Figure 3 shows 
the maser bulb assembly fastened to the top plate. 
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The assembly shown in Figure 1 is secured by means of a 
hold-down cylinder that is fastened to the aluminum disc as the 
base of the cavity. A Belville spring assembly is used to pro- 
vide uniform downward pressure at the outer circumference of the 
top plate. The tension of the spring is set so that thermal 
movement of the aluminum hold-down cylinder relative to the CER- 
VIT does not cause an appreciable change in the compressive force 
on the cavity. Figure 4 is a sketch of the hold-down assembly. 

This new concept of an adjustable all CER-VIT cavity was 
conceived November 17th. First drawings were made by Thomas E. 
Hoffman on November 23, 1970. The cavity was built by Owens- 
Illinois and delivered June 17, 1971. We silvered the cavity 
and tested its electromagnetic properties July 14-15, 1971. 

This disclosure was first written June 1, 1971. 


This invention was explained to me by the above-identified 
inventor on i ? August. . 1971. 



by RT Vessot 

334 Ocean Ave. Marblehead, Mass. 



T. E. Hoffm'a'n 

18 Franklin St., Marblehead, Mass 



Big Rock Rd., Manchester, Mass. 
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6. OTHER ACTIVITIES AND PUBLICATIONS BEYOND THE SCOPE 

OF THIS CONTRACT 


Concurrent with work directly related to the contract were other activities closely 
connected with application of the maser to very long -baseline interferometry and to 
measurement of the gravitational redshift by means of a direct maser clock in a 
satellite. 

A. The small maser originally built for NASA under Contract NAS 8-2604 was 
refurbished and extensively tested. With NASA’s permission, the maser has been 
used in very long-baseline interferometry experiments at the following sites: 

Agassiz Observatory, Harvard, Mass. 

Haystack Observatory, Westford, Mass. 

Naval Research Laboratory, Maryland Point, Md. (two trips) 

National Radio Astronomy Observatory, Greenbank, W. Va. 

Jet Propulsion Laboratory Deep-Space Station DSS-62, Robledo, Spain 

The methods employed to adapt existing stations to use the maser are described in 
a reprint from Radio Science at the end of this section. 

B. The continuing interest of the principal investigator and coinvestigator 
associated with the present NASA contract, as well as the interest of NASA’s Physics 
and Astronomy Department, has led to further study in the feasibility of testing the 
Einstein equivalence principle by measuring the gravitational redshift of the earth. 
This experiment would make use of the extraordinary frequency stability of the 
hydrogen maser in a vehicle that moves upward and downward by being either in a 
highly eccentric orbit or, as suggested by Dr. Nancy Roman of NASA, in a rocket 
probe where the clock is launched vertically to about 10, 000 n.mi. and allowed to 
impact into the ocean. 
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Vessot and Levine presented "Measurement of the Gravitational Red shift Using a 
Clock in an Orbiting Satellite" at the conference on Experimental Tests of Gravitation 
Theories at the California Institute of Technology on November 11, 1970. This was sub- 
sequently published in the proceedings of that conference and forms part of the present 
report. 


C. Further activity by the principal investigator and coinvestigator in the field 
of frequency standards has resulted in the following papers, which are also appended 
as part of this report: 

Recent Developments Affecting the Hydrogen Maser as a Frequency Standard, 
published in NBS Special Publication 343 Precision Measurement and Funda- 
mental Constants, edited by D. N. Langenberg and B. N. Taylor, August 
1971. 

Characterization of Frequency Stability, published in IEEE Transactions on 
Instrumentation and Measurement, Vol. IM-20, No. 2, May 1971; also as 
NBS Technical Note 394, issued October 1970. 
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7. CONCLUSION 


The contract has been very productive, and we have achieved the objectives set 
forth in the original work statement. The establishment of a laboratory at SAO to 
continue the research and development of hydrogen masers for space and ground appli- 
cations has ensured that activities begun in 1960 at Varian Associates would be con- 
tinued in an unbroken manner. SAO’s affiliation with Harvard University and the fact 
that the Principal Investigator and Coinvestigator are Associates of HCO have led to 
stimulative interactions with the Physics and Applied Science departments of Harvard 
University. In particular, we have enjoyed a productive and interesting liaison with 
Professor Ramsey and his graduate students; we have given several informal seminars 
to this group and have coauthored articles with them. 

In addition to the interaction of ideas, there has been interaction at the hardware 
level. Two surplus hydrogen masers in very decrepit condition have been loaned by 
SAO to the Maser Lab at Harvard, where they have been rebuilt and put into service 
as frequency references. One of these was used in the flexible-bulb experiments. We 
have freely borrowed equipment and components, saving both time and money. During 
these years of stringent support for research, this interaction has been of great value 
to the programs both at SAO and at Harvard. 

SAO and HCO have an ongoing need for maser frequency standards to conduct Very 
Long-Baseline Interferometry (VLBI) experiments. This has led us to become involved 
with the reality of making the masers work in the field and in applications where the 
stability and high spectral purity of the maser signal are pressed to their limits. 

The small maser, developed under contract NAS 8-2604, has been the workhorse 
VLBI instrument. It has also served as a lab standard and as a transfer standard to 
bring together the Harvard, SAO, and NBS (hydrogen) masers for comparison with the 
NBS embodiment of the cesium frequency. Much has been learned about making a 
useful instrument, and this knowledge has led us to the redesign of a ground-based 
maser that incorporates many of the improvements developed or discovered in the 
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contract now completed. The redesign has been done under contract with a consortium 
of five separate organizations, NRL (A), NRL (B), National Radio Astronomy Observa- 
tory (NRAO), MIT, and the Swedish National Science Foundation. For the redshift 
experiment mentioned below, two additional masers will be required for NASA. 

Our facilities at SAO have been expanded to meet the demands of these contracts. 
The new package is a completely self-contained frequency standard requiring input 
power at 110 v AC and 28 v DC. 

During the course of our contract, we have been called on to refurbish or rebuild 
masers at MIT, NRAO, NRL, and the U. S. Naval Observatory. This work, performed 
under other contracts, has resulted in keeping a number of badly needed instruments 
operational so that VLBI and other experiments could be conducted. The experience 
gained in such field activity has been valuable, since we can review after many years 
the design faults of masers put into the field as long as 7 years ago. 

We have now begun a redirection of our efforts (under NASA Contract NAS 8-27969) 
toward the goal of performing an experiment of potentially great significance to physics 
and astronomy. We plan to apply the hydrogen maser technology to measure the gravi- 
tational redshift due to the earth's mass and to test by direct measurement the equiv- 
alence principle. 

We are grateful to the National Aeronautics and Space Administration for the 
support during the last 3 years that has allowed us to continue our work on masers at 
SAO. We look forward to a continuing successful involvement with NASA as we proceed 
toward our new objectives. 
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A hydrogen maser was installed at the Harvard College Observatory-Smithsonian 84-foot 
radio telescope in Harvard, Massachusetts, for very long baseline interferometry (VLBI) ex- 
periments in October, November, and December of 1969. The maser is a compact unit (42 
by 22 by 22 inches) of relatively low power consumption, designed specifically for field use; 
it is readily transportable. A precision 5-MHz crystal oscillator that is phase locked to the 
maser provides all frequency and time references for the VLBI receivers, recorders, and clocks. 
The maser was tuned against a rubidium standard by a rapid flux-switching technique; maser 
frequency resettability was estimated to be approximately ±3 X 10 13 with an over-all accuracy 
of approximately ±1 X 10 -12 . Loran-C was continuously monitored during the VLBI experi- 
ments, and the results of the Loran-C maser comparisons are included. 


INTRODUCTION 

The hydrogen maser is generally regarded as the 
ideal time and frequency standard for VLBI; it has a 
unique combination of superior short-term phase sta- 
bility and freedom from drift. Further, the maser is a 
primary standard with unsurpassed resettability of 
frequency. In the past, the size, weight, and general 
awkwardness of the equipment have limited the use 
of the hydrogen maser to fixed installations. For ob- 
servatories with only periodic or infrequent need for 
the phase stability of the maser standard, the invest- 
ment required for the maser and its supporting elec- 
tronics could not be justified. 

The Smithsonian Astrophysical Observatory 
(SAO) hydrogen maser described in this paper was 
designed to be easily transportable and rugged 
enough to withstand the normal rigors of being moved 
from one station to another. It was transported by 
truck from Cambridge to Agassiz Station, Harvard, 
Massachusetts, for use during two VLBI experiments 
from October 14 to 15 and November 22 to De- 
cember 14, 1969. The maser was installed at Agassiz, 
tuned against a rubidium-cell reference, and op- 
erated as the master time and frequency standard for 
the observatory for the duration of the observations. 
At the completion of the experiment, the maser was 
returned to Cambridge for use in other experiments. 

HYDROGEN-MASER OSCILLATOR 

The SAO hydrogen maser originated from a pro- 
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gram to develop a clock for orbital applications for a 
precise test of the principle of equivalence [Kleppner 
et al, 1970]. A second program was funded for con- 
struction of a small ground hydrogen maser incorpo- 
rating the advanced technology of the space maser. 

The design features of the SAO hydrogen maser 
are described in detail in an earlier paper [Vessot 
et al., 1968], Figure 1 is a diagram of the major as- 
semblies of the maser. It should be noted that the 
design philosophy of this maser emphasizes good 
long-term frequency stability without the necessity 
for continuous automatic tuning. The cavity resonator 
is silver-coated Cer-Vit (Owens-Illinois trade name 
for a low-expansion coefficient partly vitrified glass) 
thermally compensated and enclosed in a two-stage 
oven. Provisions have been made for manual or 
automatic tuning that would use only a stable ex- 
ternal reference oscillator and a digital-period coun- 
ter. The tuning procedure, usually required only 
when the maser is moved to a new location, is given 
in detail in a later section. 

The short-term performance of the SAO hydrogen 
maser is summarized in Figure 2, which is a conven- 
tional «t(t) plot. The rms frequency stability as a 
function of averaging time is shown for quartz-crys- 
tal oscillators, rubidium and cesium standards, and 
the maser. The same data are shown in Figure 3 in a 
form that might be more useful for VLBI experimen- 
ters. The rms phase deviation as a function of obser- 
vation time (normalized at 1 GHz) is shown for the 
four types of frequency standards. 


1287 


1288 


LEVINE AND VESSOT 



Fig. 1. A cut away section of the SAO hydrogen maser 
showing: 1, R F discharge tube; 2, ion pump; 3, hexapole 
magnet; 4, magnetic shield; 5, solenoid; 6, tuning loop; 7, 
low frequency transition coil; 8, output coupling loop; 9, 
cavity structure; 10, quartz storage bulb; and 11, oven heat- 
ers. This unit, which constitutes the maser oscillator, is 
mounted with the electronic support systems in a steel 
cabinet 42 x 22 X 22 inches. 


MASER CLOCK SYSTEM 

The hydrogen maser is an active oscillator with an 
output of approximately —95 dbm at 1420 MHz. A 
5-MHz crystal oscillator is phase locked to the maser 
signal to provide standard frequencies at levels useful 
to the clock system. A simple dual-conversion re- 
ceiver is used to generate the phase locking signals; a 
synthesizer is also used to set the time scale. 

The entire maser clock system is shown in block 
diagram in Figure 4. The maser output at 1420-405 
MHz is fed to a balanced crystal diode mixer where 
it is hetrodyned with the 1400-MHz output of a 5- 
to 1400-MHz (or X280) multiplier. The first IF 
signal at 20.405 MHz is amplified by 40 db and then 
heterodyned with the 20-MHz output of the x4 mul- 
tiplier. The first IF signal at 20.405 MHz is ampli- 
fied by 40 db and then heterodyned with the 20-MHz 


output of the x4 multiplier. The second IF signal, 
at 405 kHz, is amplified by approximately 50 db 
and fed to one input of a balanced phase detector. 
The reference signal for the phase detector is sup- 
plied by a Hewlett-Packard 5 102 A digital synthe- 
sizer. The synthesizer is usually set at 405,794.3 
Hz to approximate the UTC time scale. The out- 
put of the phase detector is used to phase lock 
a precision 5-MHz quartz-crystal oscillator through 
a very high-gain operational integrator. Double in- 
tegration is incorporated into the phase-locked 
loop so that the phase tracking error (or velocity 
error) is very nearly zero. The time constant 
of the integrator is selected so that the closed loop 
response is critically damped and the closed loop 
bandwidth is approximately 10 Hz. The output of 
the phase-locked 5-MHz oscillator is used to drive 
the x280 and x4 multipliers and the digital synthe- 
sizer so that the entire maser clock system is phase 
coherent. A 1-MHz signal, derived from internal 
dividers in the crystal oscillator, is used to drive the 
VLBI timing system. 

The Agassiz VLBI installation includes a Hewlett- 
Packard 5065A rubidium standard. This instrument 
provides the 5-MHz quartz-crystal oscillator used in 
the hydrogen maser clock system. The rubidium-cell 
reference loop in the 5065A is usually disabled and, 
as shown in Figure 4, the quartz oscillator is slaved 
to the maser. Using the rubidium standard in this 
fashion has a number of advantages. The cost of an 
additional quartz oscillator is eliminated, a stable 
backup standard is always available on-line, and the 
tuning of the maser is simplified. The change from 
hydrogen maser operation to the rubidium standard 
requires only the turn of a switch to transfer the elec- 
tronic-frequency control of the quartz oscillator from 
the external maser circuitry to the internal rubidium 
cell. The amplitude of the clock signals at 5 MHz, 1 
MHz, and 100 kHz is not affected by the changes in 
the phase lock circuitry, nor are any interconnection 
changes required. 


TUNING THE HYDROGEN MASER 


The oscillation frequency of a hydrogen maser is 
‘pulled’ by the maser cavity resonator. The amount of 
pulling is a function of cavity detuning and the line 
width [Vanier and Vessot, 1964] of the atomic res- 
onance : 
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Fig. 2. The rms frequency stability of the 
hydrogen maser as a function of observation 
time. The characteristics of typical quartz- 
crystal, rubidium, and cesium standards are 
shown for comparison. 



where 
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is maser oscillation frequency. 

is center frequency of the atomic resonance 

line. 

is center frequency of the cavity resonance. 

is quality factor of the cavity. 

is mean velocity of the hydrogen atoms. 

is first Bohr orbit radius. 

is Planck’s constant divided by 2ir. 

is Bohr magneton. 

is volume of the cavity. 

is volume of the storage bulb. 



Fig. 3. The data of Figure 2 are replotted to express di- 
rectly the phase stability of the four types of frequency 
standards as a function of observation time. All the data 
have been normalized to 1 GHz [from Cutler and Vessot, 
1968]. 


r) is ratio of average electromagnetic field 
energy density in the storage bulb to average 
energy density over the entire cavity. 

Av, is atomic resonance line width. 

If the cavity is tuned to a frequency v c0 such that 

Vco = vh{1 - [0 .29 Wh V./QWv V t ))~ l (2) 

then the term in brackets in equation 1 vanishes and 

v = v„ (3) 

The maser oscillates at the center of the atomic reso- 
nance line. 

The cavity is tuned to v c0 by modulating the flux 
of atoms entering the storage bulb. It can be shown 
that the spin-exchange contribution to the line width 
is directly proportional to the flux, and (referring to 
equation 1), when the cavity is at the tuning point 
f c0 , then v is independent of line width and therefore 
also independent of the total hydrogen flux. 



Fig. 4. The maser clock system in block diagram form. 
The rubidium cell in the Hewlett-Packard 5065A is dis- 
abled, and only the internal quartz-crystal oscillator is used. 
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The maser can be tuned by using a frequency ref- 
erence much less stable than that of the maser itself. 
The technique used is switching of the hydrogen flux 
periodically from a very low level, near threshold, to 
a very high level and looking for corresponding 
changes in the oscillator frequency of the maser. The 
flux level must be modulated as rapidly as possible to 
minimize the phase-noise contribution from the refer- 
ence oscillator. The phase-noise power spectral den- 
sity of quartz-crystal oscillators, rubidium standards, 
and cesium standards decreases as the offset from the 
carrier frequency is increased. 

The effect of linear frequency drift in the reference 
standard can be eliminated by using an algorithm 
devised by L. C. Cutler (private communication, 
1968). The period of the beat between the maser 
and the reference oscillator is measured an odd num- 
ber of times, alternately at low and at high flux. The 
measurements at low flux are assigned a negative sign 
and the measurements taken at high flux are assigned 
a positive sign. The first and last measurements are 
divided by two, and all the measurements are then 
summed algebraically. This residue, 

Nr = -L a l +-?■) + g (4) 

Z k-2 

is a measure of the cavity detuning. It can be readily 
shown that the value of the residue is unaffected by 
linear frequency drift of the reference. 

The maser at Agassiz was tuned against the HP 
5065A rubidium standard by the procedure outlined 
in the previous paragraph. Control of the atomic hy- 
drogen flux into the maser storage bulb is achieved 
by varying the RF power level in the hydrogen dis- 
sociator. The power level can be switched between 
two preset values by a manual control on the front 
panel or by a logic signal from an external controller. 
The output of the phase detector (Figure 4) was dis- 
connected from the 5065A rubidium standard and 
temporarily connected to a digital frequency counter 
set to read 10-period averages. The 5065 A was slaved 
to its internal rubidium-cell reference to provide maxi- 
mum stability. The frequency of the digital synthe- 
sizer was offset by 1.2 Hz to yield a beat period of 
0.833 sec at the phase-detector output for a 10-period 
measurement interval of 8.333 sec. 

Figure 5 shows the results of a tuning run at Agas- 
siz taken just before the VLBI observations. The 
flux switching was done manually with the front-panel 
control on the maser. The data, eleven 8.3-sec meas- 
urements for each point on the curve, were re- 
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Fig. 5. Data from a tuning run of 
the hydrogen maser using the rubid- 
ium standard as a reference. 


corded by a digital printer. A desk calculator was 
then used to calculate N R from equation 4. (The en- 
tire process has since been automated by the addition 
of an index register and a reversible counter.) 

The precision of this tuning procedure is estimated 
to be approximately ±7 X 10* 13 . That is, the maser 
can be reset to the same oscillation frequency with a 
standard deviation of about one millihertz. The actual 
oscillation frequency, the mean of a large number of 
such measurements, is a function of the bulb wall 
shift, the second-order Doppler shift in the bulb, and 
the magnetic field. These effects can be measured or 
calculated; they are discussed in the next section. 

MASER ACCURACY AND LORAN-C 
The wall shift of the SAO hydrogen maser has been 
directly measured against the Harvard University 
reference maser, the SAO reference maser, and a 
maser at the National Bureau of Standards. The 
second-order Doppler was calculated for the known 
temperature of the storage bulb, and the magnetic 
field was estimated from the measured frequency 
of the Zeeman transition. The oscillation frequency 
of the maser [Hellwig et al., 1970] is determined 
as follows: 

Unperturbed hydrogen fre- 
quency in terms of the 

cesium second 1, 420, 405, 751. 768 Hz ± 0.003 

Conversion for Al to UTC 42.612 

Second-order Doppler —0.053 

Measured wall shift —0.020 

Magnetic field shift +0.001 

Total (v H in UTC) 


1, 420, 405, 794. 298 Hz ± 0.003 
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Fig. 6. The VLBI timing system shown controlled by the hydrogen-maser clock 
system. The Loran-C capability permits continuous monitoring of the epoch. 


The digital-synthesizer resolution was limited to 0.1 
Hz; therefore, the closest synthesizer setting was 
405,794.3 Hz. We thus introduced a known ‘error’ 
of —0.002 Hz; the VLBI timing-system clocks should 
have run 1.4 parts in 10~ 12 fast with respect to 
Loran-C. 

Loran-C timing pulses were continuously moni- 
tored at Agassiz from November 24 to December 15, 
1969. The Agassiz VLBI timing system gained 6 /tsec 
with respect to Loran-C, whereas the United States 
Naval Observatory (USNO) daily phase values pre- 
dict a gain of 1.6 psec. The difference of 4.4 pscc for 
a 21-day period corresponds to a local clock 2.42 
X 10 -12 high in frequency. Since it is known that the 
Agassiz clock was 1 .4 x 1 0 -12 low, owing to the gran- 
ularity of the synthesizer, the total average discrepancy 
between the maser and Loran-C was 

2.42 X 10- 12 + 1.4 x 10- 12 , or 3.8 X Hh 12 

The Loran-C transmissions^ are monitored and 
steered by the USNO, which establishes and main- 
tains a coordinated time scale based on the average of 
a large number of clock systems in this country and 
abroad. The unperturbed hydrogen frequency shown 
in an earlier paragraph was measured against the en- 
semble of cesium clocks at the National Bureau of 
Standards (NBS) in Boulder, Colorado. The discrep- 
ancy of 3.8 parts in 10~ 12 observed in this experiment 
lies well within the combined error limits of the 
Loran-C time scale and the precision of the hydrogen- 
cesium intercomparison at NBS, Boulder. 


VLBI TIMING SYSTEM 

The hydrogen maser clock provides the 1-MHz 
signal that drives the VLBI timing system shown in 
Figure 6. The SAO digital clock accumulates and dis- 
plays the time of day with a resolution of 1 /xsec. The 
digital clock also generates a 1-pps signal for the SAO 
digital delay generator. The delay generator is used 
to preset the Loran-C propagation delay and to de- 
termine the time differential between a portable clock 
and the system clock. 

Time of day is set into the VLBI timing system by 
means of either a portable quartz-crystal clock or 
Loran-C transmissions. The digital clock is designed 
so that it can be directly reset to the 1-pps portable- 
clock signal. Setting the digital clock to the Loran-C 
epoch is accomplished visually by means of an oscillo- 
scope triggered by the Loran rate generator. 

The ‘demand’ inputs of the timing system freeze 
the clock display at the time of a demand input pulse. 
The Loran rate generator, the delayed 1 pps, a VLF 
1-pps signal, and the VLBI event encoder can be used 
to ‘demand’ the clock. (The complete VLBI timing 
system is described in the article by R. D. Michelini 
in this issue.) 

CONCLUSION 

For VLBI observations of more than a few minutes 
or at wavelengths less than 5 cm, the hydrogen maser 
is the only satisfactory frequency standard. The SAO 
hydrogen maser has demonstrated that it can be used 
routinely and reliably for VLBI. The ease with which 
the maser can be transported makes feasible the fabri- 
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cation of complete traveling back ends for radio tele- 
scopes so that maximum utilization of existing ob- 
servatories and baselines can be achieved. 
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Measurement of the Gravitational Redshift Using a 
Clock in an Orbiting Satellite* 

R. F. C. Vessot and M. W. Levine 
Smithsonian Astrophysical Observatory 
and 

Harvard College Observatory 


I. Introduction 


The proposed experiment uses a hydrogen- 
maser clock in a satellite to measure the gravita- 
tional effect on time scales with an accuracy sub- 
stantially higher than has ever been used before. 
This is a test of the principle of equivalence, which 
asserts that there is no way of distinguishing lo- 
cally between a gravitational acceleration and an 
oppositely directed mechanical acceleration. This 
principle, first stated by Einstein (Ref. 1) as a 
generalization of the observed proportionality be- 
tween gravitational mass and inertial mass, can 
be justified only by experiments. Experiments 
were carried out by Isaac Newton using pendula of 
various materials; more recently, the proportion- 
ality has been tested to a few parts in 10^ by Roll, 
et al. (Ref. 2) using a highly refined Eotvos 
balance. 

The principle of equivalence affects radiation 
and manifests itself in the gravitational redshift, 
where a source radiating in a gr avitational field 
will appear to be shifted in frequency by a frac- 
tional amount; 


Af _ A<(> 

f 2 
c 


where A<f> is the gravitational potential difference 
between the source and the observer. 


A number of proposals have been made to test 
the equivalence principle for clocks (Refs. 3-7). 
Pound and Rebka (Ref. 8) and Pound and Snider 
(Ref. 9) used the extremely narrow linewidth of 
Fe^ radiation and absorption due to the Mossbauer 
effect in a series of experiments over a 75-ft (ap- 
proximately 23-m) vertical distance. Their re- 
sults confirmed the prediction of the- equivalence 
principle to 1 part in 100. 


The advances in space-flight technology over 
the last decade and the availability of atomic oscil- 
lators with frequency stabilities better than 1 part 
in 10^ make possible a much more sensitive test 
of the principle of equivalence applied to clocks. 
The experiment we propose is a direct test of the 
relation with an accuracy of 1 part in 10^ between 
the rates of proper clocks located at substantially 
different gravitational potentials. 


A discussion of this experiment has been re- 
ported recently by Kleppner, Vessot, and Ramsey 
(Ref. 10). 


The expression describing the fractional shift 
in frequency of a s ate llite- borne oscillator ob- 
served from the earth is given by 


This work was supported in part by contract NSR 09-015-098 from the National Aeronautics and Space 
Administration. 


54 



s 


1 ? H-'. 2 ) 


( 1 ) 


where the velocity is measured in terms of an 
inertial frame whose origin is at the center of the 
earth and whose axes are aimed at the "fixed" 
stars, and where (<j> s - 4> e ) is the gravitational po- 
tential difference between the satellite and the 
earth's surface. If we neglect the effect of the 
earth's rotation (or locate our ground station at 
the north or south pole) and if we consider the 
earth to be spherical, we obtain the time average 
of the redshift: 



where a is the semimajor axis of the orbit, and 
GM e /c^ R e = 6. 94 X 10“ 10 i s the total redshift 
that would result if the satellite were at rest at a 
very large distance from the earth. The value of 
<S(£> is shown plotted in terms of orbital radius in 
Fig. 1. It is seen that with a 24-h orbit, we obtain 
77% of the total effect due to the earth's gravity. 
The value of <(Sq^ for these orbits is 5.37 X 10*10. 

If the orbit is eccentric, S will vary periodi- 
cally, and we can describe the variations as a 
gravity- induced frequency modulation of the oscil- 
lator in the satellite. We can use two separate 
properties of the oscillator — the precision* and 
its stability! — by means of an eccentric orbit. 

Assuming a spherical earth and no other per- 
turbations, the time average value of S is inde- 



Fig. 1. Redshift and orbital period plotted versus 
semimajor axis in units of the earth radii 



Fig. 2. Extremes of redshift versus eccentricity 
(The upper and lower branches give the redshift 
at apogee and perigee, respectively. ) 


pendent of eccentricity and depends only on the 
orbital period. The average frequency of the 
satellite oscillator can be related to its previously 
determined value on the ground with an accuracy 
that is limited by the precision of the oscillator. 

On the other hand, the extent of modulation of the 
oscillator frequency due to excursions back and 
forth in the earth's gravitational potential can be 
determined with an accuracy that depends on the 
oscillator's stability. This stability, as will be 
seen later, depends on the averaging time. The 
applicable averaging time will depend on the period 
of the modulation. As in other modulated sys- 
tems, the accuracy of the determination can be 
improved by repeating the measurement over many 
cycles. 

For a 24-h orbit, the extremes in the value of 
S are plotted in Fig. 2 as a function of orbital 
eccentricity. 

There are limits imposed on the eccentricity 
of the orbit by the desirability of keeping the satel- 
lite constantly in view of the ground station so that 
the telemetry system will operate with no inter- 
ruption of carrier phase. In Fig. 3 (from Ref. 

10), we show the minimum angle of elevation h m j n 
as a function of eccentricity of the orbit for a 24-h 
orbit with inclination 28. 5°. The value of h m i n 
observed from earth stations at the equator and at 
latitude 20° are shown. On the same figure, the 
value of the diurnal variation is plotted versus 
eccentricity. Since it is desirable to keep the 
minimum elevation angle more than 15° above .the 
horizon, we will look into the characteristics of 


*By precision, we mean the ability of the oscillator to be independently restarted to oscillate at a pre- 
viously determined frequency. 

+ The stability is defined here as the two- sample variance or the Allan variance (see Ref. 17). 
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Fig. 3. Variation in redshift and minimum 
angle of elevation for 24-h orbit inclined 
at 28. 5° as seen from an earth station: 
(a) 20° (b) 0° 


an experiment performed with a 24-h orbit with 
eccentricity 0. 52 and a ground station at latitude 
20°. Under these conditions, we can expect the 
shift to vary from 3. 09 X 10" ^ to 6. 08 X 10“ 
with a 24-h period. The behavior with time of the 
redshift is plotted in Fig. 4 (from Ref. 10). Shown 
in the same figure as a first-order doppler effect 
for a carrier frequency of 2. 5 GHz. 

II. The Telemetry System 
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Fig. 5. Concept of doppler- canceling telemetry 
system 


to remove the doppler shift and to extract the de- 
sired information. This system is shown in its 
conceptual form in Fig. 5, and its operation is 
described below. 

Two signals are used, a clock signal from the 
satellite and a tracking signal that is transmitted 
from the ground, received at the satellite, and 
coherently transmitted back to the ground station. 
The ratio of the received signal frequency f" to 
the transmitted frequency fo of the tracking signal 
is given by 


From Fig. 4, we see that the frequency shift 
we want to measure is very small compared to the 
doppler shift of the telemetry carrier. However, 
by use of a phase-coherent system, it is possible 


1 ' P 2 ' C 1 2 \ J 1 ~ ^3 * *23 
1 ‘ P 1 ' £ 12/ V " P 2 ' £ 23, 



Fig. 4. Redshift and doppler shift versus time 
for a 24-h orbit with eccentricity 0. 52 and in- 
clination 28. 5° (obtained from an earth sta- 
tion at 20° latitude using a 2.5-GHz carrier) 


As discussed in Ref. 10, this expression results 
from a direct application of the special theory of 
relativity, and its accuracy is sufficient for the 
purpose of this experiment. The symbols are 
explained by Fig. 5. 

The clock signal is transmitted from the sat- 
ellite at a frequency fo that is proper to the clock 
in the satellite. This signal is received by the 
ground station at a frequency f{). The relationship 
between fp and fg is given by 


f0 /l + (2Vc )-P 2 \/l-P ; 


f 0 \l + (2* 3 /c 2 )- \l - p 2 • 


This expression results directly from the principle 
of equivalence. It also results from the general 
theory of relativity when terms in v^/c^ and 4>/c^ 
are taken to first order. 


When applied to the simple case where the 
radial accelerations are small and when the 
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Fig. 6 . Schematic block diagram of a symmetrical doppler- canceling telemetry system (Data are 

taken both in the satellite and on the ground. ) 


propagation path to and from the satellite is the 
same, Eqs. (3) and (4) will give us an expression 
for S. For this case, where « ^ 3 = - t 12 and Pi = 
P 3 and 6 S and & e are defined as (5^ • and 

P j • e 1 2 . respectively, we obtain 


Af *. f n 
out 0 , 

= l. 
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n f o 


2 \f 
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S[1 - (6 s - & )] 


+ 6 e (6 - 6 JX [l -(6 - 6 >] 


( 5 ) 


The first term is the shift S multiplied by a small 
doppler correction; the second is a residual dop- 
pler effect resulting from the ground- station 
velocity. 


The effects of changes in radial velocity of 
the satellite and the change in position of the 
ground station due to the earth's rotation depend 
on the choice of orbit and the location of the ground 
station. The correction terms that result from 
these effects can be computed from knowledge of 
the orbit and the station location. 


The telemetry system described above cancels 
the first-order doppler by dividing the frequency 


shift in the tracking (go- return) signal by two and 
subtracts this average shift from the frequency of 
the one-way clock signal. Since asymmetrical 
propagation can occur in the up-down paths, the 
doppler correction of the frequency of the clock 
signal may be in error. To avoid this problem, 
it is very desirable to operate the telemetry sys- 
tem in a symmetrical manner, taking the redshift 
data both on the ground and in the satellite. This 
is accomplished by the addition of a transponder 
on the ground that coherently retransmits the clock 
signal from the satellite. In this way, the satellite 
can obtain information for doppler cancellation and 
correct the received frequency of the ground- 
station signal. This signal is controlled by the 
ground-based c lock. The symmetrical system is 
shown schematically in Fig. 6 . All frequencies 
are derived from the "proper" frequency fp, which 
in this illustration is taken as 100 MHz. Doppler 
shifts are denoted by & and relativistic-plus- 
doppler shifts by 6 '. 


The up- link to the satellite transponder is at 
2100 MHz, and the down-link is at 2100 X N/M = 
2280. 5 MHz (N = 240, M = 221). The satellite 
clock signal down- link is at 2100 Q/P = 2290 MHz 
(Q = 241, P = 221), and the ground- station trans- 
ponder frequency is at 2100 (Q/P) (M/N) = 

2108.7 MHz. These carrier frequencies should 
be sufficiently close together to avoid dispersion 
problems in the atmosphere and ionosphere. 
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This system will allow continuous tracking of 
the maser clocks by a continuous monitoring of 
the phase of the two signals transmitted from each 
of the ground and satellite clocks. Data will be 
recorded both in the satellite and on the ground and 
will be compared to the known orbital parameters 
of the satellite. To make a redshift comparison 
compatible with the clock stability of 7 parts in 
10 15, the following requirements on orbit informa- 
tion must be met: 

Radial distance Ar - 240 m 

GM known to 1 part in 10^ 

Vector velocity of Av ~ 60 cm/s 

satellite 

Radial acceleration r ~ 10 4 m/s^ 

Note that much of the tracking information is 
available from the doppler measurements made by 
the system itself and that the limitations to the 
accuracy of the redshift determination are not 
likely to result from the tracking system but will 
depend on the stability of the clock. 

III. Corrections to the Redshift 

The data obtained will contain contributions 
from other gravitational effects. These have been 
estimated in Ref. 10. The significant ones re- 
ported here apply to the elliptic orbit described 
earlier, with the ground station at the equator: 

(1) Quadrupole moment of 
the earth 

(a) Correction to the 
ground- station 

redshift = 3. 76 X 10'*’ 

(b) Correction to the 
satellite redshift 

at the perigee of , 4 

the € = 0. 52 orbit =9.0X10 

(2) Gravitational effect of 
the sun 

(a) Correction to the 

ground- station ,7 

redshift = 1.81 X 10 

(b) Correction to the , - 

satellite redshift = 1.83 X 10 

(3) Gravitational effect of 
the moon 

(a) Correction to the 
ground- station 

redshift = 3.91 X 10 

(b) Correction to the , . 

satellite redshift = 3. 28 X 10 

(4) Effect of earth tides 

Correction to the 

ground- station , - 

redshift = 1 X 10 


Corrections for these effects can be computed 
from information relating to the geopotential of the 
earth and earth-body tides. This information is 
available to far greater accuracy than we need in 
Gaposchkin and Lambeck (Ref. 11) and Kohnlein 
(Ref. 12). 

IV. Hydrogen- Maser Clock System 

The atomic hydrogen maser has been under 
continuous development since its invention in I960 
by Goldenberg et al. (Ref. 13) and Kleppner et al. 
(Ref. 14). At present, when used to control the 
phase of a crystal oscillator, it provides the 
highest stability of any known system. The char- 
acteristics of this system are compared to typical 
characteristics of other oscillators in Fig. 7 
(Vessot, Ref. 15). Note that statistical descrip- 
tions of the Rb- gas- cell- controlled crystal oscil- 
lators and of the crystal oscillator have been ob- 
tained by removing linear drift of frequency. 

The high stability of the maser results from 
the storage-bulb technique that is used to confine 
ground- state hydrogen atoms in the upper hyperfine 
level, allowing them to interact with RF resonance 
radiation within a cavity for intervals of 1 s or 
longer. During this time, the atoms are coherently 
stimulated and deliver their energy to the cavity, 
sustaining the level of RF field if a sufficient flux 
of atoms enters the bulb. 

The device is a self-oscillator at 1 . 4 GHz 
with a linewidth given by iTTqm where Tp is total 
effective storage time of the bulb. Normally, the 
oscillator Q is in the order of 109. 

Figure 8 is a schematic diagram of the maser 
and shows the energy levels of atomic hydrogen as 
a function of magnetic field strength. A cutaway 
view of the compact maser for space applications 
is shown in Fig. 9- In brief, the maser operates 
as follows. Molecular hydrogen is fed to an RF 
discharge dissociator, and the dissociated mole- 
cules are collimated into a beam of hydrogen 
atoms directed into a vacuum system along the 
axis of a hexapole magnet. Atoms in the F = 1, 

Mp = 0, and 1 states are focused into the storage 
bulb through a small, we 11- collimated hole. Atoms 
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Fig. 7. Stability of hydrogen maser as a 
function of averaging time 
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Fig. 8. Schematic diagram of maser and ground- 
state energy levels of atomic hydrogen 


in the other two states are deflected away from 
the beam axis. The storage bulb confines the atom 
in an unperturbed way in the in-phase region of 
the RF magnetic field of a circular E-mode reso- 
nant cavity. Upper state, Mp = 0, atoms are 
stimulated to radiate their energy during their stay 
in the bulb and leave the bulb eventually to be 
scavenged by the ion pump. When the power avail- 
able from the atoms in the bulb exceeds the cavity 
losses and the power coupled out, the maser 
oscillates. 

Normally the output power is about -97 dBm 
and the signal is used to control the phase of a 
crystal oscillator. In this way, signal outputs at 
useful frequencies can be generated. More de- 
tailed descriptions of the maser and its associated 
electronics are available in the listed references. 

The stability of the maser is limited by ther- 
mal noise whose frequency components lie within 
the linewidth of the oscillator and by additive 
noise that competes with the output signal in the 
bandwidth of the phase-lock system. The expres- 
sion relating the maser stability to the noise within 
the linewidth of the oscillator and within the band- 
width of the receiver is given by Cutler and Searle 
(Ref. 16): 
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Fig. 9. Cutaway view of NASA maser 


cavity frequency variations have been found to have 
a f' * spectral density, evident from the flattening 
out of the <r( t) plot of Fig. 7 (Ref. 17). 


ct(t) 


kT / Fu B Q e 

2pl 2 2 0 
\V Q c 



1/2 


( 6 ) 


Here, ct(t) is the fractional rms frequency devia- 
tion averaged over a time interval t, kT is the 
thermal noise power per unit bandwidth, ug is the 
receiver half-bandwidth (single tuned bandpass), 

P is the power delivered to the cavity by the 
atoms, Q e is the atomic line Q, and Q c and Q e 
are the loaded cavity and external cavity Q, 
re spectively . 

The ultimate limit to the maser stability for 
long time intervals (t > 10^ s) is imposed by 
variations in the resonance frequency of the cavity 
that are usually of thermal origin. The cavity 
mistuning "pulls" the output frequency by the ratio 
of the cavity Q to the line Q times the amount of 
the mistuning. Normally this ratio is about 10" 5, 
so a 10 Hz change in the cavity causes a 10*^ Hz 
change in the output (or 7 parts in 10 14). The 


It is possible to improve long-term stability 
at the expense of short-term stability by increas- 
ing the storage time of the bulb. This is accom- 
panied by a reduction in the output power from the 
maser and a corresponding worsening in the short- 
term stability that goes as P” The tradeoff of 

short- versus long-term stability is shown in Fig. 
10 for several values of the total bulb relaxation 
rateY, which is the reciprocal of the effective 
bulb- storage time. 

Several processes limit the storage time in 
the maser bulb (see References). Recent studies 
at Harvard and at the Smithsonian Astrophysical 
Observatory (SAO) on specially prepared polytetra- 
fluorethylene (PTFE- Teflon) surfaces inside the 
storage bulbs have shown that 10^ or more colli- 
sions with the bulb surface can occur before the 
atom loses phase with the RF field in the cavity. 
With bulb sizes now in use, this means that storage 
times approaching 10 s should be possible; and 
with the usual type of cavity and thermal controls, 
values of the rms frequency stability will approach 
1 part in 10^5 for averaging times of 10^ s. 
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Fig. 10. Maser stability showing influence of total bulb relaxation 
rate at long and short averaging times 


Experiments to verify this expectation are in prog- 
ress at SAO. 

The accuracy of the output frequency from the 
hydrogen maser is limited by the ability to deter- 
mine (1) the magnetic field in the bulb, (2) the 
second-order doppler shift of the atoms in the bulb 
(which is related to the bulb temperature), (3) the 
amount of systematic residual mistuning, and (4) 
the systematic average phase shift per collision of 
the stored atoms with the walls of the bulb (wall 
shift). These effects are discussed in detail in 
Refs. 18, 19, and 20 and will be reviewed only 
briefly here. 

A. Magnetic Field 

The average field in the bulb is determined by 
observing the quenching of the maser output when 
a resonance transition is caused between the F = 1 
magnetic sublevels. These transitions have a fre- 
quency dependence of 1.4 MHz/gauss, and using 
this resonance, the field can be measured to better 
than ±3 X 10" ^ gauss. Since the field dependence 
of the 0 — 0 transition is described by = 

2750 the mean square value of the field 

over the bulb volume is required. Because gradi- 
ents in the magnetic field can be observed through 
the effect of magnetic quenching of the output sig- 
nal, we can trim the magnetic fie Id , u|ing separate 
coils in the solenoid so that <^H^>£ - is very 

small. An upper bond on the estimate of inaccur- 
acy due to this is a few parts in 10^. 

B. Second- Order Doppler 

The second-order doppler frequency shift is 
due to the thermal motion of the atoms in the stor- 
age bulb. Since each atom in the bulb makes a 


large number of collisions, the atoms are in ther- 
mal equilibrium with the bulb. The expression 
describing the shift is 


A *T = - v 0-~2 = - 1 ' 9557 T 
2mc 


where vp = 1. 4204 GHz, k is Boltzmann's con- 
stant, T is in kelvins, c is the velocity of light, 
and m is the mass of the hydrogen atom. Mea- 
surement of the bulb temperature to an accuracy 
of ±0. 25°K at 320°K allows a determination of 
Av'y/vq to about 3 parts in 10^. 

C. Cavity Mistuning 

Cavity mistuning, as described earlier, will 
"pull" the output frequency, and a systematically 
mistuned cavity can be a source of inaccuracy in 
the output frequency of the maser. A part of the 
NASA contract that lead to the development of the 
maser shown in Fig. 9 consisted of the invention 
and development of an automatic cavity tuning 
servo system. 

The exact expression for the effect of cavity 
mistuning (Ref. 21) is 

Au s = (— ff^'Qc - Bt1/2 )^ v T 


where B is a constant that depends on the bulb and 
the cavity and on the cavity-bulb geometry, T is 
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in kelvins, and V -y is the total relaxation rate of 
the atoms in the bulb. 

1 /2 

The BT term results from shifts in fre- 
quency due to hydrogen-hydrogen collisions. This 
effect is described by Bender (Ref. 22). Crampton 
et al. (Ref. 21) have shown that the shift is pro- 
portional to the linewidth, and by use of the spin- 
exchange quenching tuning method described be- 
low, the shift is eliminated. 

Spin- exchange collisions occur among the 
atoms in the storage bulb and produce pressure- 
dependent relaxation among the energy levels in 
the hydrogen atom. This relaxation process, neg- 
ligible under normal operating conditions, can be 
used to' modulate the total relaxation rate 'Y'p. If 
the quantity in parentheses in Eq. (7) is not zero, 
there will be a modulation in Av s whose magnitude 
and phase are in proportion to the magnitude and 
sign of the bracketed quantity. The automatic 
cavity tuning servo consists of a synchronous de- 
tection system to observe the variations in the out- 
put frequency due to A» s and to correct the cavity 
resonance frequency so that Av s is nulled. 

The cavity itself is made of CER-VIT and is 
compensated, by use of re-entrant posts of an alloy 
with a high-expansion coefficient, to remove the 
variations in dielectric loading of the fused silica 
bulb due to the temperature dependence of its di- 
electric constant. The temperature coefficient of 
the cavity resonance frequency has been reduced 
to 100 Hz/"C. With thermal control within limits 
of ±1 X 10"^°C, the resultant output- frequency 
variations due to the cavity can be maintained be- 
low 1 part in 10 . 

The fractional frequency stability between two 
masers operating with automatic tuning systems 
has been measured to be 2 X 10" ^ for averaging 
times of 4 X 10^ (1/2 day) (Ref. 23). Thus, for a 
single maser, the stability can be estimated to be 
Z/\jl X 10" or 7 parts in 10^5 f or the above 
averaging time. 

The accuracy of the automatic -tuning tech- 
nique depends on the type of tuning reference oscil- 
lator used and on the time allowed for th e servo to 
operate. When a second hydrogen maser is em- 
ployed as a tuning reference, less than 1 h is 
necessary for an accuracy of about 1 part in 10^ 
to be accomplished. Since no systematic effects 
on the accuracy of the maser output frequency 
have thus far been attributable to the tuning pro- 
cedure, the accuracy of the tuning should approach 
the stability of the tuning, although this remains to 
be proved experimentally. 

D. Wallshift 


difficult to reproduce the same texture of the 
Teflon from one bulb to another. During the last 
year, there has been considerable activity investi- 
gating the properties of Teflon surfaces. Zitzewitz 
at Harvard (Ref. 24) found that the phase shift per 
collision for FEP Teflon* could be made to cross 
through zero from negative to positive at about 
1 10°C, a function of increasing wall-coating tem- 
perature. Vessot and Levine (Ref. 25) at SAO 
found that for high T molecular-weight PTFE sur- 
faces sintered and rapidly quenched, a zero wall- 
shift temperature of 83°C could be obtained that 
was independent of the surface-to- volume ratio of 
the bulb. Brenner (Ref. 26) suggested using a de- 
formable bulb to vary, by a known amount, the 
surface-to- volume ratio of the bulb, thus allowing 
a determination of the wall shift of the coating in 
a particular bulb. This method was successfully 
demonstrated by Debely (Ref. 27), using a cylin- 
drical bulb whose conical end, made of thin flexible 
Teflon and stretched inside or outside the cylinder, 
would determine two easily calculable bulb volumes 
with the same surface. 


Since the zero wallshift temperature is inde- 
pendent of the collision rate, knowledge of the 
collision rate is no longer required if the correct 
temperature can be determined. Recently, we 
suggested that the flexible bulb also be used to vary 
the collision rate in order to determine the tem- 
perature at which there is no collision- rate de- 
pendence on the wallshift. We have described a 
servosystem that accomplishes this automatically 
(Vessot and Levine, Ref. 28). 

Recent measurements by SAO and the National 
Bureau of Standards of PTFE surfaces applied in 
a reasonably standardized way have given agree- 
ment to 5 parts in lO 1 ^ in the maser frequency 
(Ref. 29). 

E. Summary of Present Maser Capability 


The intrinsic accuracy of the maser can be 
estimated from the rms of the several accuracy 
limitations : 
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The precision of resettability can be esti- 
mated from the rms of the parameters that vary 
with time : 


The systematic average phase shift per wall 
collision has been the chief limitation in the in- 
trinsic accuracy of the hydrogen maser. Varia- 
tions in the wallshift determined at several labora- 
tories have led to differences of some 5 parts in 
1012 in the maser output frequency. Since the 
wallshift depends on the collision rate of the atoms 
in the bulb, the texture of the Teflon wall surface 
is an important consideration in a determination 
of the surface- to- volume ratio of the bulb. It is 
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Fig. 11. Photograph of NASA maser developed 
for satellite use 
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The fractional frequency stability of the maser 
is 7 X 10-15 for averaging time intervals from 
100 to 4 X 10 4 s. 


Fig. 12. Lightweight maser based on existing 
design concepts 


We see that for an unmodulated gravitation 
experiment, for example a circular orbit, the ac- 
curacy of the measurement will depend on the 
magnitude of the shift (see Fig. 1) and the pre- 
cision of resettability of the clock. Prior to 
launch, the satellite clock and ground clock would 
have been compared. For those clocks not in- 
cluded in this preflight calibration the experiment 
accuracy is limited by the accuracy of the clocks. 

The advantage of using an eccentric orbit re- 
sults from the excellent stability of the maser and 
the fact that the measurement can be repeated 
many times, as long as the satellite is operating. 
We must, however, be certain that there are no 
systematic frequency shifts that are synchronized 
with the orbital period. 

V. Maser Configurations, Weight, and Power 

Estimate s 

Since I960, a continuous development of hy- 
drogen masers has been under way with NASA's 
support. In 1966, NASA supported the develop- 
ment of a small, relatively lightweight maser sys- 
tem for use in spacecraft. The maser oscillator 
of this system is shown in Fig. 11; the automatic 
tuner system, phase-lock synthesizer, thermal 
controls, pressure controls, and RF dissociator 


have been successfully operated in breadboard 
form. The cavity and bulb assembly was vibration 
tested in the sine sweep mode; it successfully sur- 
vived such tests in all three axes. The complete 
structure shown in Fig. 11 weighs about 112 kg 
and requires about 30 W to operate; it is about 81 
cm high and 48 cm in diameter. 

The vacuum of space can profitably be used in 
place of the ion pump to scavenge the expended 
hydrogen. A considerable reduction in weight and 
power will result from use of multilayer, reflect- 
ing insulation under vacuum to replace the micro- 
balloon foam insulation currently used to surround 
the bell jar and ovens. The total weight of the 
maser can be reduced to about 34 kg, and its 
power to about 10 W. A sketch of this lightweight 
maser is shown in Fig. 12. The design concepts 
of this maser are similar to those of the self- 
contained unit. Some changes in operation are re- 
quired, however. Before launch, the maser will 
be operated by use of a titanium sublimation pump 
located inside the hatch cover. This system will 
enable ground comparisons and preflight testing 
to be conducted. 

During launch, the maser will be shut down; 
later, in orbit, it will be restarted and the cavity 
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retuned. Several hours of operation will be possi- 
ble before the titanium evaporated just before 
launch is used up, and this will allow preliminary 
system checks. The hatch will be opened to the 
void of space as soon as the satellite outgassing in 
the direction of the pump reaches an acceptably 
low level. 

Contamination streaming back into the maser 
will be kept to a low level by baffles, and the con- 
tamination entering the bulb, if it has not been ad- 
sorbed on the metal baffles, is very unlikely to 
attach itself to the Teflon- coated interior of the 
bulb, which will be maintained at 83°C. In view 
of the even slight possibility of wallshift changes 
due to contamination, a test of the continued accu- 
racy of the maser can be made by comparing the 
average frequency over several orbits early in the 
experiment with the average frequency at a later 
time. The stability of the maser, over 4 X 10^ s 
averaging time intervals, would be negligibly af- 
fected by bulb-contamination processes that slowly 
and monotonically could shift the maser output fre- 
quency. The data from redshift modulation will 
continue to be available. 

Other, more compact, maser configurations 
are possible with the use of dielectric or capacitive 
loading of the RF cavity to reduce its size. Ex- 
periments are in progress at SAO v on a spherical, 
fused- silica, dielectrically loaded cavity 20. 3 cm 
in diameter with an inside diameter of 12. 7 cm. 

The features of spherical geometry and the integral 
bulb make a very rugged, compact unit. The ori- 
ginal proposal to NASA for a satellite-borne maser 
(Ref. 30) included this cavity. Development and 
procurement of this cavity proceeded until 1966, 
when the emphasis on light-weight and low-power 
consumption was changed to that of pressing for 
the greatest possible precision of resettability and 
stability. The cavity was taken over by the 
USAECOM and successfully used in a maser (Ref. 
31). The recent improvements on wall coatings 
and the desire for a compact lightweight device 
suitable for a wide range of vehicles bring this de- 
sign back into prominence. Figure 13 shows the 
maser with a spherical cavity. 

VI. Conclusion 

The proposed gravitational redshift experi- 
ment described in this paper makes use of tech- 
nology currently available. Since 1964, when the 
original proposal was submitted to NASA, the 
hydrogen-maser clock has been under continuous 
development for its use in spacecraft, and the 
feasibility of a lightweight, low-power device has 
been demonstrated. 

Plans for the experiment have evolved con- 
siderably. In. 1968, the George C. Marshall Space 
Flight Center completed a Preliminary Program- 
Development Plan for a Hydrogen- Maser Relativity 
Satellite, which envisioned a 2700- lb (approxi- 
mately 1000 kg) satellite containing two masers 
requiring a total of 1200 W. This satellite would 
be launched by a Titan IIIC booster into a synch- 
ronous, circular orbit. Its objective would be the 
measurement of the redshift with an accuracy of 
about 500 parts per million. 
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Fig. 13. Hydrogen maser using a dielectrically 
loaded spherical cavity 


Since 1968, the precision and stability of the 
maser clock have been substantially improved, 
largely as a result of continuing NASA support for 
maser development. Also, the use of an eccentric 
orbit, where the satellite scans back and forth 
through large differences in the earth's gravita- 
tional potential, has made it possible to utilize 
the maser's excellent frequency stability to good 
advantage. The present concept of the satellite 
clock experiment will give much better accuracy 
at substantially lower cost. A total satellite 
weight of 500 lb (187 kg), consisting of the maser, 
its electronics, a transmitter- transponder , bat- 
teries, and solar panels, is now possible. The 
thrust- augmented Thor- Delta vehicle can launch 
this payload into an eccentric orbit. 

The previously mentioned limits to eccentri- 
city, imposed by the desire to keep the satellite 
constantly in sight of the ground station, can be 
relaxed, thus giving a larger gravitational modu- 
lation. Since a satellite with a 24-h period in a 
highly eccentric orbit cannot be observed both at 
perigee and at apogee from a single ground sta- 
tion, a different period should be used, allowing 
a ground station to track all phases of the orbit 
over a period of time, though not continuously. A 
further advantage to this highly eccentric orbit is 
that more than one ground station can independently 
participate in the experiment. The limit on eccen- 
tricity is probably imposed by the requirement that 
the satellite be in sight of a Darticular station long 
enough to allow 1 h or longer of phase tracking of 
the doppler- cor rected carrier. At the 2. 5-GHz 
carrier frequency, assuming about 10° of phase 
resolution, the system should allow each measure- 
ment to be made to the limits of the clock stability. 
At present, there are several tracking stations al- 
ready equipped with hydrogen masers that, with a 
modest amount of additional equipment, can be 
used to track the satellite maser signals. Studies 
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on weight, power, orbits, and booster capability- 
are continuing. 

We believe that an accurate test of the equiva- 
lence principle for clocks in space is a necessary 
phase of the continuing space program and that this 
test, successfully performed, will not only ad- 
vance our knowledge of science but also open new 
fields of technology that will be of value to our ex- 
ploration of the universe. 
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Recent Developments Affecting the Hydrogen Maser as a 

Frequency Standard * 
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The processes limiting the accuracy and stability of hydrogen-maser frequency-standard 
systems are described. Recent developments of the maser systems in the following areas are dis- 
cussed: 

(1) The temperature dependence of wall shifts. 

(2) The use of the zero wall-shift temperature to eliminate the wall shift. 

(3) The deformable storage bulb technique for varying the wall-collision rate. 

(4) The large-bulb maser. 

(5) Electronic phase-lock systems and cavity -tuning servos. 

Key words: Accurate frequency standards; accurate time standards; frequency stability; maser 
oscillators; oscillator phase locks; quantum electronics. 


1. Introduction 

1.1. Description of the Maser 

The atomic hydrogen maser was invented in 1960 
by Goldenberg, Kleppner, and Ramsey [1]. Since 
this time, the maser has been under continual devel- 
opment as a frequency standard. At present, when 
used to control a crystal oscillator with a phase-lock 
system, it has been found to exceed that of all other 
known systems for sampling periods ranging up to 
at least 10 4 s. 

Figure la shows a schematic diagram of the maser; 
figure lb, the energy levels of atomic hydrogen; and 
figure lc, a cutaway drawing of a compact maser 
oscillator packaged in a cylindrical form 20 inch in 
diameter and 32 inch in length. The hydrogen maser 
is the first successful maser to use a confining tech- 
nique to keep a constantly replenished group of state- 
selected atoms in a region of rf field. The stored 
atoms, interacting with their rf magnetic field at their 
resonance frequency, are stimulated to emit con- 
tinuous radiation. When this radiation is strong 
enough to sustain a coherent rf magnetic field at a 
high enough level, the system becomes self-oscillating. 
In practice, the atoms are stored in a bottle located 
so that it occupies a volume of rf magnetic field of 
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Administration, the National Bureau of Standards, the National Science 
Foundation, and the Office of Naval Research. 
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nearly constant phase and direction in a resonant 
cavity. 

Limitations in the stability and accuracy of the 
maser arise from several causes. The fundamental 
stability limitations are imposed by the oscillator 
line width, thermal noise within the line width, and 
thermal noise within the bandwidth of the equipment 
receiving the maser signal. Further stability limita- 
tions arise from systematic effects that perturb the 
oscillator frequency. 

Accuracy limitations arise from incomplete knowl- 
edge of the magnitude of four contributions to fre- 

? uency shift: (1) the magnetic field offset Av m , 
2) the second-order Doppler shift Av T , (3) the 
effect of cavity mistunings Av„ and (4) the effect 
of atomic-hydrogen collisions with the walls of the 
storage bottle (wall shift) Av w . These shifts, when 
added to the unperturbed hydrogen hyperfine fre- 
quency separation t>o, give p 0 ut = ^ 0 + hv T -\- 

&v B -\-Av a , the maser output frequency. 

At present, the most uncertain correction is the 
wall-shift effect. 

1.2. Stability Limitations 

Stability is restricted by thermal noise. The expres- 
sion giving the effect of noise-frequency components 
within the line width of the oscillator and in the 
bandwidth of the receiver is given by Cutler and 
Searle [2]: 
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Figure 1. (a) Schematic diagram of the hydrogen maser-, 
(6) energy-level diagram of atomic hydrogen. 


Here, <r(r) is the fractional rms frequency deviation 
averaged over a time interval t; JtT is the thermal 
noise power per unit bandwidth; c ob is the receiver 
half-bandwidth (single-tuned band-pass character- 
istic) ; P is the power delivered by the atoms to the 
cavity; Qt is the atomic line Q given by Qi = oio/272r; 
Q c l is the loaded cavity Q; and Q e is the external 
cavity Q. 



Generally, for r<30 s, the first term in eq (1) 
dominates and the stability function cr(r) goes as r _1 . 
At longer periods, the effect due to the second term, 
resulting from thermal noise within the oscillator 
line, should dominate. The t ~ 112 behavior thus pre- 
dicted has not yet been observed because it has been 
obscured by systematic effects that have leveled off 
<t(t) at about 7X 10“ 15 . Figure 2 shows the two- 
sample variance [3] as a function of averaging time r 
for the maser. The lower curve shows the theoretical 
behavior (as described in eq 1) that applies to the 
maser and to the system used to obtain the data. 

The single most important aspect of the storage- 
bulb maser technique is the ability to store an 
aggregate of atoms with an average velocity very 
close to zero in the standing wave of a cavity so that 
they can interact for an extended period of time, At, 
with the in-phase components of the rf magnetic field 
in the cavity. The narrow resonance line width Av 
results from the Heisenberg uncertainty principle, 
AvAt^I. In the oscillating maser, the narrow reso- 
nance line width is electrically equivalent to an 
oscillator having very high Q\ this Q, designated Qi 
in eq (1), is proportional to the storage time At. 
The net effective storage time, equal to the reciprocal 
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Figure 2. Fractional frequency stability as a function 
of averaging time. 
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of the net relaxation time yrr, results from a series of 
relaxation processes that go on in the storage bulb: 

y2T = 7& + 72i0 + 72s« + 72m, (2) 

where yb = i)A a K/4Vb is the escape rate of the atoms 
from the bulb. Here v is the average velocity of the 
atoms, Vb is the bulb volume, A a is the aperture area, 
and K is the collimator kappa factor for the entrance- 
exit hole [43. Also, y? w = (v/\)uh r is the wall-collision 
relaxation rate, X is the mean distance between colli- 
sions (=4/3 R in the case of a spherical bulb), and 
w 2 , is the probability of a dephasing collision at the 
bulb wall. Chemical recombination of hydrogen or 
chemical interaction with the wall destroys the phase 
coherence of the atomic oscillating dipole moment 
with the rf magnetic field and is a dephasing collision. 

The spin-exchange dephasing collision rate is 
T 2 «« = nVrcr. It depends on the relative atomic velocity 
v„ the spin-exchange cross section <r, and the density 
n of the atoms in the bulb. The term y 2m is the 
dephasing effect of the dc magnetic-field gradients 
in the bulb. The oscillation frequency of the coupled 
electron and proton magnetic-dipole moments has, 
in the "field independent” transition, a second-order 
magnetic-field dependence Av m = 2750 H 2 , where H is 
the magnetic field. The atom will gain or lose phase 
with the rf field as it moves in the bulb if the dc 
magnetic field varies from place to place within the 

Magnetic gradients can also cause transitions 
among the F=1 magnetic sublevels at magnetic 
fields where the collision rate of the atoms is near 
the transition frequency of the sublevels. 

The strongest limitation on the long-term stability 
encountered thus far is the effect of cavity mistuning. 
Since the maser is an active oscillator, the cavity 
resonance frequency affects the output frequency in 
proportion to the ratio of the cavity Q to the line Q. 
This "pulling” effect is described approximately by 
the expression 

Al'out'^'APcavity (.Qc/Ql) ' ( 3 ) 

In practice, there is a lower limit on the cavity Q 
set by the requirement that the maser oscillate. An 
oscillation parameter q, given by Klepper et el. £4], 
relates the cavity-bulh geometry, the cavity Q, and 
other constants. The expression for q is 

_ vvAiy-nVc J. Aot 
^ 87T/io 2 76f) VbQ I ’ 

where 

a is the spin-exchange cross section for hydrogen, 

v r is the average relative velocity of the atoms, 

ft is Planck’s constant divided by 2ir, 

Vo is the Bohr magneton, 

V c and Vb are the cavity volume and bulb volume, 
respectively, 

Q is the loaded cavity quality factor, 

V is (Hz) bulb 2 /(H 2 }cavity, and 

I tot / 1 is the ratio of the total atomic flux to the 
flux in the desired state. 

Here, H z is the Z component of the rf magnetic field. 


For oscillation, the value of q must be less than 
0.172. 

The threshold flux for oscillation, neglecting spin- 
exchange processes, is 

I th = ft V c y2T 2 / 4 ttho 2 Qv- 

The power delivered by the beam to the cavity 1 is 

P=MM/thC-2 9 2 (///th) 2 +(l-3 g )(/// th )-l]. 

(4a) 

The maximum power is 

Pmax = (o>f, 2 V c y 2T 2 /64ir^Q v ) [1 - (6 /q) + (1/q 2 ) ]. 

(4b) 

The flux at maximum power is 

/™*p = (ft V c /4wno 2 Qv) 727’ 7 ( 1 - 3 q) /4q 2 ]. (4c) 

Normally, the loaded cavity Q is about 35,000 and 
the line Q is about 1.5X10 9 , making the ratio QJQi 
about 2 X 10 -5 . Thus, a shift of the cavity frequency 
of about 100 Hz will cause a change of 2 X 10 -3 Hz 
in the output or a fractional frequency shift of about 
1.5 X10- 12 . 

There are several possible approaches for reducing 
the effect of cavity "pulling.” One is to stabilize the 
cavity by improving the mechanical and thermal 
properties of the cavity. Another is to employ a 
cavity-tuning servosystem using a stable reference 
oscillator. This servo constantly maintains the cavity 
near its proper frequency by modulating the line Q 
and correlating the shifts in frequency with the 
modulation. The cavity is kept tuned by maintaining 
the condition that no correlation exist between the 
modulation and the frequency shifts. 

A third approach, involving a fundamental feature 
of the hydrogen maser, is to extend further the 
storage time of the atoms in the bulb. At this point, 
it should be recognized that in the self-oscillating 
maser, the reduction of the atomic line width is 
accompanied by a reduction in output power. This 
results from the y 2 t 2 term in the quantity J t h in eq 
(4). The stability at short time intervals is controlled 
by the first term in eq (1), which gives the r _1 
behavior. The power depends on yn 2 , and the line 
width on 72 r _1 . With Qi = o)/2yrr, the product of the 
first term in the brackets and kT/2P is adversely 
affected and the second term is unaffected. The line 
width reduction helps directly in the cavity-pulling 
expression. The effect of varying y 2 T on the long- and 
short-term stability is shown in figure 3. For this 
family of curves, the maser receiver noise figure and 
bandwidth are kept constant, the cavity bulb geom- 
etry and Q are kept constant, and only y 2 t is varied. 
The contributions from a cavity -mistuning perturba- 
tion having a particular 1// or "flicker” frequency 
spectral density are shown by the horizontal lines 
for c t(t) f2 J, the contributions from additive white 
phase noise are shown by the r _1 slope in the lines, 
and the constant limit set by thermal noise within 
the line is shown by a line with slope r _1/2 . 


1 If the wall-relaxation is due to recombination of atoms. 
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Figure 3. Effect of varying the relaxation rate yvr i n the maser: 
yi > yi >~n. 


1.3. Accuracy Limitations 

The following four effects cause a bias in the 
frequency of the maser oscillator. They are system- 
atic shifts, as distinguished from the random effects 
described above. 

A. The magnetic-field behavior of the F=l, mp = 
0— >F = 0, mp = 0 transition shown in figure 1 is 
described by 

A„ m = 2750 (H*) b , 

where Av m is in hertz when H is in oersteds. 

The squared magnetic field (IP )* is averaged over 
the volume of the bulb. To establish the value of 
Av m , the average value of H is determined by exter- 
nally causing Amp = ±l transitions among the F = 1 
state sublevels and observing the frequency v z at 
which the maser output is quenched. The average 
field H in the bulb is obtained by the expression 
(i/)i = Kz/1.4X10 6 . Since v z is easily measurable to 
better than 5 Hz, at the normal operating value of 
10 -3 Oe, the value of (H)b is known to better than 
±3X10~ 6 Oe and Av m can be determined to 1 part 
in 10 14 . If gradients are present in the bulb, it is not 
appropriate to use (H)b 2 in place of { H 2 )b ■ However, 
the effect of magnetic gradients can be determined 
and minimized by observing the quenching of the 
maser through the y m term previously discussed. 
The error in Av m resulting from taking (H)b 2 for 
(H 2 )b is about 1 part in I0 14 , giving an overall error 
in Av m of a few parts in 10 14 . 

B. The second-order Doppler frequency shift 
results from the thermal motion of the atoms within 
the bulb. Since the atoms on the average spend about 
0.5 s in the bulb and make in the order of I0 4 colli- 
sions with the walls, the temperature of the bulb is 
taken as the temperature of the atoms. The second- 
order Doppler shift is given by 

Avt— —i>D(3kT/2mc 2 ) = — 1.9557 X10 _4 T. 

Here, k is Boltzmann’s constant, T is in Kelvins, c is 
the velocity of light, and m is the mass of the hydro- 
gen atom. Since we can measure the temperature to 


better than 0.5 K, the probable error in the second- 
order Doppler shift is in the order of 5 parts in 10 14 . 

C. The effect of cavity mistuning on the output 
frequency results from the pulling effect described 
earlier. There is also a systematic shift due to atom- 
atom collisions Q5]. One can make use of the spin- 
exchange interatomic collisions to modulate the line 
Q of the maser for tuning purposes, as described 
earlier. When this is done, the output frequency is 
shifted by the pulling process and a systematic shift 
due to interatomic collisions. Both these processes 
are proportional to the line width and are described 
by tne expression for the mistuning shift: 

A^cavity n _ 13V2 vf>a 0 2 V c \ 72 T 

v 0 32 QWvVb) tt 

which is a more exact version of eq (3). 

If we modulate 727- by varying the flux into the 
bulb, we can make the quantity in the brackets zero 
by adjusting Av ca vity The density -independent out- 
put frequency (Ay, = 0) resulting from this tuning 
process in the maser effectively removes the collision 
shift [6]. 

The effect of the cavity shift on the accuracy of 
the maser frequency depends on the type of tuning 
system used. If the device is constantly being tuned 
by a servodevice, the value of Av, will have a zero 
mean value and the value of the stability function 
o-(r) will probably vary as t ~ 112 for large values of r. 

If the system is tuned and left alone, the accuracy 
of the output frequency will depend on the precision 
of the tuning and the stability of the cavity resonance 
frequency. 

D. The strongest limitation on the accuracy of the 
maser is due to the fact that the atom, upon collision 
with the walls, suffers an average phase shift per 
collision. The collision process can be briefly de- 
scribed by a Van der Waal’s potential function with 
the usual r~ 6 dipole-dipole attractive term and a r -12 
term due to exchange forces. During its approach to 
the wall, the atom is polarized and experiences a 
r -6 type of potential, gaining energy until it en- 


30 



counters the r~ 12 repulsion due to close-range electron- 
exchange forces. It is therefore polarized during the 
time it is in the r~ 6 potential, resulting in a Stark- 
effect shift in the energy levels and a reduction in 
the ground-state hyperfine separation. The oscillating 
atomic dipole is thus retarded in phase with respect 
to the rf field in the cavity during the r -6 portion of 
its travel. During the repulsive portion of the colli- 
sion, the electron cloud surrounding the proton can 
be visualized as being compressed, and during the 
time the r~ 12 forces are acting, the interaction energy 
between the proton and the electron is stronger, 
resulting in an increase in hyperfine separation. The 
oscillating atomic moment is thus advanced in phase 
during this part of the collision. It has been found Q7] 
that the either positive or negative phase shift can 
predominate and that it is possible to establish condi- 
tions where an average zero phase shift per collision 
can be obtained. 

The wall shift Av m is usually determined by coating 
a number of bulbs of different sizes that provide a 
range of wall-collision rates that is as large as pos- 
sible. The maser output frequency is then plotted as 
a function of collision rate. The value obtained on 
extrapolating to zero collision rate and correcting 
for magnetic and Doppler effects is taken as the 
unperturbed hydrogen frequency. The wall shift is 
related to the phase shift per collision as follows: 

Av vj =(v c /2tt)$(T), 

where v c , the average collision rate, equals 35/4 R for 
a spherical bulb of radius R. The average phase shift 
per collision $(T) is a function of temperature, and 
v= (HkT/tzm ) 1/2 is the average velocity. 

Comparison of wall shifts is most often made by 
means of the quantity 

K(T) = A v w D = 35<? ( T) /4tt, 
where D is the bulb diameter. 

2. Wall-Coating Materials and 
Techniques 

The chief limitation in the accuracy of the maser 
as a frequency standard has been due to the wall 
shift. Accurate measurement of the wall shift is not 
easy and has depended on the ability of the experi- 
menter to coat many bulbs in a reproducible manner. 
The value of the wall shift thus obtained can be 
applied with confidence only to the bulbs actually 
measured. Variations between laboratories in the 
measurement of this parameter have led to differ- 
ences of some 5 parts in 10 12 in the maser frequency. 

In most cases, the material that has been used for 
coating the bulbs has been the FEP-120 Teflon dis- 
persion made by duPont. This is a polytetrafluor- 
ethylene-hexafluoropropylene copolymer supplied in 
an aqueous dispersion; it is coated on the inside of 
the bulb, allowed to dry, and then sintered or fused 
to form a coating. Various coating processes have 
been used. However, all are essentially the same and 


involve burning off the wetting agent by heating in 
air and then heating to a temperature where the 
coating is fused to form a slick, clear coating. Adhe- 
sion to the inside of the bulb can be improved by 
roughing the fused silica surface with glass shot and 
abrasive grit until a ground-glass appearance iB 
obtained. 

Early measurements of the FEP material were 
made by Vanier and Vessot [83- The wall shift that 
they gave for this material and that pertaining to 
the particular set of bulbs used in the measurement 
was until recently generally accepted as an intrinsic 
property of the material. This conclusion was wrong, 
since differences in the surface texture of the material 
can cause variations in the wall shift from one bulb 
to another. 

Recently, there has been considerable activity in 
the study of the wall shift, using various types of 
Teflon bulb coatings. Zitzewitz [9J at Harvard found 
that Teflon showed a temperature-dependent shift 
that went from negative to positive, as shown in 
figure 4. Here, the data are given in terms of phase 
shift per collision versus temperature. The breaks in 
the curve, identified by a, f), and y, are due to transi- 
tions in the internal short-range ordering of the 
molecules analogous to those in glass. Data are shown 
for two types of Teflon, the FEP copolymer and poly- 
tetrafluoroethylene homopolymer (PTFE duPont 
Teflon 42). The zero phase shift per collision tem- 
perature is about 107 °C. 



Figure 4. Temperature dependence of phase shift per 
collision for PTFE and FEP Teflon. 
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At the Smithsonian Astrophysical Observatory 
(SAO), studies on PTFE Teflon have been made by 
Vessot and Levine £10] using Teflon 42 under 
various conditions of curing. In these experiments, 
bulbs having integral collimators and smooth interior 
surfaces were coated with a single layer of Teflon 
dispersion. The dispersion consisted of 50 percent 
Teflon 42, 40 percent water, and 10 percent Triton 
x-100 (a wetting agent made by Rohm and Haas). 
The coatings were drained and dried for several hours 
by a slight flow of dry nitrogen introduced into the 
bulb through the collimator. The bulbs were sintered 
for 30 min at 360 °C in a horizontal air-flow oven 
with a constant flow of dry oxygen entering the bulb. 
The bulbs were then removed from the oven and 
quickly cooled in a stream of room-temperature air. 
It was found that the zero wall-shift temperature was 
reduced to 83 °C and that the wall-relaxation rate 
appeared to be very low. Very good reproducibility 
of phase shift per collision as a function of tempera- 
ture was obtained for several bulbs. 

The temperature slope of the phase shift per colli- 
sion for the quickly cooled (quenched) PTFE was 
found to be 0.15±0.004 /irad/°C. When the bulb 
was annealed in oxygen at a rate of 0.05 °C/min 
from 340 to 295 °C, the temperature slope decreased 
to 0.09±0.01 /urad/°C and the intercept temperature 
went up to 142 °C. Zitzewitz gives values of 0.1294 
and 0.1248 for the slope and a zero wall-shift tem- 
perature of 107 °C. For a recent experiment, he re- 
ports confirmation of the lower (83 °C) temperature 
at which the phase shift per collision is zero for a 
quenched PTFE-coated bulb. 

With the same PTFE sample as that used by the 
SAO and Harvard groups, Hellwig at the National 
Bureau of Standards (NBS) made wall-shift meas- 
urements. His data taken at 24 °C gave a value of 
the wall-shift parameter K very close to the SAO 
value at the same temperature: 

Knbs = 0.528 Hz cm 

at 24 °C. 

ATsao =0.522 Hz cm 

To relate the values of the unperturbed hydrogen 
frequency among the laboratories at Harvard, SAO, 
and NBS and to give this frequency in terms of the 
l33 Cs hyperfine frequency that serves as the definition 
of the second, the SAO group undertook to make a 
traveling maser frequency -standard experiment. This 
work, reported by Hellwig et al. [11], gives the 
following new value for the unperturbed hydrogen 
hyperfine separation frequency in terms of the NBS 
cesium standard: 

no= 1 420 405 751.768±0.003 Hz. 

This value is in disagreement with that obtained by 
Vessot et al. [12] in 1965, who gave 

1-0=1 420 405 751.7864±0.0017 Hz. 

The reasons for the difference are difficult to assess. 
Retrospective corrections for possible wall-shift errors 
in hydrogen and for possible errors in the frequency 


of the cesium devices have failed to close the gap 
between the two numbers and their error estimates. 

In the present determination of the value of j> 0 » 
the agreement in the unperturbed hydrogen fre- 
quency among all three laboratories is within 
±0.002 Hz, or about ±1.4 parts in 10 12 . 

In summary, the PTFE wall-coating material 
appears to offer several advantages in performance 
over the FEP preparation: 

A. The molecular weight of PTFE is larger by a 
factor of about 10 than that of FEP. The number of 
Teflon molecular-chain end groups for a given volume 
of material is thus reduced. The presence of atoms 
other than fluorine on the walls will cause different 
values of average phase shift per collision. 

B. The relaxation rate associated with hydrogen 
collisions appears to be lower for PTFE than for 
FEP, especially between 50 and 100 °C. From pre- 
liminary measurements of relaxation rate, it appears 
possible to extend the storage time to about 5 s. This 
would substantially reduce the linewidth of the 
maser. The substance can be made highly crystalline, 
a property that reduces surface porosity and should 
further extend the storage time of the maser. 

C. Structure does not change near the tempera- 
ture at which the phase shift per collision is zero. 

D. From recent data determined by different 
experimenters, wall-shift measurements appear more 
reproducible with PTFE than with FEP. 

3. Elimination of the Wall Shift 

The existence of a reasonably low operating tem- 
perature at which the wall shift goes to zerO and at 
which the wall-relaxation rate is small makes possible 
a new method of measuring the wall shift. 

Since the wall shift can be independent of the 
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collision rate only when the average phase shift per 
collision is zero, it is not necessary to know the 
collision rate to obtain the zero shift temperature. 
One need only operate a maser with bulbs having 
different collision rates and determine the tempera- 
ture at which the output frequency is the same. By 
obtaining the output frequency at this temperature 
the only further requirement for eliminating the wall 
shift is that the wall-coating material be the same 
for all bulbs. 

Figure 5 shows wall-shift data taken at SAO [10] 
with a 4 in-diam bulb and a 7 in-diam bulb and 
plotted in terms of A v a /T ll K The division by T w 
removes the effect of the temperature dependence of 
the velocity of the atoms and makes the data linear 
and the intercept easier to compute. The intercept 
of the two resulting straight lines is well defined and 
gives a value of the unperturbed hydrogen frequency 
with an error (0.2dhl.l) X10~ 3 Hz when related to 
the value obtained in the Harvard-NBS-SAO com- 
parison reported earlier [11]. 

4. The Flexible Bulb 

To overcome the uncertainties associated with the 
use of several separate bulbs to determine the wall 
shift, Brenner [13] suggested the use of a single 
flexible storage bulb that would allow the collision 
rate to be changed by a known amount without 
changing the nature of the surface. Recently, Debely 



Figure 6. Variable-collision-rate storage bulb. 


[14] at Harvard conducted experiments on a maser 
using an open-end cylinder with a flexible cone at the 
other end that could be pushed into or pulled out of 
the cylinder to change the storage volume. This 
storage vessel is shown in figure 6. The movable part 
was made of a thin film of the FEP-120 copolymer 
that was also used to coat the bulb and the collima- 
tor. The ratio of maximum to minimum volume in 
this experiment was 1.299±0.004. By tuning the 
maser in each of the positions, Debely obtained the 
wall shift of a reference maser, with an uncertainty 
of ±0.0023 Hz. In his design, he took pains to see 
that the annular joints of the Teflon coating of the 
bulb to the movable cone did not move or expose a 
different surface when the cone was moved. 

The flexible-bulb technique can be used to estab- 
lish the temperature at which the phase shift per 
collision is zero. This removes the uncertainty in the 
constancy of the wall shift with time and also over- 
comes the problem of using separate bulbs. By 
systematically changing the volume of the bulb 
(retuning the maser cavity at each position) and 
comparing the output frequency in each position 
against a stable frequency reference, one can adjust 
the temperature so that no change is observed in the 
output frequency. This procedure could be made 
automatic [15] by using servotechniques to main- 
tain the proper temperature and promises to help 
greatly in ridding the maser of its most serious 
limitation on accuracy. 

5. Multiple-Region Maser With 
Reduced Wall Shift 

Since the wall-shift and relaxation rates in the 
maser depend on the collision rate of the atoms in 
the storage vessel, an important advantage can be 
gained by simply enlarging the storage bulb. How- 
ever, limitations on bulb diameter of the self-oscil- 
lating maser are imposed by the cavity Q and the 
filling factor combination Vt/V c that appears in the 
oscillating parameter q (see eq 4) . A new concept of 
the maser was introduced in 1968 by Uzgiris and 
Ramsey [16] to overcome this limitation. A sche- 
matic of the new maser is shown in figure 7. The 
storage vessel consists of a 60 in Teflon-coated- storage 
bulb connected to two rf interaction regions, each 
consisting of an rf cavity and a bulb, as in the usual 
maser. One cavity picks up the radiation from the 
atoms in its bulb, and this signal is amplified and 
sent to the other cavity to prestimulate the atoms 
in the other bulb. Since free molecule flow prevails 
in the three connected bulbs, the atoms can be put 
into a radiant state to emit at an enhanced rate in 
the low-level cavity. The loop is thus closed and the 
system will oscillate. The role of the amplifier is to 
make up for the lack of cavity Q and of the filling 
factor in the expression for q. The voltage gain R in 
the amplifier is given by the ratio of the rf magnetic 
field in the stimulating cavity to the rf field in the 
detecting cavity. The value of q for this maser is 
simply the value given earlier for the usual maser 
multiplied by (1/R) ( Vb/Vd), where Vb is the large 
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Figure 7 . Schematic of multiple-region hydrogen maser. 


bulb volume and Vo is the high-level cavity-bulb 
volume. The threshold flux f t h in this case is the 
usual value times a factor 1 /R(Vb/Vd) 2 - Uzgiris 
calculates that the minimum gain required is about 
70 dB and found that about 80 dB is required to 
make the maser oscillate. 

The cavity pulling effect is complicated by the 
fact that both cavity resonant frequencies and the 
phase of the amplifier response must be included in 
the pulling expression along with a very small flow- 
dependent term that is not line width dependent. 
Except that the small effect of this term (in the order 
of 1 part in 10 13 ) must be accounted for, the tuning 
procedure is similar to that ordinarily used. The ratio 
of high-flux line-width to low-flux line width is about 
1.7, as in conventional masers. 

With the bulb geometry used in the experiment, 
the confinement time is calculated to be 10 s (y6 = 
0.1). At low -flux operation, a value of 72 r of about 
0.33 s -1 was observed, and the data indicate an excess 
relaxation rate of about 0.2 s~ l over the sum of the 
spin-exchange and bulb escape rates. This excess is 
most likely due to wall relaxation and is estimated 
to correspond to a probability of 10 -4 per collision. 

The results on the wall shift confirm the prediction 
that it would be 10 times lower than that encountered 
in the usual maser operating with the same wall 
materials and at the same temperature. A value of the 
shift of (— 2.3±0.8) X10 -3 Hz was obtained. At 
present, the multiple -region maser with a large stor- 
age vessel is being modified to incorporate the de- 
formable bulb technique discussed in section 4 for 
measuring the wall shift of the maser coating that is 
actually in use. With this modification, the large-bulb 
maser should combine the advantage of a much 
smaller wall shift with the deformable technique for 
measuring that small shift. 

6. Wall Relaxation 

The storage-bulb technique depends on the ability 
of the bulb surface material to bounce the hydrogen 


atom without causing it to lose phase with the rf 
magnetic field in the cavity. Processes that lose the 
atom through chemical reactions or that change its 
atomic state and, thus, cause a loss of magnetization 
(Tj processes), and collision processes that cause the 
oscillating atom to lose phase coherence with the 
rf field (To processes) both contribute to the overall 
limitation on the storage time. 

Early work by Berg [17] at Harvard gave values 
of the wall-relaxation probability w r of about 1 X 10~ 6 
per collision for FEP Teflon. Most of the masers 
constructed since then have used the FEP copolymer; 
and later measurements of the relaxation on this 
material [18] give values of tv r of about 1 X 10 -4 per 
collision. The relaxation process has been found to 
be due mostly to recombination of the atomic 
hydrogen. 

In recent work on PTFE, the value of relaxation 
probability per collision is approximately 1 X 10~ 5 , 
and the slow annealing process described earlier has 
yielded bulbs that give even more optimistic data. 
The difficulty in these recent results lies in the 
estimation of the bulb storage time, 7 &. Most of these 
bulbs have storage times ~0.7 s (7&~1.4 s _1 ). The 
wall-relaxation rates have been estimated by meas- 
uring the exponential decay 72 t of pulsed radiation 
after shutting off the atomic-hydrogen source and 
calculating y w = yiT — yt- For a 7 in bulb, values of 
7„, = 0.1 ±0.1 s -1 have been obtained. The sensitivity 
of the measurement can be improved by using bulbs 
with very long storage times, thus making wall 
relaxation the dominant effect. Work along these 
lines is in progress at SAO to determine if the 
annealed surfaces offer any improvement in the 
relaxation rate. 

If the present data on PTFE are correct, it should 
be possible to extend the storage time by a factor of 
3 to 5 beyond the limits obtained by using FEP 
Teflon. Values of 727- = 0.25 s _1 or less should be 
possible in 7 in-diam bulbs. The benefits of this type 
of Teflon on the large bulb maser would also be very 
significant. 
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7. Electronic Systems 

7.1. Cavity-Tuning Systems 

a. Maser Cavity Resonator 

The hydrogen masers constructed to date have 
all used a resonator operating in the cylindrical 
TEon mode. This mode is convenient because, since 
only azimuthal wall currents flow, electrical contact 
between the cylinder and the end plates is not 
required. 

Silver-coated CER-VIT 2 is used for cavity con- 
struction at SAO [19], while solid aluminum is pre- 
ferred by the Goddard Space Flight Center [20]. 
The CER-VIT cavities are tuned by means of a 
voltage-variable capacitor mounted within the cav- 
ity and coupled to the rf magnetic field by a small 
loop. The aluminum cavities are thermally tuned; 
the relatively high temperature coefficient of expan- 
sion is used to advantage to permit tuning by small 
adjustments of cavity temperature. 

b. Cavity Tuning 

The oscillation frequency of the maser is pulled 
by tuning the resonator, as shown by eq (3). 

There are essentially two problems associated with 
cavity pulling: (1) the cavity must be set to the 
correct frequency at the initiation of maser opera- 
tion, and (2) the cavity must be held at this fre- 
quency indefinitely in the presence of environmental 
perturbations. The first problem is related to accu- 
racy; the second, to stability. The setting of the 
cavity frequency can be accomplished by manual 
tuning techniques or by either of the two automatic 
tuners described in the next section. A series of 
careful experiments indicates that error in maser 
frequency attributable to cavity tuning is ±7X10 -13 
in a 90 percent confidence interval [11]. 

Once the maser is initially tuned, there remains the 
problem of maintaining this condition. One approach 
is to build the resonator out of materials of ultra-low 
thermal expansion, such as CER-VIT, which has a 
coefficient of expansion less than 5X10 _8 /°C, and 
adequately control the temperature. Alternatively, 
the cavity can be servoed to resonance by continuous 
operation of an automatic tuning system. 

c. Automatic Tuning Systems and Techniques 

The hydrogen maser is tuned by modulating the 
line Q and adjusting the cavity resonator so that 
there are no correlated frequency shifts. Several 
automatic systems for timing masers have been 
devised; two of these are described below: 

A. Random-walk tuning [20]. The masers used 
in this system are equipped with beam shutters so 
that the hydrogen flux and therefore the line Q can 
be modulated mechanically. 

The maser frequency is compared to the reference 
oscillator, a stable crystal oscillator; a comparator 


1 Owens-Illinois trade name for a partly vitrified glass with a low expansion 
coefficient. 


circuit determines if the maser frequency was meas- 
ured either high or low with respect to the reference. 
The comparator generates a "plus” pulse if the maser 
is higher in frequency than the reference and a 
"minus” pulse if it is lower. The plus and minus 
pulses are added algebraically in a reversible counter. 
When the comparator pulses accumulate to either a 
positive or a negative preset limit, a pulse of the 
appropriate sign is generated. These pulses are ac- 
cumulated in a cavity-register circuit. A digital-to- 
analog converter translates the binary number stored 
in the cavity register to a cavity -control voltage that 
is used to servo the cavity to the correct frequency. 

This random-walk automatic tuner, when used 
with an extremely low-noise crystal reference oscil- 
lator (tr(r) =4X10 -13 ), will correct an initial maser- 
frequency error of 3 X 10 -12 to approximately 2 X 10~ 13 
in 20 h. The ultimate tuning error is claimed to be in 
the order of 5 X 10 -14 . 

B. Fast automatic tuning [19, 21]. The fast auto- 
matic tuner modulates the line Q by changing the rf 

f )ower level into the hydrogen dissociator. The flux 
evel must be modulated as rapidly as possible to 
minimize the phase-noise contribution from the 
reference oscillator. 

The effect of linear frequency drift in the reference 
can be eliminated by using an algorithm devised by 
L. S. Cutler (private communication, 1968). The 
period of the beat between the maser and the refer- 
ence is measured an odd number of times, alternately 
at low and at high flux. The measurements at low 
flux are assigned a negative sign; those at high flux, 
a positive sign. The first and the last measurements 
are divided by 2, and all the measurements are then 
summed algebraically. This residue, 

Nr = ~ J (“!+«>•)+ £ ( — l)*a*, 

k - 2 

is a measure of the cavity detuning and is independ- 
ent of linear frequency drifts in the reference oscil- 
lator. 

In practice, each beat measurement is approx- 
imately 6.6 s, and a complete cycle consists of 21 
measurements, 11 at low flux and 10 at high flux. 
The period measurements are taken on a reversible 
counter, counting down for the low-flux measure- 
ments and up for the high-flux. The addition of one- 
half the 1st and 21st measurements is accomplished 
automatically by dividing the clock rate by 2 for 
these two measurements. The binary number in the 
counter at the end of the cycle (the residue) is trans- 
lated by a digital -to-analog converter into a voltage 
that is used to correct the cavity frequency. 

The precision of this tuning procedure is approx- 
imately ±7X10 -13 after 2 to 3 h of operation. Nor- 
mally, the maser is initially tuned and then left 
unperturbed except for periodic tuning checks. Con- 
tinuous automatic tuning is usually not desired, 
because of the average reduction in fine Q required 
by the modulation process. A further difficulty in 
continuous tuning can arise from the changes in 
maser output-signal level with flux modulation. Un- 
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Figure 8. Phase-lock system. 


less extraordinary precautions are taken in the design 
of the phase-lock circuitry, significant phase shifts 
will accompany amplitude changes. These phase 
shifts will not affect the automatic tuner, because 
they occur between operating cycles of the auto- 
tuner; but they will degrade the performance of 
external clocks. 

Long-term frequency-stability data between two 
masers each operating with autotuning servos have 
been taken £22]. A value of <t(t ) = 2Xl0~ 14 for r — 
12 h has been obtained. 

7.2. Phase-Lock Systems 
a. Description 

The hydrogen maser is an active oscillator with an 
output of approximately —96 dBm at 1420.4 MHz. 
A crystal oscillator is phase-locked to the maser 
signal by a simple dual-conversion system, shown in 
figure 8. To provide standard frequencies at useful 
levels, the synthesizer is used to set the time scale. 

The output of the maser at 1420.4 MHz is hetero- 
dyned in a low-noise balanced mixer with the output 
of the frequency multiplier at 1400 MHz. The 20.4- 
MHz output of the first mixer is amplified and 
heterodyned in the second mixer with a 20-MHz 
signal derived from the crystal oscillator. The re- 
sulting 405-kHz signal is amplified and used as one 
input to the phase detector. The reference signal is 
provided by the digital frequency synthesizer, which 
is referenced to the quartz oscillator. The dc output 
of the phase detector is used to phase-lock the quartz - 
crystal oscillator through a high-gain operational 
integrator. Double integration is incorporated into 
the phase-lock loop so that the phase-tracking error 
is very nearly zero. The time constant of the inte- 
grator is selected such that the closed-loop response 
has a bandwidth of about 20 Hz and is critically 
damped. 


b. Recent Developments 

Before 1968, the frequency of the crystal oscillator 
shown in figure 8 was usually 1 or 5 MHz. Recently, 
a high-quality 100-MIlz quartz-crystal oscillator has 
been substituted for the lower frequency slave 
oscillator. 

The use of the 100-MIlz oscillator offers two 
important advantages. First, the additive phase 
noise at 1400 MHz is reduced by approximately 
26 dB for the same oscillator-drive level. This per- 
mits the use of a much wider phase-lock -loop band- 
width and improves the overall signal-to-noise ratio. 
Typically, the 10- to 20-Hz bandwidth used with the 
5-MHz slave oscillator can be increased to 100 to 
150 Hz for the 100-MHz quartz oscillator. Second, 
the multiplication ratio to 1400 MHz is reduced from 
280 to 14; this is of particular importance because the 
multiplier chain is a source of phase noise. 

Usually, the phase noise is very small at 1400 MHz 
since at this frequency the local oscillator is most 
directly compared with the maser. The phase of the 
crystal oscillator, however, is steered to make up for 
phase changes in the multiplier, and if the output 
of the crystal oscillator is used, phase noise will be 
observed. Generally, the spectral density of the 
output-frequency fluctuations exhibits a flicker of 
frequency (1//) behavior. It has been observed that 
for a typical X280 multiplier driven by a 5-MHz 
oscillator, the cr(r) function is unaffected at 1 s, but 
in the range from 1 to 100 s, there is degradation of 
stability. The high-frequency crystal oscillator sub- 
stantially improves the performance in this range, 
and the value of <r (100 s) is limited by the cavity 
stability. 

8. Conclusion 

Present stability performance of the maser fre- 
quency-standard system exceeds that of all other 
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known frequency standards; however, the accuracy 
of the device has justifiably been open to question, 
mostly owing to problems involving the wall shift. 
The improvements that can be brought about as a 
result of recent developments when engineered into 
operational devices will make a great deal of differ- 
ence in the acceptability of the hydrogen maser as 
both an accurate and a stable frequency standard. 
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Abstract — Consider a signal generator whose instantaneous time averages implied in the definitions are normally not available; 
output voltage Vtt) may be written as thus estimates for the two measures must be used. Estimates of 

S V (J) would be obtained from suitable averages either in the time 
F(0 — [ V 0 + e(0] sin [2irv 0 t + ys(0] domain or the frequency domain. An obvious estimate for a\{j) is 


where V 0 and are the nominal amplitude and frequency, res- 
pectively, of the output. Provided that t(t) and y>(t) = ( d<p/(dt ) 
are sufficiently small for all time t, one may define the fractional 
instantaneous frequency deviation from nominal by the relation 

_ <p(0 

2irv 0 

A proposed definition for the measure of frequency stability is 
the spectral density S„(/) of the function y(t) where the spectrum 
is considered to be one sided on a per hert 2 basis. 

An alternative definition for the measure of stability is the 
infinite time average of the sample variance of two adjacent averages 
of y(t); that is, if 


1 

Vt - ~ / y(i) dt 

T J t k 

where r is the averaging period, it+i = ti + T, k = 0, 1, 2 • • • , to is 
arbitrary, and T is the time interval between the beginnings of 
two successive measurements of average frequency, then the 
second measure of stability is 


<T, 


to 



where ( > denotes infinite time average and where T = r. 

In practice, data records are of finite length and the infinite 
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Parameters of the measuring system and estimating procedure 
are of critical importance in the specification of frequency stability. 
In practice, one should experimentally establish confidence limits 
for an estimate of frequency stability by repeated trials. 
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Glossary of Symbols 

Bias function for variances based 
on finite samples of a process 
with a power-law spectral density. 
(See [13].) 

A real constant defined by (70). 

Real constants. 

A real, deterministic function of 
time. 

Expected value of the squared 
second difference of x(t) with lag 
time r. See (80). 

Fourier frequency variable. 
High-frequency cutoff of an ideal- 
ized infinitely sharp cutoff low-pass 
filter. 

Low-frequency cutoff of an ideal- 
ized infinitely sharp cutoff, high- 
pass filter. 

A real function of time. 

Positive real coefficient of /" in a 
power series expansion of the spec- 
tral density of the function y(jt). 
Integers, often a dummy index of 
summation. 

Positive integer giving the number 
of cycles averaged. 

Positive integer giving the number 
of data points used in obtaining a 
sample variance. 

A nondetcrministic function of time. 
Autocovariance function of y(t). 

See (58). 

Positive real number defined by 
r = T/t. 
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An intermediate term used in 
deriving (23). The definition of S 
is given by (64). 

One-sided (power) spectral density 
on a per hertz basis of the pure real 
function g(t). The dimensions of 
S,(j) are the dimensions of g 2 (t)/f. 
A definition for the measure of fre- 
quency stability. One-sided (power) 
spectral density of y(t) on a per 
hertz basis. The dimensions of 
S y (j) are Hz -1 . 

Time interval between the begin- 
nings of two successive measure- 
ments of average frequency. 

Time variable. 

An arbitrary fixed instant of time. 
The time coordinate of the begin- 
ning of the fcth measurement of 
average frequency. By definition, 
<* + , = <* + T, k = 0,1,2 •••. 
Dummy variable of integration; 
u = irfr. 

Instantaneous output voltage of 
signal generator. See (2). 

Nominal peak amplitude of signal 
generator output. See (2). 
Instantaneous voltage of reference 
signal. See (40). 

Peak amplitude of reference signal. 
See (40). 

Voltage output of ideal product 
detector. 

Low-pass filtered output of product 
detector. 

Real function of time related to the 
phase of the signal F(<) by x(t) s 
[<p( 0 ]/( 2 tv o ). 

A predicted value for x(t). 
Fractional frequency offset of F(f) 
from the nominal frequency. See (7). 
Average fractional frequency offset 
during the fcth measurement in- 
terval. See (9). 

The sample average of N successive 
values of y k . See (76). 
Nondeterministic (noise) function 
with (power) spectral density given 
by (25). 

Exponent of / for a power-law 
spectral density. 

Positive real constant. 

The Kronecker 5 function defined 
ji, if r = 1 
by Mr - 1) - -jo othemi<e 

Amplitude fluctuations of signal. 
See (2). 


a) 


*(0 

cl(N, T, r) 

K(N, T, r)> 

<r„(r) 

<rl(r) 

T 

m 

<p(t) 

+l(T, r) 
u = 2ir/ 

T HE measurement of frequency and fluctuations in 
frequency has received such great attention for 
so many years that it is surprising that the con- 
cept of frequency stability does not have a universally 
accepted definition. At least part of the reason has been 
that some uses are most readily described in the fre- 
quency domain and other uses in the time domain, as 
well as in combinations of the two. This situation is 
further complicated by the fact that only recently have 
noise models been presented that both adequately de- 
scribe performance and allow a translation between the 
time and frequency domains. Indeed, only recently has 
it been recognized that there can be a wide discrepancy 
between commonly used time domain measures them- 
selves. Following the NASA-IEEE Symposium on Short- 
Term Stability in 1964 and the Special Issue on Fre- 
quency Stability in the Peoceedings of the IEEE, 
February 1966, it now seems reasonable to propose a 
definition of frequency stability. The present paper is 
presented as technical background for an eventual IEEE 
standard definition. 

This paper attempts to present (as concisely as prac- 
tical) adequate, self-consistent definitions of frequency 
stability. Since more than one definition of frequency 
stability is presented, an important part of this paper 


Exponent of t. See (29). 
Instantaneous frequency of F(t). 
Defined by 


v(t) 
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Nominal (constant) frequency of 
V(t). 

The Fourier transform of n(t). 
Sample variance of N averages 
of y(t), each of duration t, and 
spaced every T units of time. 
See (10). 

Average value of the sample vari- 
ance <r*(AT, T, r). 

A second choice of the definition for 
the measure of frequency stability. 
Defined by = (<r’(V = 2, 

T = T, r)). 

Time stability measure defined by 
<rl(r) = r 2 al(r). 

Duration of averaging period of 
y(t) to obtain y t . See (9). 
Instantaneous phase of F(<). De- 
fined by 4>(/) = 2irv 0 t + <p(t). 
Instantaneous phase fluctuations 
about the ideal phase 2irv 0 t. See (2). 
Mean-square time error for Doppler 
radar. See (82). 

Angular Fourier frequency variable. 


I. Introduction 
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(perhaps the most important part) deals with transla- 
tions among the suggested definitions of frequency sta- 
bility. The applicability of these definitions to the more 
common noise models is demonstrated. 

Consistent with an attempt to be concise, the refer- 
ences cited have been selected on the basis of being of 
most value to the reader rather than on the basis of 
being exhaustive. An exhaustive reference list covering 
the subject of frequency stability would itself be a 
voluminous publication. 

Almost any signal generator is influenced to some ex- 
tent by its environment. Thus observed frequency in- 
stabilities may be traced, for example, to changes in 
ambient temperature, supply voltages, magnetic field, 
barometric pressure, humidity, physical vibration, or 
even output loading, to mention the more obvious. While 
these environmental influences may be extremely im- 
portant for many applications, the definition of fre- 
quency stability presented here is independent of these 
causal factors. In effect, we cannot hope to present an 
exhaustive list of environmental factors and a prescrip- 
tion for handling each even though, in some cases, these 
environmental factors may be by far the most im- 
portant. Given a particular signal generator in a partic- 
ular environment, one can obtain its frequency stability 
with the measures presented below, but one should not 
then expect an accurate prediction of frequency stability 
in a new environment. 

It is natural to expect any definition of stability to 
involve various statistical considerations such as sta- 
tionarity, ergodicity, average, variance, spectral density, 
etc. There often exist fundamental difficulties in rigorous 
attempts to bring these concepts into the laboratory. It 
is worth considering, specifically, the concept of sta- 
tionarity since it is a concept at the root of many statis- 
tical discussions. 

A random process is mathematically defined as sta- 
tionary if every translation of the time coordinate maps 
the ensemble onto itself. As a necessary condition, if one 
looks at the ensemble at one instant of time t, the dis- 
tribution in values within the ensemble is exactly the 
same as at any other instant of time t'. This is not to 
imply that the elements of the ensemble are constant 
in time, but, as one element changes value with time, 
other elements of the ensemble assume the previous val- 
ues. Looking at it in another way, by observing the 
ensemble at some instant of time, one can deduce no 
information as to when the particular instant was chosen. 
This same sort of invariance of the joint distribution 
holds for any set of times t lt 1 2 , • • •, t n and its transla- 
tion ti + r, fa + t, — , t n + t. 

It is apparent that any ensemble that has a finite 
past as well as a finite future cannot be stationary, and 
this neatly excludes the real world and anything of 
practical interest. The concept of stationarity does vio- 
lence to concepts of causality since we implicitly feel 
that current performance (i.e., the applicability of sta- 
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tionary statistics) cannot be logically dependent upon 
future events (i.e., if the process is terminated some time 
in the distant future). Also, the verification of station- 
arity would involve hypothetical measurements that are 
not experimentally feasible, and therefore the concept of 
stationarity is not directly relevant to experimentation. 

Actually the utility of statistics is in the formation 
of idealized models that reasonably describe significant 
observables of real systems. One may, for example, con- 
sider a hypothetical ensemble of noises with certain 
properties (such as stationarity) as a model for a par- 
ticular real device. If a model is to be acceptable, it 
should have at least two properties: first, the model 
should be tractable ; that is, one should be able to easily 
arrive at estimates for the elements of the models; and 
second, the model should be consistent with observables 
derived from the real device that it is simulating. 

Notice that one does not need to know that the device 
was selected from a stationary ensemble, but only that 
the observables derived from the device are consistent 
with, say, elements of a hypothetically stationary en- 
semble. Notice also that the actual model used may 
depend upon how clever the experimenter-theorist is in 
generating models. 

It is worth noting, however, that while some texts 
on statistics give “tests for stationarity,” these tests are 
almost always inadequate. Typically, these tests de- 
termine only if there is a substantial fraction of the 
noise power in Fourier frequencies whose periods are of 
the same order as the data length or longer. While this 
may be very important, it is not logically essential to 
the concept of stationarity. If a nonstationary model 
actually becomes common, it will almost surely be be- 
cause it is useful or convenient and not because the 
process is “actually nonstationary.” Indeed, the phrase 
“actually nonstationary” appears to have no meaning 
in an operational sense. In short, stationarity (or non- 
stationarity) is a property of models, not a property of 
data [1]. 

Fortunately, many statistical models exist that ade- 
quately describe most present-day signal generators; 
many of these models are considered below. It is obvious 
that one cannot guarantee that all signal generators are 
adequately described by these models, but the authors 
do feel they are adequate for the description of most 
signal generators presently encountered. 

II. Statement of the Problem 

To be useful, a measure of frequency stability must 
allow one to predict performance of signal generators 
used in a wide variety of situations as well as allow 
one to make meaningful relative comparisons among 
signal generators. One must be able to predict perform- 
ance in devices that may most easily be described either 
in the time domain, or in the frequency domain, or in 
a combination of the two. This prediction of perform- 
ance may involve actual distribution functions, and thus 
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second moment measures (such as power spectra and 
variances) are not totally adequate. 

Two common types of equipment used to evaluate the 
performance of a frequency source are (analog) spectrum 
analyzers (frequency domain) and digital electronic 
counters (time domain). On occasion the digital counter 
data are converted to power spectra by computers. One 
must realize that any piece of equipment simultaneously 
has certain aspects most easily described in the time 
domain and other aspects most easily described in the 
frequency domain. For example, an electronic counter 
has a high-frequency limitation, an experimental spectra 
are determined with finite time averages. 

Research has established that ordinary oscillators dem- 
onstrate noise, which appears to be a superposition of 
causally generated signals and random nondeterministic 
noises. The random noises include thermal noise, shot 
noise, noises of undetermined origin (such as flicker 
noise) , and integrals of these noises. 

One might well expect that for the more general cases 
one would need to use a nonstationary model (not sta- 
tionary even in the wide sense, i.e., the covariance sense) . 
Nonstationarity would, however, introduce significant dif- 
ficulties in the passage between the frequency and time 
domains. It is interesting to note that, so far, experi- 
menters have seldom found a nonstationary (covariance) 
model useful in describing actual oscillators. 

In what follows, an attempt has been made to separate 
general statements that hold for any noise or perturba- 
tion from the statements that apply only to specific mod- 
els. It is important that these distinctions be kept in 
mind. 


sent the one-sided spectral density of the (pure real) 
function g(t) on a per hertz basis; that is, the total 
“power” or mean-square value of g(t) is given by 

f S.(j) dj. 

Jo 

Since the spectral density is such an important con- 
cept to what follows, it is worthwhile to present some 
important references on spectrum estimation. There are 
many references on the estimation of spectra from data 
records, but worthy of special note are [2] — [5 ] . 

IV. Definition of Measures of Frequency Stability 
(Second-Moment Type) 


A. General 

Consider a signal generator whose instantaneous out- 
put voltage V [t) may be written as 

V(0 = [Vo + e(0] sin [2 W + *>(0] (2) 

where T 0 and v u are the nominal amplitude and fre- 
quency, respectively, of the output and it is assumed 
that 

kol ~ 


and 


|2mJ " * 

for substantially all time t. Making use of (1) and (2) 
one sees that 


To 


<p( 0 
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III. Background and Definitions 

To discuss the concept of frequency stability imme- 
diately implies that frequency can change with time and 
thus one is not considering Fourier frequencies (at least 
at this point). The conventional definition of instantan- 
eous (angular) frequency is the time rate of change of 
phase; that is 


2 mo - - m 


( 1 ) 


where $(() is the instantaneous phase of the oscillator. 
This paper uses the convention that time-dependent 
frequencies of oscillators are denoted by r(£) (cycle fre- 
quency, hertz), and Fourier frequencies are denoted by 
u> (angular frequency) or / (cycle frequency, hertz) where 
w = % rf. In order for (1) to have meaning, the phase $(t) 
must be a well-defined function. This restriction imme- 
diately eliminates some “nonsinusoidal” signals such as 
a pure random uncorrelated (“white”) noise. For most 
real signal generators, the concept of phase is reasonably 
amenable to an operational definition and this restric- 
tion is not serious. 

Of great importance to this paper is the concept of 
spectral density, S g (J). The notation S„(f) is to repre- 


4>(<) = 2irv a t + ip(t) 

(5) 

II 
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Equations (3) and (4) are essential in order that <p{t) 
may be defined conveniently and unambiguously (sec 
measurement section). 

Since (4) must be valid even to speak of an instantan- 
eous frequency, there is no real need to distinguish 
stability measures from instability measures. That is, 
any fractional frequency stability measure will be far 
from unity, and the chance of confusion is slight. It is 
true that in a very strict sense people usually measure 
instability and speak of stability. Because the chances of 
confusion are so slight, the authors have chosen to con- 
tinue in the custom of measuring “instability” and speak- 
ing of stability (a number always much less than unity). 

Of significant interest to many people is the radio fre- 
quency (RF) spectral density Sy{f). This is of direct 
concern in spectroscopy and radar. However, this is not 
a good primary measure of frequency stability for two 
reasons. First, fluctuations in the amplitude «(f) contrib- 
ute directly to S r (/); and second, for many cases when 
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«(£) is insignificant, the RF spectrum Sv(f) is not 
uniquely related to the frequency fluctuations [6]. 

B. General: First Definition of the Measure of Frequency 
Stability — Frequency Domain 

By definition, let 



(7) 


where <p(t) and f 0 are as in (2). Thus y{t) is the in- 
stantaneous fractional frequency deviation from the nom- 
inal frequency v 0 - A proposed definition of frequency 
stability is the spectral density S„ (/) of the instantaneous 
fractional frequency fluctuations y{t). The function S„(/) 
has the dimensions of Hz -1 . 

One can show [7] that if <S P (/) is the spectral density 
of the phase fluctuations, then 

SM ' 

Thus a knowledge of the spectral density of the phase 
fluctuations S v (f) allows a knowledge of the spectral 
density of the frequency fluctuations jS„ (/) , the first def- 
inition of frequency stability. Of course, S„(/) cannot 
be perfectly measured— this is the case for any physical 
quantity; useful estimates of S„(/) are, however, easily 
obtainable. 


C. General: Second Definition of the Measure of Fre- 
quency Stability — Time Domain 

The second definition is based on the sample variance of 
the fractional frequency fluctuations. In order to present 
this measure of frequency stability, define y k by the 
relation t 


1 r“ +T 

V* = “ / y(t) dt 


<p(tj + t) — <p(t k ) 
2 TTVqT 


(9) 


where <* +I = <* + T, k = 0, 1, 2, • ■ • , T is the repetition 
interval for measurements of duration r, and t 0 is arbitrary. 
Conventional frequency counters measure the number of 
cycles in a period t; that is, they measure f 0 t(1 + y k ). 
When r is 1 8 they count the number of f 0 (1 + Vt)- 
The second measure of frequency stability, then, is 
defined in analogy to the sample variance by the relation 

(cl(N, T, r)> - £ (*. y*J) , (10) 


where ( g ) denotes the infinite time average of g. This 
measure of frequency stability is dimensionless. 

In many situations it would be wrong to assume that 
(10) converges to a meaningful limit as N — ► ». First, 
of course, one cannot practically let N approach infinity 
and, second, it is known that some actual noise processes 
contain substantial fractions of the total noise power in 


the Fourier frequency range below one cycle per year. 
In order to improve comparability of data, it is important 
to specify particular N and T. For the preferred definition 
we recommend choosing N = 2 and T = r (i.e., no dead 
time between measurements). Writing (o\(N = 2, T = r, r)) 
as cr*(r), the Allan variance [8], the proposed measure of 
frequency stability in the time domain may be written as 

°»( T ) = (ID 


for T = t. 

Of course, the experimental estimate of <r\{r) must be 
obtained from finite samples of data, and one can never 
obtain perfect confidence in the estimate; the true time 
average is not realizable in a real situation. One estimates 
<r*(r) from a finite number (say, m) of values of <rj( 2, r, r) 
and averages to obtain an estimate of <rj(r). Appendix I 
shows that the ensemble average of <r*(2, r, r) is convergent 
(i.e., as m — * <») even for noise processes that do not have 
convergent (<r*(JV, r, r)) as N — »• <». Therefore, <r“(r) has 
greater utility as an idealization than does (<rj(<», T , r)) 
even though both involve assumptions of infinite averages. 
In effect, increasing N causes &l(N, T, r) to become more 
sensitive to the low-frequency components of S u (f). In 
practice, one must distinguish between an experimental 
estimate of a quantity (say, of v*(t)) and its idealized 
value. It is reasonable to believe that extensions to the 
concept of statistical ("quality”) control [9] may prove 
useful here. One should, of course, specify the actual 
number m of independent samples used for an estimate 
of ,»( r). 

In summary, therefore, S,(J) is the proposed measure of 
(instantaneous) frequency stability in the (Fourier) 
frequency domain and o % y (r) is the proposed measure of 
frequency stability in the time domain. 

D. Distributions 

It is natural that people first become involved with 
second moment measures of statistical quantities and only 
later with actual distributions. This is certainly true with 
frequency stability. While one can specify the argument 
of a distribution function to be, say (y k+l — y k ), it makes 
sense to postpone such a specfication until a real use has 
materialized for a particular distribution function. This 
paper does not attempt to specify a preferred distribution 
function for frequency fluctuations. 

E. Treatment of Systematic Variations 

1 ) General: The definition of frequency stability <r’(r) 
in the time domain is useful for many situations. However, 
some oscillators, for example, exhibit an aging or almost 
linear drift of frequency with time. For some applications, 
this trend may be calculated and should be removed [8] 
before estimating <rl(r). 

In general, a systematic trend is perfectly deterministic 
(i.e., predictable) while the noise is nondeterministic. 
Consider a function g(t), which may be written in the form 
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g(t) = c(t) + n(t) (12) 

where c(t) is some deterministic function of time and n(t), 
the noise, is a nondeterministic function of time. We will 
define c(t) to be the systematic trend to the function g(t). 
A problem of significance here is to determine when and 
in what sense c(t) is measurable. 

2) Specific Case — Linear Drift: As an example, if we 
consider a typical quartz crystal oscillator whose fractional 
frequency deviation is y(t), we may let 

g(t) = f t y(t). (13) 

With these conditions, c(t) is the drift rate of the oscil- 
lator (e.g., 10~ 10 /day) and n{t) is related to the fre- 
quency “noise” of the oscillator by a time derivative. 
One sees that the time average of g{t) becomes 


1 /*<•+ T i <*<• + T 

Tfi J g(t) dt = c x - f - J n(t ) 


dt 


(14) 


where c(t) = c x is assumed to be the constant drift rate 
of the oscillator. In order for to be an observable, 
it is natural to expect the average of the noise term to 
vanish, that is, converge to zero. 

It is instructive to assume [8], [10] that in addition 
to a linear drift, the oscillator is perturbed by a flicker 
noise, i.e., 


\h- if\ 0 < / < ft 
10 , . / > ft 

where h _! is a constant (see Section V-A-2) and thus, 

m - a® 

lo, i > h 

for the oscillator we are considering. With these assump- 


S,(J) 


(15) 


tions, it is seen that 

1 r‘ ,+T 

lim ■= I n(t ) dt = k(0) = 0 

r—oo 1 J t % 


(17) 


and that 


lim s variance 

r-»«o l 


[!£‘ +r n(t)df]} = 0 (18) 


where «(/) is the fourier transform of n(t). Since S„(0) 
= 0, k (0) must also vanish both in probability and in 
mean square. Thus, not only does n(t) average to zero, 
but one may obtain arbitrarily good confidence on the 
result by longer averages. 

Having shown that one can reliably estimate the drift 
rate Ci of this (common) oscillator, it is instructive to 
attempt to fit a straight line to the frequency aging. 
That is, let 


g(t) = y( 0 


(19) 


and thus 


( 20 ) 


where c*, is the frequency intercept at t = t 0 and Ci is 
the drift rate previously determined. A problem arises 
here because 


S..Q) = S,(j) (21) 


! Ti 




< variance - 

^ / n'(0 dt 

ilk 

[= OO 

(22) 


for the noise model we have assumed. This follows from 
the fact that the (infinite N) variance of a flicker noise 
process is infinite [7], [8], [10]. Thus, c<> cannot be 
measured with any realistic precision, at least, in an 
absolute sense. 

We may interpret these results as follows. After ex- 
perimenting with the oscillator for a period of time one 
can fit an empirical equation to y{t) of the form 

2/(0 = c 0 + tc, + n'(0, 

where n’ (t) is nondeterministic. At some later time it is 
possible to reevaluate the coefficients c 0 and C\. Accord- 
ing to wffiat has- been said, the drift rate should be 
reproducible to wdthin the confidence estimates of the 
experiment regardless of when it is reevaluated. For c„, 
however, this is not true. In fact, the more one attempts 
to evaluate Co, the larger the fluctuations are in the 
result. 

Depending on the spectral density of the noise term, 
it may be possible to predict future measurements of 
c 0 and to place realistic confidence limits on the predic- 
tion [11]. For the case considered here, however, these 
confidence limits tend to infinity when the prediction 
interval is increased. Thus, in a certain sense, c 0 is 
“measurable” but it is not in statistical control (to use 
the language of the quality control engineer [9]). 


V. Translations Among Frequency Stability 
Measures 

A. Frequency Domain to Time Domain 

1) General: It is of value to define r = T/t\ that is, 
r is the ratio of the time interval between successive 
measurements to the duration of the averaging period. 
Cutler has shown (see Appendix I) that 


T, r)> 


(N - 1) /o 


[sin 2 (w /r)] 
(tt/t ) 2 


sin 2 (irrf N t) \ 
N 2 sin 2 (nr/ r)J 

(23) 


Equation (23) in principle allows one to calculate the 
time-domain stability (<r 2 v (N, T, r)} from the frequency- 
domain stability &„(/). 

2) Specific Model: A model that has been found use- 
ful [8], [ 10]— [ 13] consists of a set of five independent 
noise processes z„(<), n = —2, —1, 0, 1, 2, such that 


g(t) = c 0 + C& - t 0 ) + n’{t) 
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2/(0 = z-i(t) + + z o (0 + Zi(0 + z 2 (0 (24) 

and the spectral density of z„ is given by 

5..(f) - K ° ^ (25) 

10, / > U,n = -2, -1,0, 1,2, 

where the /i„ are constants. Thus, S„(/) becomes 

<S„(/) = A- j/ 2 + A_,/ 1 + /io + hif + h 2 f, (26) 

for 0 < / < /* and <S,(/) is assumed to be negligible beyond 
this range. In effect, each z„ contributes to both S u (j) and 
( a 2 ,(N , T, t)) independently of the other z„. The con- 
tributions of the z„ to ( o\(N , T, r)) are tabulated in 
Appendix II. 

Any electronic device has a finite bandwidth and this 
certainly applies to frequency-measuring equipment also. 
For fractional frequency fluctuations y(t ) whose spectral 
density varies as 

S,(J)~r, «>- 1 (27) 

for the higher Fourier components, one sees (from 
Appendix I) that ( <rl(N , T, r)) may depend on the exact 
shape of the frequency cutoff. This is true because a 
substantial fraction of the noise “power” may be in these 
higher Fourier components. As a simplifying assumption, 
this paper assumes a sharp cutoff in noise “power” at the 
frequency /» for the noise models. It is apparent from the 
tables of Appendix II that the time domain measure of 
frequency stability may depend on /* in a very important 
way, and, in some practical cases, the actual shape of the 
frequency cutoff may be very important [7]. On the 
other hand, there are many practical measurements 
where the value of /* has little or no effect. Good practice, 
however, dictates that the system noise bandwidth /* 
should be specified with any results. 

In actual practice, the model of (24)-(26) seems to fit 
almost all real frequency sources. Typically, only two or 
three of the fi-coefficients are actually significant for a 
real device and the others can be neglected. Because of 
its applicability, this model is uoed in much of what 
follows. Since the z„ are assumed to be independent noises, 
it is normally sufficient to compute the effects for a 
general z„ and recognize that the superposition can be 
accomplished by simple additions for their contributions 
to S V (J) or (al(N, T, r)>. 

B. Time Domain to Frequency Domain 

1) General: For general ( <rl(N , T, t)) no simple pre- 
scription is available for translation into the frequency 
domain. For this reason, one might prefer £„(/) as a 
general measure of frequency stability. This is especially 
true for theoretical work. 

2) Specific Model: Equations (24)-(26) form a realistic 
model that fits the random nondeterministic noises found 
on most signal generators. Obviously, if this is a good 
model, then the tables in Appendix II may be used 
(in reverse) to translate into the frequency domain. 


Allan [8] and Vessot [12] showed that if 
S„(f) = J ha ^’ 0 ~ 1 - 


2tt/ a » 1 


where a is a constant, then 
T, t)> - 


for N and r = T/r held constant. The constant a i& 
related to a by the mapping shown 1 in Fig. 1. If (28) 
and (29) hold over a reasonable range for a signal gen- 
erator, then (28) can be substituted into (23) and evaluated 
to determine the constant h a from measurements of 
(<r*(fV, T, r)). It should be noted that the model of (28) 
and (29) may be easily extended to a superposition of 
similar noises as in (26). 

C. Translations Among the Time-Domain Measures 

1) General: Since {a\{N, T, r)) is a function of N, T, 
and t (for some types of noise /» is also important), it in 
very desirable to be able to translate among different 
sets of N, T, and r (j k held constant). This is, however, 
not possible in general. 

2) Specific Model: It is useful to restrict consideration 
to a case described by (28) mid (29). Superpositions of 
independent noises with different power-law types of 
spectral densities (i.e., different a) can also be treated by 
this technique, e.g., (26). One may define two “bias 
functions,” B , and B 2 by the relations [13] 

B t (N, r, M ) = ( ,» (2t T t)) (30)- 


D /_ \ _ (<r 2 (2, T , t)) . . 

Bi{r ’ m) = <^(2, r, r)> (31) 

where r s= T/r and /i is related to a by the mapping of 
Fig. 1. In words, is the ratio of the average variance- 
for N samples to the average variance for two samples 
(everything else held constant), while B 2 is the ratio of 
the average variance with dead time between measure- 
ments (r p* 1) to that of no dead time (r = 1 and with 
N = 2 and r held constant). These functions are tabulated 
in [13]. Figs. 2 and 3 show a computer plot of 
Bi(N, r = 1, u) and B a (r, #«)• 

Suppose one has an experimental estimate of (o-J(A' I , 
T» rO) and its spectral type is known, i.e., (28) and (29) 
form a good model and p is known. Suppose also that one 
wishes to know the variance at some other set of measure- 
ment parameters N 2 , T 2 , r 2 ■ An unbiased estimate of 
(crl(N 2 , T 2 , t 3 )) may be calculated by 


1 It should be noted that in Allan [8], the exponent a cor- 
responds to the spectrum of phase fluctuations while variances 
are taken over average Irequency fluctuations. In the present 
paper, a is identical to the exponent a - f 2 in [8], 
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Fig. 2. Function B t (N, r = 1, /*). 


T„ r 2 )> 



1*2, ii)B 2 ( r 2 . m) 
.Bi(N i , f, , n)B 2 (r l , n)_ 


(<rl(N t , 


T u 


T\)), 


(32) 


-where ri = Tj/ti and r 2 = T 2 /t 2 . 

3) General: While it is true that the concept of the 
"bias functions J3i and B 2 could be extended to other 
processes besides those with the power-law types of 
spectral densities, this generalization has not been done. 
Indeed, spectra of the form given in (28) [or super- 
positions of such spectra as in (26) ] seem to be the 
most common types of nondeterministic noises encoun- 
tered in signal generators and associated equipment. For 
other types of fluctuations (such as causally generated 
perturbations), translations must be handled on an in- 
-dividual basis. 


VI. Applications of Stability Measures 

Obviously, if one of the stability measures is exactly 
dhe important parameter in the use of a signal generator, 
-the stability measure’s application is trivial. Some non- 
'irivial applications arise when one is interested in a dif- 


Fig. 3. Bias function B 2 (r, /i). 

ferent parameter, such as in the use of an oscillator in 
Doppler radar measurements or in clocks. 

A. Doppler Radar 

1) General: From its transmitted signal, a Doppler 
radar receives from a moving target a frequency-shifted 
return signal in the presence of other large signals. These 
large signals can include clutter (ground return) and 
transmitter leakage into the receiver (spillover). In- 
stabilities of radar signals result in noise energy on the 
clutter return, on spillover, and on local oscillators in 
the equipment. 

The limitations of subclutter visibility (SCV) rejec- 
tions due to the radar signals themselves are related to 
the RF power spectral density S v (/) . The quantity typi- 
cally referred to is the carrier-to-noise ratio and can be 
mathematically approximated by the quantity 

Sy(f) 

r svV'w 

■ j ° 

The effects of coherence of target return and other 
radar parameters are amply considered in the literature 
[14] — [17] . 

2) Special Case: Because FM effects generally pre- 
dominate over AM effects, this carrier-to-noise ratio is 
approximately given by [6] 

» *S,(|/ - p 0 |). (33) 

S v (f)dr 

Jo 

for many signal sources provided \j — v 0 \ is sufficiently 
greater than zero. (The factor of \ arises from the fact 
that S p (f) is a one-sided spectrum.) Thus, if / — to is 
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a frequency separation from the carrier, the carrier-to- 
noise ratio at that point is approximately 

- vq I) = | (^—Js.(\l - f„|). (34) 

B. Clock Errors 

1) General: A clock is a device that counts the cycles 
of a periodic phenomenon. Thus, the reading error. x(t) 
of a clock run from the signal given by (2) is 



(35) 


and the dimensions of x{t) are seconds. 

If this clock is a secondary standard, then one could 
have available some past history of x(t), the time error 
relative to the standard clock. It often occurs that one 
is interested in predicting the clock error x(t) for some 
future date, say t 0 + t, where t 0 is the present date. 
Obviously, this is a problem in pure prediction and can 
be handled by conventional methods [3]. 

2) Special Case: Although one could handle the predic- 
tion of clock errors by the rigorous methods of prediction 
theory, it is more common to use simpler prediction 
methods [10], [11]. In particular, one often predicts a clock 
error for the future by adding to the present error a 
correction that is derived from the current rate of gain 
(or loss) of time. That is, the predicted error £(t 0 + r) 
is related to the past history of x(t) by 

x(/o + r) = x(/ 0 ) 4- r[ x(tn) ~ * (to ~ r) ] • (36) 


It is typical to let T = r. 

Thus, the mean-square error of prediction for T = r 
becomes 


([x(<0 + r) — £(l 0 + r)] 2 > 

= ([x(t 0 + t) - 2x(l 0 ) + X(t 0 ~ r)] 2 ), (37) 

which, with the aid of (11), can be written in the form 
<[*00 + r) - *(fo + r)] 2 ) = 2rV u (r). (38) 


y(t). One of the most common techniques is a heterodyne 
or beat frequency technique. In this method, the signal 
from the oscillator to be tested is mixed with a reference 
signal of almost the same frequency as the test oscillator 
in order that one is left with a lower average frequency 
for analysis without reducing the frequency (or phase) 
fluctuations themselves. Following Vessot et al. [18],. 
consider an ideal reference oscillator whose output signal 
is 

V,{t) = For sin 2irv 0 t (40) 

and a second oscillator whose output voltage F(<) is 
given by (2): F(<) = [F 0 + *(<)] sin [2vv 0 t + y>(0]- Let 
these two signals be mixed in a product detector; that is, 
the output of the product detector v(t) is equal to the 
product yV(t) X F,(f), where y is a constant (see Fig. 4). 

Let v(t), in turn, be processed by a sharp low-pass filter 
with cutoff frequency /< such that 

0 </»</!<*,. (41)' 

One may write 
yV(l)-V r (t) 

= y F o r ( F 0 + e)[sin 2irv 0 f][sin (2irv 0 t + <p)] 

= v(t) = y ( ^ ^ ( 1 + -L-)[cos tp — COS (470-0^ 4 <p)]. 

(42) 

Assume that cos [v>(0] has essentially no power in Fourier 
frequencies / in the region / > /(. The effect of the low-pass 
filter then is to remove the second term on the extreme 
right of (42) ; that is 

v\t) = 7 (l 4- -pr) COS <p(t) . (43)' 

This separation of terms by the filter is correct only if 
|[^(0/2 tti/ o ]| « 1 for all t (4). 

The following two cases are of interest. 

Case I: The relative phase of the oscillators is ad- 
justed so that \f{t)\ <K 1 (in-phase condition) during: 
the period of measurement. Under these conditions 


One can define a time stability measure u 2 (r) by 

<r 2 (r) = tV 2 (t). (39) 

Clearly, however, the actual errors of prediction of clock 
readings are dependent on the prediction algorithm used 
and the utility of such a definition as u 2 (t) is not great. 
Caution should be used in employing this definition. 

VII. Measurement Techniques for Frequency 
Stability 

A. Heterodyne Techniques (General) 

It is possible for oscillators to be very stable and 
values of <r v {r) can be as small as 10 -14 in some state-of- 
the-art equipment. Thus, one often needs measuring tech- 
niques capable of resolving very small fluctuations in 


v\l)^\V 0T V, + |F„ r «(0, (44> 

since cos <?(t ) ~ 1. That is to say one detects the ampli- 
tude noise e(t) of the signal. 

Case II: The relative phase of the oscillators is ad- 
justed to be in approximate quadrature; that is 

<p\t)=<p(t) 4- 1 (45)' 

where |/(t)| <SC 1. Under these conditions, 

cos <p(t) = sin ip'(t) Pa y>'(t) (46) 

and 

v'(t) = l V ar V 0 <p'(t) + l V 0r <p'(LMt). (47) 
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Fig. 4. Heterodyne scheme. 


If it is true that |[«(£)/F 0 ]| «. 1 for all t (3), then 
(47) becomes 

v'(t) V 0r Vov'(t)-, (48) 

that is, t/(<) is proportional to the phase fluctuations. 
Thus, in order to observe /(£) by this method, (3) and 
(4) must be valid. For different average phase values, 
mixtures of amplitude and phase noise are observed. 

In order to maintain the two signals in quadrature for 
long observational periods, the reference oscillator can 
be a voltage-controlled oscillator (VCO) and one may 
feed back the phase error voltage as defined in (48) to 
control the frequency of the VCO [19]. In this condition 
of the phase- locked oscillator, the voltage v'(t) is the 
analog of the phase fluctuations for Fourier frequencies 
above the loop cutoff frequency of the locked loop. For 
Fourier frequencies below the loop cutoff frequency of 
the loop, v'(t) is the analog of frequency fluctuations. 
In practice, one should measure the complete servo-loop 
response. 

B. Period Measurement 

Assume one has an oscillator whose voltage output 
may be represented by (2). If |[<(t)/V 0 ]| 1 for all 

£ and the total phase 

4>(£) = 2 7rv 0 t + <p{t) (5) 

is a monotonic function of time (that is, |[<p(£)/2jtv 0 ]| < 1), 
then the time t between successive positive going zero 
crossings of V(t) is related to the average frequency during 
the interval r. Specifically 

- = r 0 (l + y.). (49) 

T 

If one lets r be the time between a positive going zero 
crossing of 7(£) and the Mth successive positive going 
zero crossing, then 

— = r 0 (l + V.). (50) 

T 

If the variations At of the period are small compared to 
the average period t 0 , Cutler and Searle [7] have shown 


that one may make a reasonable approximation to 
(<rl(N, T, t 0 )) using period measurements. 

C. Period Measurement With Heterodyning 

Suppose that ?(t) is a monotonic function of time. 
The output of the filter of Section VII-A (43) becomes 

v'(t) « y cos <p(t) (51) 

if | [e (f) /Vo] | <K 1. Then one may measure the period 
t of two successive positive zero crossings of v'(t). Thus 

~ = "0 [27*1 (52) 

and for the Mth positive crossover 
M 

f = Vo 10 . 1 . (53) 

The magnitude bars appear because cos <p(t) is an even 
function of <?(t). It is impossible to determine by this 
method alone whether <p is increasing with time or de- 
creasing with time. Since y„ may be very small 
or 10 -12 for very good oscillators) , t may be quite long 
and thus measurable with a good relative precision. 

If the phase <p(t) is not monotonic, the true y„ may be 
near zero but one could still have many zeros of cos <p(t) 
and thus (52) and (53) would not be valid. 

D. Frequency Counters 

Assume the phase (either or <p) is a montonic func- 
tion of time. If one counts the number M of positive going 
zero crossings in a period of time r, then the average fre- 
quency of the signal is M/r. If we assume that the signal 
is F(£) as defined in (2), then 

^ = r„(l + y n ). (54) 

If we assume that the signal is t /(£) as defined in (48), 
then 

— = v 0 |y„|. (55) 

T 

Again, one measures only positive frequencies. 

E. Frequency Discriminators 

A frequency discriminator is a device that converts 
frequency fluctuations into an analog voltage by means 
of a dispersive element. For example, by slightly detuning 
a resonant circuit from the signal V(t) the frequency 
fluctuations (1/2 ir)]^(£) are converted to amplitude fluc- 
tuations of the output signal. Provided the input amplitude 
fluctuations [t{t))/V 0 are insignificant, the output ampli- 
tude fluctuations can be a good measure of the frequency 
fluctuations. Obviously, more sophisticated frequency 
discriminators exist (e.g., the cesium beam). 
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From the analog voltage one may use analog spectrum 
analyzers to determine S v (f) , the frequency stability. By 
converting to digital data, other analyses are possible 
on a computer. 

F. Common Hazards 

1) Errors Caused by Signal-Processing Equipment: The 
intent of most frequency stability measurements is to 
evaluate the source and not the measuring equipment. 
Thus, one must know the performance of the measuring 
system. Of obvious importance are such aspects of the 
measuring equipment as noise level, dynamic range, 
resolution (dead time), and frequency range. 

It has been pointed out that the noise bandwidth / A is 
very essential for the mathematical convergence of certain 
expressions. Insofar as one wants to measure the signal 
source, one must know that the measuring system is not 
limiting the frequency response. At the very least, one 
must recognize that the frequency limit of the measuring 
system may be a very important, implicit parameter for 
either <r 2 (<) or S y (j). Indeed, one must account for any 
deviations of the measuring system form ideality such as 
a “nonflat” frequency response of the spectrum analyzer 
itself. 

Almost any electronic circuit that processes a signal 
will, to some extent, convert amplitude fluctuations at the 
input terminals into phase fluctuations at the output. 
Thus, AM noise at the input will cause a time-varying 
phase (or FM noise) at the output. This can impose im- 
portant constraints on limiters and automatic gain control 
(AGC) circuits when good frequency stability is needed. 
Similarly, this imposes constraints on equipment used for 
frequency stability measurements. 

2) Analog Spectrum Analyzers (Frequency Domain) : 
Typical analog spectrum analyzers are very similar in 
design to radio receivers of the superheterodyne type, and 
thus certain design features are quite similar. For exam- 
ple, image rejection (related to predetection bandwidth) 
is very important. Similarly, the actual shape of the 
analyzer’s frequency window is important since this af- 
fects spectral resolution. As with receivers, dynamic 
range can be critical for the analysis of weak signals in 
the presence of substantial power in relatively narrow 
bandwidths (e.g., 60 Hz). 

The slewing rate of the analyzer must be consistent 
with the analyzer’s frequency window and the post-detec- 
tion bandwidth. If one has a frequency window of 1 Hz, 
one cannot reliably estimate the intensity of a bright 
line unless the slewing rate is much slower than 1 Hz/s. 
Additional post-detection filtering will further reduce the 
maximum usable slewing rate. 

3) Spectral Density Estimation from Time Domain 
Data: It is beyond the scope of this paper to present a 
comprehensive list of hazards for spectral density estima- 
tion; one should consult the literature [2]-[5]. There 


are a few points, however, which are worthy of special 
notice: a) data aliasing (similar to predetection band- 
width problems); b) spectral resolution; and c) con- 
fidence of the estimate. 

4) Variances of Frequency Fluctuations o-’(t): It is not 
uncommon to have discrete frequency modulation of a 
source such as that associated with the power supply 
frequencies. The existence of discrete frequencies in S y (J) 
can cause <t 2 (t) to be a very rapidly changing function 
of r. An interesting situation results when r is an exact 
multiple of the period of the modulation frequency (e.g., 
one makes r = 1 s and there exists 60-Hz frequency 
modulation on the signal). In this situation, a\(r = 1 s) 
can be very optimistic relative to values with slightly 
different values of r. 

One also must be concerned with the convergence 
properties of <rl(r) since not all noise processes will have 
finite limits to the estimates of <r*(r) (see Appendix I). 
One must be as critically aware of any “dead time” in the 
measurement process as of the system bandwidth. 

5) Signal Source and Loading: In measuring frequency 
stability one should specify the exact location in the 
circuit from which the signal is obtained and the nature 
of the load used. It is obvious that the transfer character- 
istics of the device being specified will depend on the load 
and that the measured frequency stability might be 
affected. If the load itself is not constant during the 
measurements, one expects large effects on frequency 
stability. 

6) Confidence of the Estimate: As with any measurement 
in science, one wants to know the confidence to assign to 
numerical results. Thus, when one measures S y (j) or <r*(r), 
it is important to know the accuracies of these estimates. 

o) The Allan Variance: It is apparent that a single 
sample variance <r’(4, r, r) does not have good confidence, 
but, by averaging many independent samples, one can 
improve the accuracy of the estimate greatly. There is a 
key point in this statement, “independent samples.” For 
this argument to be true, it is important that one sample 
variance be independent of the next. Since </J(2, r, r) is 
related to the first difference of the frequency (11), 
it is sufficient that the noise perturbing y(t) have “inde- 
pendent increments,” i.e., that y(t) be a random walk. 
In other words, it is sufficient that S,(j) ~ /“ 2 for low 
frequencies. One can show that for noise processes that 
are more divergent at low frequencies than /” 2 , it is 
difficult (or impossible) to gain good confidence on 
estimates of <rj(r). For noise processes that are less 
divergent than /~ 2 , no problem exists. 

It is worth noting that if we were interested in 
<rl(N = oo, r, t), then the limit noise would become 
S y (j) ~ f instead of f~ 2 as it is for <r 2 (2, r, t). Since most 
real signal generators possess low-frequency divergent 
noises, (<r\(2, r, r)) is more useful than <?\(N = <*, r, r). 

Although the sample variances <z 2 (2, r, r) will not be 
normally distributed, the variance of the average of m 
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independent (nonoverlapping) samples of <rf(2, r, r) 
(i.e., the variance of the Allan variance) will decrease as 
1/m provided the conditions on low-frequency divergence 
are met. For sufficiently large m, the distribution of the 
m sample averages of ov(2, r, t) will tend toward normal 
(central limit theorem). It is thus possible to estimate 
confidence intervals based on the normal distribution. 

As always, one may be interested in r values approach- 
ing the limits of available data. Clearly, when one is 
interested in r values of the order of a year, one is severely 
limited in the size of m, the number of samples of trj(2, r, r). 
Unfortunately, there seems to be no substitute for many 
samples and one extends r at the expense of confidence in 
the results. “Truth in packaging” dictates that the sample 
size m be stated with the results. 

6) Spectral Density: As before, one is referred to the 
literature for discussions of spectrum estimation [2]-[5]. 
It is worth pointing out, however, that for S v (j) there are 
basically two different types of averaging that can be 
employed: sample averaging of independent estimates 
of <S„(/), and frequency averaging where the resolution 
bandwidth is made much greater than the reciprocal data 
length. 


VIII. Conclusions 

A good measure of frequency stability is the spectral 
density S„(j) of fractional frequency fluctuations y(t). 
An alternative is the expected variance of N sample 
averages of y(t) taken over a duration r. With the begin- 
ning of successive sample periods spaced every T units 
of time, the variance is denoted by a l(N, T, r). The 
stability measure, then, is the expected value of many 
measurements of <r 2 ,(N, T, r) with N = 2 and T = r; 
that is, (r). For all real experiments one has a finite 
bandwidth. In general, the time domain measure of 
frequency stability v*(r) is dependent on the noise band- 
width of the system. Thus, there are four important 
parameters to the time domain measure of frequency 
stability. 

N Number of sample averages (N = 2 for preferred 
measure). 

T Repetition time for successive sample averages 
(T = r for preferred measure). 

t Duration of each sample average. 

U System noise bandwidth. 

Translations among the various stability measures for 
common noise types are possible, but there are significant 
reasons for choosing N — 2 and T = r for the preferred 
measure of frequency stability in the time domain. This 
measure, the Allan variance, ( N = 2) has been referenced 
by [12], [20]— [22] and more. 

Although S V (J) appears to be a function of the single 
variable /, actual experimental estimation procedures 
for the spectral density involve a great many parameters. 
Indeed, its experimental estimation can be at least as 
involved as the estimation of 
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Appendix I 

We want to derive (23) in the text. Starting from (10) 
we have 


(<rl(N, T, r)> 

tfhS^-ht4) 

X (yl) - b X X <M,> 

n- 1 iV ,-i ,„} 

-w^{xjr dr C d, '< y(, ' Mn> 

t AT M /• I i t r /» t j + r | 

- X X J ( dt" J' dl' <y(02/(n>j (56) 

where (9) has been used. Now 

{y{t’)y{t")) = R v {t' - l") (57) 

where U„(t) is the autocorrelation function of y(t) and 
is the Fourier transform of <S V (/) , the power spectral 
density of y(t). Equation (57) is true provided that 
y(t) is stationary (at least in the wide or covariance 
sense), and that the average exists. If we assume the 
power spectral density of i/(t), S„(/) has low and high 
frequency cutoffs /, and j h (if necessary) so that 




f S,(f) df 

Jo 

exists, then if y is a random variable, the average does 
exist and we may safely assume stationarity. 

In practice, the high-frequency cutoff f h is always 
present either in the device being measured or in the 
measuring equipment itself. When the high-frequency 
cutoff is necessary for convergence of integrals of S v (f) 
(or is too low in frequency), the stability measure will 
depend on j h . The latter case can occur when the measur- 
ing equipment is too narrow-band. In fact, a useful 
type of spectral analysis may be done by varying /» 
purposefully [18]. 

The low-frequency cutoff ft may be taken to be much 
smaller than the reciprocal of the longest time of inter- 
est. The results of calculations as well as measurements 
will be meaningful if they are independent of fi as / ( 
approaches zero. The range of exponents in power law 
spectral densities for which this is true will be discussed 
and are given in Fig. 1. 

To continue, the derivation requires the Fourier trans- 
form relationships between the autocorrelation function 
and the power spectral density 

S„(f) = 4 [ K u (t) cos 2 t/t aV 
Jo 

R v (t) = ^ SM) COS 2jr/r df. ^ 

Using (58) and (57) in (56) gives 
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dt" 


■ f’ +r dt' cos 2x/(f' - t") if: r df SM) 

Jim iV i-1 j_| Jo 

■ f ' dt" f dt' cos 2ir/(r - *")} 

= (A T - ])r 2 {|l [/ 0 dj SM) ~W$~\ 

-htt[l' d <^ (2eos2rm, - i) 

- cos 2ir/[7'(i - i ) + r] - cos 2 t t/[T(; — i) — r])J j. 

(59) 

(The interchanges in order of integration are permissible 
here since the integrals are uniformly convergent with 
the given restrictions on S„(/).) The first summation in 
the curly brackets is independent of the summation in- 
dex n and thus gives just 




The second term is then 

-£(/.’ “i w sin ’ * £ § 005 2 ' TO - 4 (62 > 

(The interchange of summation and integration is justi- 
fied.) We must now do the double sum. Let 
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(60) 


The kernel in the second term in the curly brackets 
may be further simplified 

2 cos 2ir/7’(j — i) — cos 2r](T(j — i) + r) 


Fig. 5. Region of summation for * and k for N = 4. 
This may be written as 

S- W + 2B.[iF-Ii]!> (66) 

where Re [ C7] means the real part of V and d/dx 
is the differential operator. The series is a simple geo- 
metric series and may be summed easily, giving 

S -AT + 2Re{[j V -l£]^ 

[ l - e iNl - N(l - e ix ) 

\ 4 sin 2 x/2 


- i). 
(61) 

S = N + 2 Re 

= iV + 2 Re 


sin 2 Nx/2 

(62) 

sin 2 x/2 


(67) 


Combining everything we get, after some rearrangement, 
T, r)> 


j — i = k 

2wjT = x . (63) 

Changing summation indices from i and ; to i and k 
gives for the sum 

AT AT N N-i 

S = X3 X cos x U — i) - X X cos kx. (64) 

i-l j -I t-1 Jk-l-i 

The region of summation over the discrete variables i 
and k is shown in Fig. 5 for IV = 4. 

The summand is independent of i so that one may inter- 
change the order of summation and sum over i first. 
The summand is even in k and the contributions for 
k < 0 are equal to those for k > 0, and so we may pull 
out the term for k = 0 separately and write 

( ,V-1 N-k \ N 

X cos kx X 1 ) + X 1 
= 2(li (N - k) cos kx) + N. (65) 



sin 2 irriNr 
N 2 sin 2 7r r/r_ 


( 68 ) 


where r = Tfr. This is the result given in (23). 

We can determine a number of things very easily from 
this equation. First let us change variables. Let x/r = u, 
then 


T, r)) 

N f " „ ( u \ sin 2 u f _ sin 2 Nru \ 

~ (N - l)xr J 0 v Wr/ u 2 \ W 2 sin 2 ruf 

(69) 

The kernel behaves like u as u — * 0 and like w~ 2 as u — * 00 . 
Therefore {<r‘ v {N, T, r)) is convergent for power law 
spectral densities, S t (j) = hj a , without any low- or high- 
frequency cutoffs for —3 < a < 1. Using (69) for power 
law spectral densities we find 

T, t)) = T~°- l h a C a , -3 < a < 1 

= T*h a C a , n = -a - 1 
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and 


C„ = 


N 


(N - 1)tt 




du u a 


sin' u 
u 2 


/i_^L 

\ A' 2 si 


Ara_j 


sin ru 


(70) 


This is the basis for the plot in Fig. 1 in the text of u 
versus a. For a > 1 we must include the high-frequency 
cutoff / A . 

For N = 2 and r = 1 the results are particularly sim- 
ple. We have 


(<r 2 (2, t, r)) = r~ a 'h a [ du u° 2 sin 4 u (71) 
?r“ Jo 

for power law spectral densities. For N = 2 and gen- 
eral r we get 


<*:< 2 , t, r)> 



Let y(t) be a sample function of a random noise pro- 
cess with a spectral density S„(/). The function y(t) is 
assumed to be pure real and S„(/) is a one-sided spectral 
density relative to a cycle frequency (i.e., the dimensions 
of S v (f) are that of y 2 per hertz). (For additional infor- 
mation see Appendix I, [7], [8], [18].) 

Let x(t) be defined by the equation 


*«> - 1 - »«)■ 


(73) 


Define the following. t 0 is arbitrary instant of time and 
<„ + . = t n + T, n = 0, 1, 2, • • • , (74) 

x(t„ + t) - x(L) 


1 f 

tin = - / 2/(0 d < = 


1 * 

(^)a = it 

* * n “ 1 


(75) 

(76) 


and let be a high-frequency cutoff (infinitely sharp) 
with 2 it/at 3> 1. 


1 — cos 2u~ cos 2ru- 


cos 2u(r+ 1) 
2 


cos 2u(r— 1) 
2 


u 


Definition: 

(<rl(N, T, r)) = (~T[ g (P. - ^)v) 2 )- (77) 


7TT 

The first form in (72) is particularly simple and is also 
useful for r = 1 in place of (71). 

Let us discuss the case for «> 1 in a little more de- 
tail. As mentioned above we must include the high-fre- 
quency cutoff /a for convergence. The general behavior 
can be seen most easily from (68). After placing the 
factor t -2 outside the integral and combining the factor 
f~ 2 with S v (f) we find that the remaining part of the 
kernel consists of some constants and some oscillatory 
terms. If 2 u-/at 2> 1 it is apparent that the rapidly oscil- 
lating terms contribute very little to the integral. Most 
of the contribution comes from the integral over the 
constant term causing the major portion of the t de- 
pendence to be the r~ 2 factor outside the integral. This is 
the reason for the vertical slope at p = —2 in the 
fx versus « plot in Fig. 1 in the text. 

One other point deserves some mention. The constant 
term of the kernel discussed in the preceding paragraph 
is different for r = 1 from the value for r ¥= 1. This is 
readily seen from (72) for N = 2; for r = 1 the constant 
term is 3/2 while for r ^ 1 it is 1. This is the reason 
for 8 fc (r — 1), which appears in some of the results of 
Appendix II. In practice, 8 fc (r — 1) does not have zero 
width but is smeared out over a width of approximately 
(2 tt/at) - 1 . If there must be dead time r ^ 1, it is wise to 
choose (r — 1) (2 tt/at) _ 1 or (r — 1) <SC ( 2irf h T )~ 1 but 

with 27 t/at » 1. In the latter case, one may assume 
r ~ 1. 


S, 


(u\ sj 
\tt) 


sin 2 u sin 2 ru 


(72) 


(78) 


(79) 


Special Case: 

(*1(2, T, r)) = ( ( £z - ~ ^ )- 

Special Case: 

<rl(r) = (vy(2, t, t)) 

= ( [x(t n + 2 t) - 2x(t 0 + r) -f r(/ 0 )] 2 \ 

Definition: 

D\( r) m ([x(l 0 + 2 t) - 2x(/ 0 + r) + x(« 0 )] 2 ). (80) 
Consequence of Definitions: 

D\( r) = 2rV 2 (r) ^ 2ol(r). (81) 

Definition: 

il(T, r) - ([x(to +T + r) - x(t 0 + T) 

- x(t 0 + r) + x(<o)] 2 ). (82) 

Consequence of Definitions: 

il(T, r) = 2rV„(2, T, r)). 

Special Case: 

tl(r, r) = Dl( r). 

Random Walk y 

s.(0 = Y (s,(fl = 

r « - , 0 </</». 

r 


(83) 

(84) 
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Quantity 

Relation 


Wl( X, T, T)> 

h - 2 ' — 12^ r| [r(Ar + 1} ~ !!• r 

> 1 
(85) 

(°v(Af, t, t)) 

r = 1 

(86) 

<^(r) 

ft— , N = 2,r = l 

(87) 

#(r) = 2«r 2 (r) 

. 2(2jt) 2 (t| 3 

'•’-2 .1 

b 

(88) 

*l(T, r) 

fe-2- ( — (3r - i), r > 1 

■ (89) 

Flicker y 



SM - *f (s.U) - 



r = 7’/ r, 0 </</». 


Quantity 

Relation 


Wi(N, T, r)> 

h - r N(N~ l)£ (Ar n) 



■ (~2(nr) 2 ln (nr) + (nr + l) 2 ln (nr -f- 1) 


+ (nr — l) 2 ln |nr — 1 1) 

(90) 

(<rl(N , r, t)> 

t AT ln N , 

_ 1 • (r - L) 

(91) 

<^(r) 

A_,-2 ln 2, (N = 2,r = 1) 

(92) 

D, 2 (r) = 2o 2 (t) 

/i_, -4r 2 ln 2 

(93) 

tf(T, r) 

h. t ■ r 2 ( — 2r 2 ln r + (r + l) 2 ln (r + 1) 


+ (r — l) 2 In |r — 1|] 

(94) 


~li-i -2r 2 (2 + ln r), r » 1 

~A_, -2T 2 (2 - Inr), r«l. 

(95) 

White y ( Random Walk x) 


s.(f) = Ao 5,0) = ^5 



r = T/r, 0 </</». 


Quantity 

Relation 


(°l(N, T, r)) 

r > 1 



Ao-w + 1) Irl’ 1 , ATr < 1 

(96) 

(<rl(N, t, r)) 

r = 1 

(97) 

a*( T ) 

fH’-r. N = 2, r = 1 

(98) 

D](r) = 2 <r 2 (r) 

A 0 ' | r I 

(99) 


tl(T, r) io-M, r > 1 

h 0 T, T < I. (100) 

Flicker x 

S.U) = K \)\{s z (f) = 

r = T/r, 2irf h T » 1, 27T UT » 1, o < / < 1„. 

Quantity Relation 

<«;(», T, ,)) {2 + In (2*/,r) 

+ A-<AT=Tj § <" “ "> 

'"[y£b]}. '»> («•« 

(»:(Y, T. t)> 777(2^ [ 2 + ln ~ > 

r = 1 (102) 

<rl ( r) h, { 3[2 + In (2r/*r)J - In 2}, 

iV = 2, r = 1 (103) 

Dl(r) = 2<r 2 (r) + ln (2rf k r)] - In 2} (104) 

*l(T, r) hr~ji [2+ ln (2>0 *t)], r » 1 

V^{ 3[2 + ln(2rf»r)] - In 2}, 

r = 1 (105) 


hr -~ 5 (2 + In (2^7’)], r«l. 

White x 


S.(f) = Kf ( 5,0 ) = ^3) 


r = 

T/r; a*(r - 1) = j 1, lfr 1 
10, otherwise 



2irj K r » 1, 0 </</*. 


Quantity 

Relation 


(ol(N, T, r)> 

■ N + 5 t (r - 1) 2 f t 
N(2 tt) 2 t 2 

(106) 

( a \(N , r, r)> 

h . N + l M r 1 

.** N(2v} 2 r 2 ’ r “ 1 

(107) 

<rj(r) 

N = 2,r = 1 

(2ir) r 

(108) 

7>’(r) = 2<r 2 (r) 

, 6/* 
h2 :(2irf 

(109) 

**(T, r) 

hi- [ 2 + 5 ,(r - 1)1^5 

(2jt) 

(110) 
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